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INTRODUCTION 


HEN a person is immersed in a sound wave 
he experiences the sensation of hearing. 
Psychologists usually speak of this sensation as 
the sound. In this sense the sound has three 
general characteristics, namely, loudness, pitch 
and quality or timbre. It is only musical sounds 
which have a definite pitch. The term ‘‘quality”’ 
contains all the other auditory characteristics not 
included in the terms “loudness” and “pitch.” 
In recent years there have been attempts to 
separate quality into the other aspects of the 
sound such as volume, extension, brilliance, etc. 
It seems to me that quality and timbre are still 
very vague terms and in spite of the recent inter- 
esting work they are still far from being reduced 
toa quantitative basis. This is probably because 
there is such a great variety of things included in 
the term ‘‘quality.’’ On the other hand, pitch and 
loudness have been reduced to a definite quan- 
titative basis, although the former has a limited 
application, namely, to musical sounds, while the 
latter can be applied to every type of sound. 
I will discuss experimental and analytical 
methods used to reduce loudness toa quantitative 


basis and point out the bearing of the experi- 
mental results upon the mechanism of hearing. 
Before proceeding with this discussion I shall 
review some of the quantitative terms which will 
be used. 


DEFINITIONS 


In dealing in a quantitative way with the 
sound we will use the terms intensity J, acoustic 
pressure , and their relation to the intensity 
level 8; the loudness level L, and the loudness NV; 
the spectrum level B, and noise spectrogram ; the 
masking M, and noise audiogram. 


Intensity J 

As used in this paper, the intensity J will be 
used with 10-'* watts per square centimeter as a 
unit. In other words, J=1 for the standard 
reference level of 10-'® watts per square centi- 
meter. For example, if J is 10°, then when a free 
wave passes the point under consideration, 10~7 
watts of power, or 1 erg per second, pass through 
each square centimeter of the wave front. 


Pressure p 


As the wave passes, the acoustic pressure for 
unit intensity varies with temperature and 
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atmospheric pressure. These variations are very 
small compared to the range of values used here, 
so that for most practical calculations it can be 
assumed that when 7 is unity, then p can be 
taken equal to unity. This also corresponds to 
the standard threshold for a tone having a fre- 
quency of 1000 c.p.s. and corresponds approxi- 
mately to 0.0002 dyne per square centimeter. 
With these conventions then 


I=p’. (1) 
Intensity level 6 


The intensity level 8 in decibels is then defined 
by 


B=10 log J=20 log p. (2) 
Loudness level LZ 


For purposes of making comparisons of loud- 
ness of various sounds, a reference sound has been 
adopted by the American Standards Association 
as a reference standard. This same standard has 
recently been recommended for adoption to the 
International Standards Association by an inter- 
national acoustical conference. This reference 
sound is a free progressive wave having a fre- 
quency of 1000 cycles per second to be received 
by the listener facing the source of this tone. Its 
intensity could be defined by giving either J, p, 
or 8, but it has been agreed internationally that 
it shall be defined by giving the intensity level 
in decibels. The loudness level of the reference 
tone expressed in phons is the intensity level of 
this tone expressed in decibels. The loudness 
level of any other sound is determined by com- 
paring it to the reference tone, adjusting the 
latter until the two sound equally loud. Then the 
loudness level of the reference tone becomes also 
the loudness level of the sound being measured 
and is expressed in phons. 


Loudness N 


Loudness is the magnitude of the auditory sen- 
sation produced by the sound. A scale for repre- 
senting loudness has been found such that units 
on this scale agree with common experience in 
the estimates made of the sensation magnitudes. 
This scale is definitely related to the loudness 
level scale. Units on the scale are called loudness 
units, abbreviated LU. 
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Spectrum level B 


When a sound is not musical but is what is 
ordinarily called a noise then the vibration rates 
are scattered throughout the audible frequency 
range and it can be considered as a continuoys 
spectrum. An important noise having a cop. 
tinuous spectrum is the thermal noise which Was 
used in the experiments to be described later, If 
AT is the intensity or power per square centimeter 
carried by the frequency band between f and 
f+Af, then the spectrum level B at the frequency 
f is defined by 


B=10 log (AJ/Af). (3) 


The actual intensity level 8 of any band of fre. 
quencies between f; and f2 is given by 


f2 
p=10 log [ 107/19gf, (3A) 


fi 
Noise spectrogram 


The values of AJ for a given Af can be obtained 
experimentally by means of a noise analyzer, 
Consequently the values of B at each frequency 
can be obtained. A curve plotted with the fre- 
quency f as abscissae and B as ordinate is called 
the spectrogram. For purposes of loudness cal- 
culations it is not necessary to make Af smaller 
than 100 c.p.s. 


Overtone structure 


If the sound is a musical tone it can be defined 
by giving the intensity level of the fundamental 
and of each overtone. Strictly speaking, one 
should give also the phases between the com- 
ponents, but it has been found that for these 
steady sounds phases produce only a negligibly 
small effect on the auditory sensation and in this 
discussion phase changes are excluded. The over- 
tone structure is usually given in terms of the 
frequency and intensity levels of the component. 
Masking M 

Whenever one sound is being heard it reduces 
the ability of the listener to hear another sound. 
This phenomenon is called masking. In testing 
the masking properties of a sound usually pure 
tones are used as the masked sound. The number 
of decibels that the threshold level of a pure tone 
is shifted due to the presence of the noise is called 
the masking M at the frequency corresponding to 
that of the pure tone. 
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Noise audiogram 


The values of M for each frequency are deter- 
mined and the corresponding values of f and M 
are plotted to form a curve. The resulting curve 
is called the noise audiogram. 

So it is seen that a noise may be defined by 
giving either its overtone structure, the noise 
spectrogram or the noise audiogram. 


DEVELOPMENT OF A TRUE LOUDNEss SCALE 


A true loudness scale must be one such that 
when the number of units on this scale is doubled, 
the magnitude of the sensation as experienced by 
typical observers will be doubled, or when trebled 
the loudness sensation is trebled, etc. The exist- 
ence of such a scale is dependent upon the sup- 
position that consistent judgment tests of 
loudness can be obtained. Experimental tests 
indicate that such judgments are sufficiently 
consistent to warrant choosing such a scale. 
To find the relation between the loudness level 
and the loudness on the proposed new scale, 
various types of experiments were performed by 
investigators in our laboratory and elsewhere. 

The general scheme of these experiments is to 
fnd experimentally pairs of values of loudness 
level L; and Le phons, which correspond to pairs 
of values of the loudness NV, and Ne LU. The 
first type of experiment is concerned with a com- 
parison of the loudness experienced by a typical 
observer when listening first with one ear and 
then with both ears. Since we are dealing with 
the typical observer, the two ears will be alike in 
sensitivity. Consequently, if loudness is propor- 
tional to the total nerve energy sent to the brain, 
then it follows that the loudness experienced 
when using two ears will be double that when 
using one ear. To obtain such data the sound is 
frst produced at any convenient intensity and 
listened to with only one ear, and its loudness 
level measured. Let us call this value L;. Without 
changing the intensity of the sound the observer 
listens with both ears and the loudness level is 
measured. Let us call its value Lo. If N; is the 
number of loudness units corresponding to listen- 
ing with one ear and N, to that corresponding to 
listening with both ears, then it is evident that 
N; must be just double N;. We will again listen 
with one ear, but increase the intensity of the 
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‘sound until a loudness level of Lz is obtained. If 
now again the intensity of the sound is held 
constant and we listen with both ears, we will 


- find a loudness level Ly. This must correspond to 


a loudness N, which is four times N;. And so we 
might continue throughout the entire intensity 
range. If we now choose JN, arbitrarily equal to 
some convenient number, we will have a set of 
corresponding values of loudness units and loud- 
ness levels from which we can draw a curve 
giving the relation which we are seeking. 

The experiments as outlined could not be 
carried out because our typical observer was 
hypothetical. To obtain results corresponding to 
that which would be obtained by such a typical 
observer, one must deal with a large group of 
observers, taking the average value as that cor- 
responding to a typical observer. For this reason, 
pairs of loudness levels ZL; and Le were obtained 
starting with Z; at various values throughout the 
audible range. The points indicated by the open 
square in Fig. 1 are such observed points. A solid 
curve is drawn to give the best fit of the observed 
points corresponding to N;, 2Ni, 4N,; and 8N,. 
From this curve the proper values of Li, Le, Ls, 
Ls, etc., can be obtained to build the scale in the 
manner described above. 

In the second type of experiment the following 
procedure is used. A 1000-cycle reference tone is 
balanced against a pure tone having a frequency 
far removed from 1000 until they are equally 
loud. Then the two tones are combined to form 
a complex tone having two components. Let L; 
be the loudness level of each component alone 
and Lz the loudness level of the combination. 
Since the components were adjusted to produce 
equal loudness, it is assumed that each is pro- 
ducing the same nerve stimulation and that when 
they act together they will produce twice the 
loudness. This hypothesis is reasonable provided 
in the mechanism of hearing the two sets of 
nerves do not interfere with each other. It can 
be shown by experiments on the masking of one 
tone by another that such interference is negli- 
gibly small when their frequencies are widely 
separated. Consequently, balancing such a two- 
component tone against the reference tone should 
give pairs of values of ZL; and Lz which can be 
used in the same way as indicated by the one-vs.- 
two ear method described above, that is, L2 cor- 
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Fic. 1, Experimental values of two loudness levels Z; and Lz such that the sound 
corresponding to the latter is twice as loud as that corresponding to the former. 


responds to a loudness number which is just England have taken a very large number of 
twice that corresponding to L;. Such observed observations by this method, not only on the 
points are shown in Fig. 1. 

In the third type of experiment a typical ob- 
server is asked to listen to a sound, first at a 
loudness level L; and then to the same sound 
when its loudness level has been raised to a value 
Lz such that to him the loudness seems to have 
doubled. The experiment is repeated, starting 
with another initial value of L,, and thus a series 
of pairs of values of ZL; and L2 are obtained which 
have the same relation to our loudness scale as 
those obtained by the first two methods. A much 
greater amount of data is required by this 
method because the probable error of a single 
judgment is very large. However, this method is 
related more directly to the scale we are seeking LOUDNESS LEVEL,L,, IN PHONS 
than the two preceding ones. Four different in- Fic. 2. Relation between loudness level in phons and the 
vestigators! in different parts of this country and loudness in LU. 
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the 1000-cycle tone this was taken as the loudness level L. 
_1]In these experiments the loudness level was not meas- _ For other pure tones the loudness level was computed from 
ured but only the number of db above the threshold. For our data on the loudness level of pure tones. 
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Fic. 3. Experimental values of two loudness levels L; and Lio such that the sound 
corresponding to the latter is ten times louder than the former. 


1000-cycle reference tone but on many types of 
tones. The data taken by Laird have been 
omitted since they were so far out of line from 
those obtained by the other three groups of 
investigators, namely, Ham and Parkinson, 
Firestone and Geiger, and Churcher-King- 
Davies. Their observations are given in Fig. 1 by 
the points indicated. 

It will be seen that the observed points ob- 
tained by these diverse methods fit consistently 
on one curve, particularly well in the range of 
loudness usually experienced in everyday life; 
namely, from L2=40 to L:=80 phons. It was 
suficiently consistent to justify choosing the 
resulting loudness scale as a true loudness scale. 

Using this curve in Fig. 1, pairs of values of L; 
and L, were taken to construct the curve of Fig. 2 
in the following way. It will be seen that a stair- 
case has been erected upon the curve with the 
corners of the steps lying on a 45 degree straight 
line. The value of N for L =0 was taken as unity. 
If, then, the sound is adjusted so that the 
threshold is reached for either ear listening alone, 
then when listening with both the loudness level 
is between 2 and 3 phons and the loudness is 


doubled, or 2. This corresponds to the bottom 
step. Now if the sound is adjusted so that for 
each ear alone the loudness level is 2 phons, and 
then if both ears are used the loudness is raised 
to 4 LU. This corresponds to the second step. 
So it is seen that the width of the steps shows the 
change in the loudness level necessary to double 
the loudness. The values of the loudness N 
indicated on each step can now be plotted against 
the corresponding L to give the desired relation 
between L and N. This is shown by the curve of 
Fig. 2. 

This same scale can be developed from data 
obtained by balancing a three-component tone 
against the reference tone or by judging when a 
threefold increase in loudness has been produced. 
Or we could use any other number besides 2 or 3 
and get the same scale if it is a real loudness 
scale. An extreme case would be to use a ten- 
component tone and use judgment tests for 1/10 
and ten times as loud. Tests of this kind have 
been made which I shall now describe. 

A complex tone? was generated having ten 


? The description of this generator is given in the paper by 
Fletcher and Munson, J. Acous. Soc. Am. 5, 82-108 (1933). 
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harmonics of a fundamental frequency of 530 
c.p.s. The components were each adjusted so 
that when listened to separately they had the 
same loudness level L: phons. They were then 
combined into a complex tone, and it was found 
that the loudness level was raised to a value Lio 
phons. So in this way pairs of values of loudness 
levels L; and Li were obtained starting with L, 
at different intensities of sound. If the patches 
of nerve stimulation corresponding to each com- 
ponent do not overlap then the complex tone 
should be ten times as loud as each component. 
Such pairs of intensities were obtained using a 
group of observers, and are shown by the points 
in Fig. 3. Ham and Parkinson, and Firestone and 
Geiger made judgment tests to find a pair of 
loudness levels L; and Lio such that the tone for 
L,; was one-tenth the loudness of that for Ly». 
Their results are shown by the points A and V, 
respectively. The points are somewhat scattered 
but I think it is remarkable that the two sets of 
points determine the single curve within the 
observational error, for such judgment tests are 
extremely difficult to make. It is surprising that 
the agreement is as good as indicated. The judg- 
ment tests and the ten-tone balance tests were 
all made in different laboratories and before any 
loudness scale of the type discussed here was 
constructed. 

Now the curve of Fig. 3 can be used to con- 
struct a loudness scale. In this case the width of 
the steps of the staircase corresponds to the 
phons increase necessary to make the sound ten 
times as loud. The corresponding values of N and 
L can be used to determine the curve of Fig. 2. 
The solid curve was drawn so that points on it 
would yield the same loudness scale as shown in 
Fig. 2. It is seen that such a curve does indeed 


fit the observed data. The following equation 
105/37 105/31 Qz/10 
N=——_-______ 
(I+ 105/?)?/3 (104/19-+ 105/2)2/8 


proposed by Knauss represents approximately 
the relation shown in the curve of Fig. 2. This 
reduces to 

N=I=104/" (5) 


for values of L below 20, and 


N=105/3(I)"/3 = 105/21049= 46 X104/™ (6) 


for values of L above 40. This last relation fit, 
the curve accurately, the other two only 
approximately. 

It is important to remember that the values of 
the intensity J in the above equation refer only 
to the 1000-cycle reference tone but the values 
of L can be for any type of sound. The curve of 
Fig. 2 has been taken as defining the relation 
between the loudness N expressed in LU’s ang 
the loudness level L expressed in phons. The 
loudness of any other tone, either simple or 
complex, can be determined by comparing it to 
this reference tone, and adjusting the intensity of 
the latter until the two sound equally loud. The 
loudness of the reference tone thus adjusted js 
then the loudness of the unknown sound being 
measured. 


THE DETERMINATION OF THE POSITION Coop. 
DINATE x OF THE AUDITORY NERVE Patcy 
Wuicu Is Most EXciTeD BY A TONE 

OF FREQUENCY f 


Besides direct anatomical measurements there 
are two methods of determining this relation, 
The first of these was proposed by Wegel and 
Lane in 1923. It is based upon measurements of 
the minimum perceptible increments of pitch. 
A second one is given in this paper and is based 
upon masking measurements of thermal noise. It 
is generally agreed that the auditory nerve ter- 
minals are spread out along a winding spiral line 
so that a single coordinate x will define any 
particular group. This is usually considered to 
be the array of nerve endings along the basilar 
membrane. This membrane is about 30 mm long 
and 1/2 mm wide. It is not necessary to know the 
actual distances or the distribution of the nerve 
terminals if we define x as the percent of the 
total nerve endings that the maximum stimula- 
tion passes over as the frequency of a stimulating 
tone passes from the lowest audible pitch to the 
pitch corresponding to a frequency f. The posi- 
tion coordinate x thus defined is very convenient 
and will be used for loudness calculations. The 
position coordinate then must remain between 
0 and 100. The problem then is to find the rela- 
tion between f and x. Let this relation be 
expressed by 


x=(f). (7) 
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Fic. 4. Relation between the position of maximum nerve 
stimulation and the frequency of the stimulating tone. 
Data obtained from minimum perceptible differences in 
pitch. 


The basis of the first method rests upon the 
assumption that for each minimum perceptible 
step on the pitch scale the position of maximum 
stimulation changes over the same number of 
nerve terminals. 

Putting this assumption into an equation 


Ax = Ab = (dx/df)Af=constant (8) 


where Ax is the percent of the total nerve endings 
passed over in a just perceptible change in fre- 
quency Af. Since the Ax is assumed to be always 
the same, we have the equation 


dx 1 
—=—A (9) 
df Af 

where A is a constant. Integrating and deter- 
mining the constant by the condition that x= 100 
when the frequency goes from the lower limit / 
to the upper limit u, we have 


x = 100-——_-. (10) 


anal 
- 


The values of Af at different frequencies are 
determined from experiments. 

There are difficulties in applying this method 
as these values of Af depend upon the intensity 
level of the tone as well as the frequency. The 
relation @(f) then also depends upon whether 
you go through the pitch range on a loudness 
level contour, a sensation level contour or an in- 
tensity level contour. It also depends upon what 
limits you choose for u and /. Also there is consid- 


erable uncertainty’ that pitch is determined 
solely by a change in the position of maximum 
stimulation for the low pitched tones. This 
problem was discussed by Steinberg‘ in a paper 
presented before the Acoustical Society. Not- 
withstanding the difficulties mentioned above, 
the curve for #(f) which he gave as being the 


“most probable is given in Fig. 4. He showed that 


this agreed very well with results from animal 
experimentation. 

In the second method results are obtained 
from masking experiments. Electrical thermal 
noise is generated in an electrical conductor due 
to the haphazard movement of the electrons. 
This electrical noise can be amplified by vacuum 
tubes and converted into a sound by means of a 
telephone receiver or loudspeaker. Such an 
electrical disturbance has the property that the 
electrical power in any frequency band is propor- 
tional to the width of that band. This makes it 
possible to use thermal noise for building up a 
continuous spectrum of any desired form. The 
electrical noise going from a resistance to the 
amplifier has the same intensity for equal band 
widths regardless of where the bands are located 
on the frequency scale. However, by choosing the 
transmission characteristics of the amplifier and 
telephone receiver for converting the electrical 
noise to acoustic noise, one can obtain any desired 
type of sound spectrum from such thermal noise. 

Let us assume that a spectrum of thermal 
noise is impressed upon the ear, defined by its 
spectrogram; that is, by giving B for each fre- 
quency throughout the audible range. Then we 
wish to find the stimulation of the nerve patches 
at the position x. The intensity AJ contained in 
the small frequency band between f and f+Af 
goes mostly to stimulate the patches between x 
and x+Ax. There is of course some power going 
to patches outside this region. Since there are 
no abrupt changes in the spectrum, the patches 
outside this small band at x will already be 
stimulated by power from other frequency bands 
outside the small band at f. So it seems reason- 
able to consider that the major part of the 
stimulation in the element between x and x+Ax 
is from the frequency region between f and f+Af, 


3 The solution of the problem when Ax is not constant 
will be dealt with in another paper. 
4 J. C. Steinberg, J. Acous. Soc. Am. 8, 176 (1937). 
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Fic. 5. Spectrogram of noise for producing constant stimu- 
lation of the nerve endings at all positions. 


and that all the loudness from the patches 
located between x and x+Ax will be due to the 
intensity coming from the region between f and 
f+Af. When the threshold level is reached the 
patch of nerves is excited at the position x where 
the maximum excitation occurs. The acoustic 
intensity AJ due to the thermal noise between f 
and f+Af is then 108°Af. The intensity, when 
confined to a narrow band of frequencies that is 
necessary to excite the threshold of hearing, is 
100/10, where 8» is the intensity level theshold 
for a pure tone of frequency f. This assumes that 
the intensity is the same whether it is concen- 
trated in a single frequency or spread throughout 
the band between f and f+Af. This assumption 
has been justified by experimental determination 
of the threshold for small bands of thermal noise. 

For convenience let us call one percent of the 
nerve endings a patch. Then the values of x cor- 
responding to the whole numbers from 0 to 100 
locate each of these patches. Let the power of 
agitation of a patch be J, and let Jp be the power 
of agitation of the patch for it alone to produce 
the threshold of hearing. It will be assumed that 
this power can be spread over all the patch or 
concentrated in a small part of it and still have 
the same value for the threshold. This is true 
because at the threshold the loudness is numeri- 
cally equal to the intensity. For example, if the 
intensity is spread over twice as many patches 
the intensity per patch is halved, and conse- 
quently the loudness per patch is halved, but the 
total loudness remains the same because there 
are twice as many patches excited. Then the 
following proportion will hold 
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Fic. 6. Difference between the intensity level of the 
masked tone at each frequency f and the spectrum leve| 
of the thermal noise at the same frequency. 


J dx 108/%df dx Jo 
=— or —=[ — }104—Ao)/10 (11) 
Jo 1040/10 df \J,z ; 


This would be sufficient to determine the rela- 
tion between f and x if we knew how J,/J, 
varied with frequency ; or conversely if we knew 
dx/df we could determine the variation of J/J, 
for any spectrum B since 8» is known. Whenever 
the nerve patch is agitated a tone which js 
normally sensed by this patch will be masked. 
It seems reasonable to assume that the amount 
of masking produced can be taken as a measure 
of its agitation, or at least to assume that when 
the masking is uniform for all frequencies all of 
the patches have the same relative agitation, 
that is, the same value of J,/J. The spectrogram 
for a thermal noise which will give such a uniform 
masking has been experimentally determined. It 
is given in Fig. 5. It was shown in the paper by 
Fletcher and Munson® that the shape of such a 
spectrum for producing uniform masking is 
independent of the total noise level used. Since 
for such a spectrum J,/J) is a constant inde- 
pendent of the frequency, we have 


Jo\ 
(>) f 10 (8-80) /0g fF, (12) 
J: = 


The values of the functions B and £» at various 
frequencies can be taken from experimentally 
determined curves. The quantity J,/Jo can be 
determined from the requirements that x must 
be 100 when f= ~, which gives 


, 1 ng 
(=) = 10@-s01gf (13) 
Jo) 100 





°H. Fletcher and W. A. ie J. Acous. Soc. Am. 9, ! 


(1937). 
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Fic. 7. Relation between the position of maximum nerve 
stimulation and the frequency of the stimulating tone. 
Data obtained from masking effects of thermal noise. 


This equation then enables one to calculate the 
relative power of agitation per patch in terms 
of the noise spectrogram. Combining (12) and 
(13) we have 


f 
f 10 (B—60)/10g f 
a’ (14) 


f 10 (8-80)/10g F 
0 


It will be seen that adding a constant to B which 
is independent of frequency will not change the 
values of x. Therefore any one of the masking 
contours shown in Fig. 15 of the paper just cited 
can be used for the calculation. 
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Let the difference between the intensity level 
Bm of the masked tone and the spectrum level B 
of the thermal noise doing the masking be 
designated K or 

Bu—-B=K. (15) 


This difference is constant at any one frequency 
and independent of B except at very low values 
of B. Now the masking M is given by 


M=Bn—Bo. (16) 
Substituting and rearranging we have 
B-—Bo=—-K+M and M=B-—£,+K. (17) 


If this is substituted in (14) we obtain, since M 
is a constant independent of frequency, 


f 
J 10-* gf 


a (18) 


f 10-K/ogf 


This is the equation which is most suitable for 
obtaining the relation between x and f for it does 
not contain threshold values 8». The experi- 
mental values obtained for K are given in Fig. 6. 

It is important to notice that the values of K 
can be obtained without knowing the amplitude 
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Fic. 8. Chart illustrating how the relation between position coordinate x and the 


frequency f were obtained from experiments on masking of thermal noise. 
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frequency characteristics of the amplifiers and 
head phones used for transferring the electrical 
noise to acoustic noise. In a noise analyzer the 
band width is kept constant as the noise spec- 
trum is heterodyned through the audible fre- 
quency region. With such a device for producing 
the noise, it is only necessary to measure the 
ratio 7 of the currents or voltages at the telephone 
receiver terminals for the thermal noise and for 
the tone being masked. It can easily be shown 
that Eq. (18) can then be written 


f ‘af 


x= 100. (18A) 


J “Pd Lf 


The results of integrating Eq. (18) or (18A) by 
graphical methods are shown in Fig. 7. The curve 
from Fig. 4 is shown by the dotted line. There is 
a fairly good agreement. If pitch changes for low 
frequencies are not sensed entirely by displace- 
ment of maximum stimulation, then the two 
curves should be slightly separated and in the 
direction shown. However, the general agreement 
of these curves with anatomical data gives us 
considerable confidence in their validity. 

If S is the stimulation per patch expressed in 
decibels above its threshold, then by Eq. (13) 


Js ae 4 
S=10 log —= 10( tog —+— 
Jo 100 10 


+log f 10-Kiaf), (19) 
0 


The last term was evaluated graphically using 
the curve of Fig. 5 and was found to be 1.95. 


























TABLE I. 
Lo =86 db ] Lo =56 db 
L Obs.R. | Catc.R. ||' zc | Obs.R. | Calc. R. 

86 1.00 1.00 | 56 ~|~ 1.00 1.00 
82 0.68 0.75 | 54 0.93 0.85 
80 0.53 0.59 52 0.75 0.71 
76 0.41 0.45 49 0.55 0.58 
71 0.26 0.32 | 46 0.41 0.47 
66 0.20 0.21 41 0.25 0.32 
56 0.13 0.10 36 0.11 0.18 
46 0.08 0.05 | . 
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Consequently the above equation becomes 
S=M-0.5. | (20) 


This is an important relation because it states tha 
the stimulation per patch differs by only 0.5 @ 
from the masking of a tone having a frequency f 
corresponding to the patch position x. This is not 
exactly true for very low values of M below 15 
db or for very high values of M above 70 db 
because, although the shape of the curve giving 
K remains constant, its whole position is shifted 
downward so that S becomes 3 or 4 db less than 
M at these extreme limits. 

This same result can be obtained without using 
a spectrum producing a uniform masking. Let us 
assume that the stimulation S and the masking 
M differ only by a constant C independent of 
frequency. It is not even necessary to assume 
that C does not vary with different levels of ¥. 
Substituting these values, namely, 


M-S=C and M=B-£8,4+K, 
into Eq. (13) it becomes 
dx/df = 10 (8—-80—S)/10 — 4Q-K+C_ (18’) 


If this is integrated we arrive at Eq. (18) with 
value of C=0.5 db. So it is seen that either 
method leads to a method of evaluating the 
function ®(f) from measurements of the differ- 
ence K between the intensity level 6,, of a masked 
tone in any frequency region and the spectrum 
level B of the thermal noise in the same frequency 
region. 

The relations involved in this analysis are 
illustrated by the chart of Fig. 8. The frequency 
f and the position coordinate x are plotted 
linearly along the two lower horizontal lines. The 
position coordinate line is divided into 20 equal 
parts so that the positions correspond to 0, 5, 
10, 15, 20, etc. These positions of x are connected 
to the corresponding frequencies —0, 250, 390, 
520, 660, etc. Let us assume that the spectrum 
level B is the same for all frequency regions. 
Then the amount of power concentrated upon 
any 5 percent nerve patch, compared to the 
power concentrated upon the same patch by a 
single frequency tone which is just masked by 
the noise, is proportional to the distance between 
the corresponding frequencies. So it is seen that 
at high frequencies a much larger amount of 
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Fic. 9. Relation between the loudness per patch of nerve 
endings and the stimulation per patch in decibels. 


relative power is concentrated upon a 5 percent 
patch than at low frequencies. Consequently, a 
greater intensity level 8,, for the masked tone is 
required, as is shown by the curve at the top of 
Fig. 8. This curve is the same as that in Fig. 6 
except that the frequencies are plotted linearly 
instead of logarithmically. Since for this case B 
was assumed to be constant, the ordinates of this 
curve show the relative levels for the masked 
tone. 

One might interpret the curve of Fig. 8 in 
another way by assuming the noise spectrum 
to be such that a constant intensity level 8,, for 
the masked tone is produced. Then the relative 
acoustic power in the masked tone and the 
thermal noise concentrated upon each 5 percent 
patch must be the same for all patches. Conse- 


quently the quantity 10-2/!°Af must be constant 
for all frequency bands. Where the band Af is 
large therefore the value of 10-2"° must be 
small. In other words, if the ordinates of Fig. 8 
are considered as negative values, then they give 
the spectrum levels which will just mask a 
constant intensity level 8, of the masked tone. 
For such a condition the acoustic power in each 
of the bands shown is the same. For a large Af 
the power per cycle must be small, as is seen to 
be the case. 

Now that we know how to calculate the posi- 
tion coordinate x from f, and the loudness N from 
the loudness level L, we are in a position to 
describe several methods of calculating the 
loudness of a tone or a noise both when listened 
to by a normal ear and also by a deafened ear. 


COMPARISON OF CALCULATED AND OBSERVED 
JUDGMENTS OF FRACTIONAL LOUDNESS 


The loudness function will be first applied to 
the calculation of fractional reductions in loud- 
ness. If Lo is the initial loudness level and L the 
loudness level for reducing the loudness to the 
fractional amount R of the original loudness, then 


R=N/No=1062- 20/30 (21) 


for values of ZL greater than 40 db. For lower 
levels the N/ No can be calculated directly from 
Eq. (6). For example, here is a sample of such 
data taken by Ham and Parkinson.® If the tests 
are carried out with a 1000-cycle tone, the loud- 
ness levels L and Ly become simply the decibels 
above the standard threshold pressure 0.0002 
dyne per square centimeter. The table shows 
results upon this reference tone. (Table I.) 

The observed values were obtained by asking 
a large number of observers to set the attenuator 
so as to produce the tone at levels which seem to 
them to be reduced by the fractional amounts 
given. The fact that sets of data of this type do 


6 Ham and Parkinson, J. Acous. Soc. Am. 3, 511 (1932). 
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Fic. 10. Room noise spectrogram. 


agree with this loudness scale gives us consider- 
able confidence that this is a true loudness scale. 
Considerable data of this sort have been taken 
and the calculated values agree very well with 
those which are measured. 


CALCULATION OF THE LOUDNESS FROM THE 
NoIsE SPECTRUM 


When a patch containing 1 percent of the 
total nerve endings in the organ of hearing is 
stimulated S db above its threshold value it sends 
out a certain loudness to the brain which we will 
call N,. The value of S except’ at very low or 
very high values can be obtained by combining 
(17) and (20), or 


S=B-—6)+K-—-0.5. (22) 


The value of NV, should be dependent only upon 
S and should be independent of the position of 
this 1 percent of nerve fibers. In other words, it 
should be independent of the frequency. Con- 
sequently 

N.=Q(S) (23) 


where Q is some function which we can determine 
from our experiments on masking. The total 
loudness N from 100 percent of the nerve endings 
is then given by the equation 


100 
N= ; N.dx. (24) 


To determine the functional relationship indi- 
cated by (23) a noise spectrum of the shape 
shown in Fig. 5 is chosen. If such a spectrum is 
impressed upon the ear it will produce a constant 


7A more involved method of determining S at all 
values is given in the paper by Fletcher and Munson al- 
ready cited. 
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Fic. 11. Stimulation spectrum produced by room noise. 


stimulation per patch for all values of x. Con. 
sequently, since the nerve endings at various 
positions are stimulated alike, the loudness N, 
from each element is independent of the value of 
x. For such a spectrum we have the following 
relation 








100 
NEN, f dx=100N,. (25) 
0 


For various values of stimulation per unit 
length S the corresponding values of N were 
determined experimentally as outlined above. 
From these values of N the value of N, was 
determined from (25) for each value of S. Then 
the corresponding values of NV, and S are plotted 
in a curve to give the desired relation. It is shown 
in Fig. 9. 

As seen from Egg. (20) and (17), S—B is equal 
to K—6 )—0.5 and consequently is always the 
same function of the frequency for all values of 
B and therefore can be given in the form of a 
curve, or in a table. For convenience the cor- 
responding values of f, x, Bo, K and S—B are 
given in Table II. 

By means of this table, the value of S at any 
of the frequencies given is obtained by simply 
adding the values of B from the spectrum of the 
sound to the numbers given in the fifth line. The 
procedure, then, for calculating the loudness V 
from Eq. (24) is to find from the noise spec- 
trogram for each frequency f the corresponding 
value of B. From Table II find the corresponding 
stimulation level S. With this value enter the 
abscissae of Fig. 9 to find the corresponding N:. 
From Fig. 7, or from Table II, the value of x cor- 
responding to f is found. The point (x, Nz) is 
then plotted. This is done with enough fre- 
quencies to determine a curve going from x=0 
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LOUDNESS, 


to x=100. The area under this curve gives the 
yalue of loudness N which is sought. 

The application of this formula is limited to 
such noises as have essentially continuous spec- 
tra. If we use the masking ordinate determined 
experimentally, then it is shown below that a 
formula can be obtained which is applicable to 
ll classes of sounds which we have tested, both 
those having a continuous spectrum and also 
those having discrete frequencies. 


CALCULATION OF THE LOUDNEss OF RooM 
NOISE FROM Its SPECTRUM 


In trying to simulate conditions under which 
the telephone is used an artificial room noise has 
been created which has been used in some experi- 
ments at Bell Telephone Laboratories. The 
spectrogram of this noise is shown in Fig. 10. 
In Fig. 11 the corresponding stimulation spec- 
tum is shown with S as ordinates and x as 
abscissae. This is interesting because it shows 
directly what part of the nerve endings is 
stimulated most by such a room noise. By using 
Eq. (24) and the curve in Fig. 9, the loudness 
level for this noise spectrum was calculated to 
be 37 db. 

If the spectrum shown in Fig. 10 is amplified 
or attenuated, then similar calculations can be 
made from the resulting spectrum. In Fig. 12 are 
shown the values of the calculated loudness level 
for such spectra. The abscissae give spectrum 
levels at the frequency 1000 c.p.s. The crosses on 
the curve give the observed results as obtained 
by our laboratory technique. Here again it is 
seen that the agreement is well within the ob- 
servational error. 
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Fig. 12, Comparison of calculated and observed loudness 
levels for room noise. 
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CALCULATION OF LOUDNESS FROM THE 
MASKING AUDIOGRAM 


It has been shown that for every value of the 
stimulation per unit length S there is a cor- 
responding value of M. Consequently the func- 
tion NV, can be written as a function of (M) or 


N.=F(M). (26) 


The function F(M) can be obtained directly from 
the corresponding values of N,, and M from the 
same loudness data from which Q(.S) was deter- 
mined. The relation between F (M ) and M deter- 
mined in this way is shown in Fig. 13. Since S 
and M are approximately equal above 20 db 
these functions should be approximately equal 
above this value and are found to be so. Conse- 
quently, Eq. (25) may be expressed in the fol- 
lowing form 


N= f F(M)dx. (27) 


By means of this formula the loudness can be 
calculated directly from the masking audiogram 
determined experimentally. It can be applied to 
any type of sound in which the masking audio- 
gram is known. It should be applicable to pure 



































aay gman ae ol 
oseneeueen 
2 ==SSS== 

Fa 8 I ik a ow 
a a a a 





10 





LOUDNESS PER PATCH, F (mM) 








poi. - 


Fic. 13. The relation between the loudness per patch and 
the masking. 
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Fic. 14. Comparison between the observed loudness level of 
pure tones and that calculated from masking data. 


tones, or any combination of them. It should be 
emphasized that Eq. (27) is the only loudness 
formula which can be universally applied. These 
statements should probably be somewhat modi- 
fied to cover only those sounds which are con- 
tinuous a sufficient time to obtain an audiogram. 
When the sound is varying, the loudness is 
varying. There is no doubt that Eq. (27) still 
applies but the problem of determining the 
average loudness to correspond to judgment tests 
is a difficult one. When applied to a discrete fre- 
quency spectrum it becomes a tedious job to 
obtain an accurate audiogram. 


CALCULATION OF LOUDNESS OF PURE TONES 
FROM THEIR MASKING SPECTRUM 


A severe test of the general validity of Eq. 
(27) is the application to masking audiograms 
of pure tones. The function F(M) was deter- 
mined by a thermal noise having a continuous 
spectrum. If it is applicable to a tone having a 
single frequency, then it seems reasonable to 
suppose that it will be applicable to a large class 
of sounds. 
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Calculations’ from the masking audiograms of 
several pure tones were made. Typical results are 
shown in Fig. 14. These results are for pure tones 
having frequencies 200, 1000 and 4000 cycles per 
second. The dotted curves are the ones obtaingj 
from experiment and the crosses give the cal. 
culated points obtained from the masking audio. 
gram shown on the left-hand side of the chart 
It is seen that the agreement is very good. |f 
these audiograms are plotted with x as abscissa. 
and S as ordinates, they will then show directly 
the stimulation per patch at the various x pog. 
tions. If this is done it will be seen that the 209. 
cycle tone covers a much greater percent of the 
total nerve endings than the 4000-cycle tone, 


CALCULATION OF THE LOUDNESS OF CompLgy 
TONES FROM THE LOUDNESS OF EACH 
COMPONENT 


When the sound is composed of definite 
discrete frequency components a different ap. 
proach to the problem resulted in finding a 
workable formula which is better adapted for 
calculating the loudness of such tones. The 
loudness values for pure tones have been deter- 
mined experimentally for frequencies and inten- 
sities over the entire audible range. The loudness 
of any combination of these pure tones is cal- 
culated from the loudness of its components. 

When tones which are sufficiently separated 
in frequency are impressed upon the ear it has 
been seen that the number corresponding to the 
loudness experienced is obtained by adding the 
loudness numbers of each component tone. The 
effect of frequency separation is shown in Fig. 15. 
This shows the results of tests* on the combina- 
tion of ten tones adjusted so that each tone 
separately has the same loudness as every other 
component tone. It is seen that when the 
separation is 500 cycles, the combined loudness 
is always ten times that of a single component 


except perhaps for the very high intensities. As _ 


the frequency of separation becomes smaller, 
the intensity where this relation holds also 
becomes smaller. When the separation becomes 
as low as 50 cycles this relation fails to hold even 


for faint sounds. For a 50-cycle separation, the 


8 Experimental apparatus for producing such combina- 
tions is described in the paper by Fletcher and Munson, 
reference 2. 


TOTAL LOVONESS 'N LU 


OLTams of 
esults are 
ure tones 
cycles per 
Obtained 
> the cal. 
Ng audio. 
he chart. 
| good. If 
abscissae 
v directly 
US X posi- 
t the 200. 
ont of the 
e tone, 


COMPLEX 
“ACH 


~ definite 
“rent ap- 
finding a 
upted for 
nes. The 
en deter- 
nd inten- 
loudness 
es is cal- 
nents. 

separated 
‘ar it has 
ng to the 
lding the 
one. The 
n Fig. 15. 
combina- 
ach tone 
ery other 
then the 
loudness 
mponent 


sities. As ; 


smaller, 
Ids also 
becomes 
old even 
tion, the 


. combina- 
1 Munson, 


LOUDNESS, MASKING AND 
ten tones each having a loudness of 10,000 have 
a combined loudness of only a little more than 
twice that of a single component. If ten equally 
loud components are spaced equally between 
2,000 and 2,500 they will also be approximately 
equally intense tones. Consequently, from a 
loudness standpoint, the contribution is equiva- 
lent to a single tone of 2250 cycles per second 
which has been increased tenfold in intensity. 
A tenfold increase in the high intensity range 
corresponds to (10)!, or 2.15 times, which agrees 
with the loudness observations. In other words, 
for this case each component gave only 0.215, or 
about one-fifth the loudness that it would have 
given if allowed to operate on the nerve endings 
without the interference of the adjacent tones. 

In the paper cited above it was shown how to 
calculate a fraction called 6, which was always 
between zero and unity and which gave the 
fractional reduction in the loudness contribution 
for each component due to the interference to its 
neighboring components. In other words, if N is 
the loudness number corresponding to the 
complex tone and JN; the loudness number for 


its k’'th component if acting separately, then 


N= ON;. (28) 


The formula for calculating } is complicated 
and the reader is referred to the original paper 
cited above. It is sufficient here to state that 0 is 
unity for frequency separations greater than 500 
cycles and that all components within a band 
width less than 100 cycles can be lumped together 
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Fic. 15. The relation between the loudness of a complex 
tone having ten equally loud components, the loudness of 
each component, and the frequency separation between 
the components. 
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and treated as a single component. The values of 
N;, that is, the loudness of pure single frequency 
tones, have been experimentally determined for 
the frequencies and intensities throughout the 
audible range. The curves in Fig. 16 give a sum- 
mary of the results obtained. The numbers on 
the curves give the frequencies of the pure tones. 
For example, for all pure tones having frequencies 
between 800 and 2000 cycles the same relation 
holds between the intensity and the loudness as 
that given for the reference tone. It will be 
noticed that for the lowest pitched tone, namely 
the 50-cycle tone, the intensity variation at the 
ear which is required in order to reach the 
threshold of hearing is 500,000 times that re- 
quired for the reference 1000-cycle tone. It is 
important, however, to notice how rapidly this 
tone increases in loudness as its intensity in- 
creases. These low frequency tones stimulate 
practically the entire basilar membrane, even 
for only small intensities above the threshold, as 
contrasted with the 1000-cycle tone which 
stimulates only a very small fraction of the 
basilar membrane for low intensities. At very 
high intensities both of these tones stimulate 
almost the entire basilar membrane. 

The following example illustrates the use of 
these curves for complex tones which have com- 
ponents widely separated in frequency. 

Consider a tone having a fundamental fre- 
quency of 500 cycles per second and four har- 
monics, all equally intense, having an intensity 
of 10,000, that is, having an intensity level of 
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Fic. 16. The relation between the intensity level and the 
loudness of pure tones. 
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Fic. 17. Comparison between the calculated and ob- 
served loudness levels of pure tones for an ear having 
obstructive deafness. 


40 db. The loudness of the components can be 
determined from the curves of Fig. 16 and are 
found to be equal to 700, 1000, 1000, 1000, 1000. 
The combined loudness then is 4700. This loud- 
ness corresponds to a loudness level of the 
reference tone of 62 phons. For a more detailed 
application of Eq. (28) the reader is referred to 
the paper already cited. 

When the components are so close together 
that the sound may be considered as a continuous 
spectrum the above method becomes very com- 
plicated to apply and usually gives calculated 
results which are too high. It has been shown 
that such noises are better treated by loudness 
formulae using the spectrum and Eq. (24). 


CALCULATION OF THE LOUDNESS OF A SOUND 
FOR A PARTIALLY DEAFENED EAR 


First let us consider a case of pure conductive 
deafness—one in which the nerve terminals 
function normally. Let A be the function of fre- 
quency defining the audiogram of the deafened 
ear. Then if M is the masking for a normal ear 
(that is, an ordinate in the noise audiogram) the 
loudness experienced by the deafened ear is 
given by 


100 
n={ F(M —A)dx. (29) 

0 
If, however, the case is one of nerve deafness, 


then the loudness for the deafened ear should be 
given by 
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Fic. 18. Comparison between the calculated and ob- 
served loudness levels of pure tones for an ear having 
nerve deafness. 


N= [ (F(M) — F(A))dx. (30) 
In case of mixed deafness 
N= | " (F(M—Ac)—F(Ay))dx. (31) 


In a recent paper® Steinberg and Gardner dis- 
cussed these relations and found indeed that 
these equations gave calculated results which 
agreed with experimental observations of par- 
tially deafened persons. 

In the tests, observers having one deafened 
and one normal ear, listened to a tone first with 
one and then with the other ear. The intensity 
level of the tone on the normal ear was adjusted 
until it sounded as loud as a fixed level of tone 
on the deafened ear. By choosing different values 
of fixed level, the hearing of the deafened ear 
throughout the intensity range above 8, was 
studied. 

Figure 17 shows a typical case of conductive 
deafness where the deafened ear hears with 
reduced loudness at all intensity levels. The 
ordinates used in these curves and those of Figs. 
18 and 19 are levels above threshold. Since here 


9 J. C. Steinberg and M. B. Gardner, J. Acous. Soc. Am. 
9, 11 (1937). 
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we are dealing with pure tones, the relation 
between loudness level and level above threshold 
can be obtained from the curves of Fig. 16, the 
zero loudness levels being interpreted as the 
threshold levels. The upper charts show the left 
and the right ear audiograms. The points in the 
lower charts show the levels on the deafened and 
normal ears for equal loudness. If both ears 
were normal, the points would fall on the 45° 
dashed line through the origin. The displacement 
of the points to the right of the dashed line gives 
the hearing loss for the deafened ear which is 
seen to be the same at all levels. The solid line 
was calculated from Eq. (29). With this type 
of deafness the loudness perceived is the same as 
that perceived by a normal ear with the sound 
attenuated in db the same as the displacement 
of these straight lines, or equal to A, the ordinate 
of the audiogram at 250 c.p.s. 

Figure 18 is typical of a case of nerve deafness 
where the deafened ear tends to hear sounds 
above the deafened threshold with normal loud- 
ness. In this figure, the audiograms are shown 
by the left-hand charts. In both cases, a threshold 
loss of 60 db is shown for a 2000-cycle tone. 
However, when the level is raised 30 db above 
the deafened threshold, this tone is heard with 
practically normal loudness. The figures at the 


100 


250~ MASKING 








LEVEL ABOVE NORMAL fFHRESHOLD ON UNMASKED EAR 





right show the experimental data by circles and 
crosses. The solid curves were calculated from 
Eq. (30), the values of A being taken from the 
audiogram of Fig. 18 and the values of M from 
the masking curve of a 2000-cycle tone. The 
values of NV obtained were converted to loudness 
level L, and then knowing the frequency of the 
tone, to the level above threshold. The agreement 
of the solid lines with the observed points is a 
striking confirmation of the fundamental theory 
underlying all these calculations. 

Cases of deafness were also observed in which 
both types occurred in various proportions. The 
curves of levels for equal loudness on such 
deafened and normal ears were combinations of 
the curves shown in Figs. 17 and 18. The oc- 
currence of these two types of deafness can be 
explained from our picture of loudness previously 
described. 

It is of interest to note that one would expect 
the hearing loss caused by a masking sound to be 
of the type resulting from nerve atrophy, if the 
nerve fibers which are activated by the masking 
sound are ineffective in contributing to the 
loudness of a tone heard in the presence of the 
masking sound. To test this, one ear of an ob- 
server was masked by a wide band of thermal 
noise, extending over most of the audible range. 
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Fic. 19. Comparison between calculated and observed loudness levels of pure tones for an ear 
immersed in thermal noise. 
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Fic. 20. Relation between the position coordinate x and the 
distance x’ in millimeters along the basilar membrane. 


A pure tone heard in the unmasked ear was 
adjusted to sound as loud as a tone of the same 
frequency heard in the masked ear. The results 
are shown in Fig. 19, where it will be seen that 
the tones heard in the masked ear increase in 
loudness very much like tones heard in the nerve 
deafened ear of Fig. 18. The solid lines showing 
the loudness of various tones when masked 40 db 
were calculated and are seen to be in good 
agreement with the observed data. 


BEARING OF THESE LOUDNESS RELATIONS 
UPON QUESTIONS CONCERNING THE 
HEARING PROCESSFS 


During the discussion of the material dealt 
with in the previous sections it has been neces- 
sary to frequently refer to the processes involved 
in hearing. It should be emphasized here again 
that it was assumed that a sound which produces 
a constant agitation level above threshold level 
at all nerve patch positions produces also a 
constant masking for pure tones of any fre- 
quency. By means of this assumption it was 
shown that the experiment on masking by 
thermal noise leads to a method of determining 
the relation between f and x which was given 
in Fig. 7. The position coordinate x was not 
taken as distance along the basilar membrane, 

















TABLE III. 
POSITION x! Count x 
Upper Middle plus Apical | 0-7 | 4200| 0-14.5 
Lower Middle 7-12 | 5600 | 14.5-33.8 
Upper Basal 12-19} 8500 | 33.8-63.1 
Lower Basal 19-30 | 10700 | 63.1-100 
Total 0-30 | 29000 0-100 
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Fic. 21. Relation between the position of maximum re- 
sponse along the basilar membrane and the frequency of 
the stimulating tone. 


as has usually been done, but it is a coordinate 
which designates a particular set of nerves. 

In the interval between x=zero and x=1, one 
percent of the nerve endings which are stimulated 
by low tones is meant; at the position between 
x=1 and x=2 is meant one percent of the nerve 
endings adjacent to the first patch. In other 
words, the x values refer to nerve patches in such 
a way that for equal intervals dx there are an 
equal number of nerve endings. If the nerves are 
distributed uniformly on the basilar membrane 
then the x as used in this discussion corresponds 
to one percent of the total length of the basilar 
membrane. 

If, however, the nerve endings are not dis- 
tributed uniformly, then the length dx must be 
stretched or contracted depending upon whether 
the nerve density is less or greater than the 
average. Let x’ be the distance from the helico- 
trema in millimeters. 

Guild!® has made measurements of the number 
of ganglion cells along the basilar membrane. 
He divided the cochlea into four parts, namely, 
(1) upper middle and apical, (2) lower middle, 
(3) upper basal, (4) lower basal. He made a 
count of the total number of ganglion cells in 
each of these parts. They are given in column 3 
of Table III. Steinberg has estimated the dis- 
tance in millimeters along the basilar membrane 


10 Guild, ‘Correlations of Histologic Observations and the 
Acuity of Hearing,’ Acta Oto-Laryng. 17, 207 (1932). 
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Fic. 22. Relation between the relative intensity of a 
thermal noise producing uniform masking and the loudness 
it produces. 


for each of these parts to be that given in the 
column under x’. The corresponding values of x 
are then easily calculated as the percent of the 
total cells up to the position corresponding to x’. 
They are given in the last column. 

It is seen that these data on the number of 
ganglion cells allow us to calculate only three 
points besides the two end points. These values 
of x and x’ are plotted on Fig. 20 and a smooth 
curve is passed through the plotted points. 
Using the relation shown by this curve then the 
values of x and f from the solid curve of Fig. 7 
were transferred to values of x’ and f. The values 
obtained in this way are shown in Fig. 21. The 
relation between f and x from measurements on 
the minimum perceptible differences in pitch 
discussed above is shown by the dotted line in 
this figure. The close agreement of these two 
curves is very good evidence of the validity of the 
interpretation given to the experiments on mask- 
ing and loudness. 

The curve of Fig. 9 gives the loudness when 
the entire membrane is stimulated uniformly. In 
Fig. 22 it is replotted, shifting the abscissa so 
that for the loudness N=1 the intensity J=1. 
It will be seen that there is a range of loudness 
values of about 1,000,000 from the threshold to 
the highest level that the ear can stand. At low 
levels the most sensitive nerve endings begin to 
send discharges to the brain. As the intensity 
increases more discharges occur from these 
sensitive nerve endings and also more endings are 
excited. The curve of Fig. 22, according to my 
interpretation, is a combination of these two 
effects. Possibly the steep part of the curve is 
due principally to the first process and the 


constant slope part is due to the second process. 

It is generally supposed that each nerve fiber 
cannot discharge more rapidly than 1000 times 
per second. The 30,000 nerves can then carry 
30,000,000 discharges per second. The number V 
which corresponds to the highest intensity that 
the ear can stand is about 1/30 this value, which 
seems to indicate that at the threshold there are 
about 30 nerve discharges per second. If loudness 
were proportional to the number of active fibers, 
as some have supposed, then our loudness 
numbers could reach only 30,000 instead of 
1,000,000. This, it seems to me, is evidence, 
although not conclusive, that the total number 
of nerve discharges going to the brain is pro- 
portional to the corresponding number on our 
loudness scale. 


SUMMARY 


It has been shown how to define loudness and 
loudness level in a quantitative way. Definite 
procedures are given for determining experi- 
mentally the loudness level of any sound heard 
by any person. For a typical observer a true 
loudness scale is developed. The relation of the 
scale to the loudness level scale is determined 
experimentally. The scale has been found to 
be very useful for calculating loudness from the 
noise spectrogram, the noise audiogram, or the 
overtone structure of the sound. 

The relation between the masking and the 
loudness produced by a sound has been quan- 
titatively determined and a formula deduced 
from this relation which has proved useful for 
calculating the loudness. This formula may be 
applied with equal success to a normal ear and 
also to a deafened ear. Evidence has been given 
that the masking expressed in decibels produced 
upon any pure tone is equal directly to the 
agitation of 1.1 percent of the total nerve endings 
expressed in decibels above the threshold value 
for such a patch and at the position where such 
a tone would be sensed. These loudness relations 
throw light upon some of the important processes 
involved in hearing. In particular the data from 
the masking effects of thermal noise were used 
to calculate the relation between the position of 
maximum stimulation on the basilar membrane 
and the frequency of the tone producing the 
stimulation. 
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Limits of Audition for Bone Conduction 


NorMAN A. WATSON 
University of California at Los Angeles 
(Received January 31, 1938) 


The laboratory facilities and apparatus described in an earlier paper were used to study the 
normal intensity and frequency limits of audition for bone conduction. The condition of 
the observer’s head to give greatest acuity for pure tones, and the optimal total force of applica- 
tion and area of the vibrator button were found to be essentially the same as for speech sounds, 
Acuity was found to vary with the position of application of the vibrator; for example, the 
maximum variation with position of the 1000-cycle threshold was 18 db. A threshold for open 
canals was determined for the frequency range 80-2000 cycles, without interference from stray 
air radiation (18 normal observers). Above 2000 cycles air radiation interfered with open 
canal tests to such an extent that a special technique was devised to obviate its effects. The 
occluded threshold curve, for a single normal individual, lay below his open canal curve over 
the range 80-2000 cycles, and also at 10,000 cycles. The usable intensity range for open canals 
was tentatively determined; it isa maximum of 80-90 db at 1000-2000 cycles. The total audible 
frequency range for open canals was found to be at least 25-17,000 cycles. 


INTRODUCTION 


HE normal frequency and intensity limits of 

hearing by air conduction have been deter- 
mined accurately by many observers. Fletcher 
and Munson! have delineated the auditory 
sensation area for binaural listening in a free 
sound field (for young observers with normal 
hearing). Several investigators have determined 
threshold of audibility curves for sound produced 
in or near the external auditory canal and 
listened to monaurally. A discussion of their 
various methods and an average curve have 
been given by Sivian and White.? Wegel*® has 
determined the threshold of feeling for monaural 
listening. These curves of Sivian and White and 
of Wegel have been extrapolated by Knudsen* 
to give an auditory sensation area for monaural 
air conduction for young persons with normal 
hearing. However, comparatively little work has 
been done on the determination in absolute units 
of .the limits of audition for bone conduction, 
although relative thresholds of hearing, suitable 
for testing of hearing, have been obtained by 


1H. Fletcher and W. A. Munson, J. Acous. Soc. Am. 5, 
82 (1933), Fig. 4; also see H. Fletcher, Bell System Tech. J. 
13, 239 (1934), Fig. 1. 

2L. J. Sivian and S. D. White, J. Acous. Soc. Am. 4, 
288 (1933), Fig. 10. 

?R. Wegel, Proc. Nat. Acad. Sci. 8, 155 (1922); see also 
H. Fletcher, Speech and Hearing, (D. Van Nostrand, 
New York, 1929), p. 141, Fig. 70. 

*V. O. Knudsen, Audiometry (The Sonotone Corpora- 
tion, Elmsford, New York, 1937), p. 6, Fig. 2. 
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several investigators. Two threshold of hearing 
curves for bone conduction in absolute units 
have been determined over limited frequency 
ranges: one by Knudsen and Jones,° in terms of 
the amplitude of vibration of the vibrator 
button ; and one by von Békésy® in terms of the 
effective amplitude of vibration of the forehead 
skin compressed over the frontal bone near the 
point where the vibrator button was applied. 


The objectives of the present investigation . 


have been to determine a threshold of hearing 
curve for bone conduction over as wide a 
frequency range as possible without interference 
from stray air radiation, to estimate the useful 
intensity range for tones within that frequency 
range, to determine the total audible frequency 
range, and to plot a tentative auditory sensation 
area for bone conduction. 


LABORATORY FACILITIES, APPARATUS AND TECH- 
NIQUE OF TESTING 


All the measurements reported in this paper 
were made in the small soundproof test room of 
the acoustical laboratory of the University of 
California at Los Angeles, which has been de- 
scribed in an earlier paper.” Thus the tests were 
free from the influence of external noise and of 


®V.O. Knudsen and I. H. Jones, Arch. Otolaryng. 13, 
489 (1931), Fig. 3. 
‘ ®°G. von Békésy, Ann. d. Physik (5), 13, 111 (1932), 
ig. 10. 
7N. A. Watson, J. Acous. Soc. Am. 9, 99 (1937). 
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adventitious noise arising within the test room. 
The induction type electrodynamic vibrator® 
described in the previous paper was used for 
most of the tests. The driving system for the 
vibrator is indicated by the block diagram of 
Fig. 1. It consists of a General Radio beat 
frequency oscillator, a meter to measure the 
input current, a voltage divider, attenuators and 
amplifiers, one of which contains resonance type 
filter networks. The frequency response charac- 
teristics of the vibrator and the driving system 
are known so that appropriate corrections can be 
made for their deviation from uniformity. Both 
the vibrator and driving system are essentially 
linear over the intensity range used in the tests. 
The largest overtone in the vibrations of the 
bone conductor for frequencies above 100 cycles 
is 35 db below the fundamental, and most of the 
harmonics are 50 db or more below. The second 
harmonic of 80 cycles is only 11 db below the 
fundamental; thus too much weight cannot be 
given the measurements at that frequency. 

The method of determining a threshold was as 
follows: The observer was seated in the test 
room with the vibrator against his head and a 
signal button in his hand. The operator sat at 
the controls in the adjoining room.’ By means of 
them he could change the intensity of the test 
tone or interrupt it momentarily. The interrup- 
tion was brought about by breaking the filament 
supply line to the last stage of the second ampli- 
fer. The tubes in this stage were not of the 
heater- type; thus the filaments cooled and 
heated again very rapidly so that the tone 
regained its original intensity almost immedi- 
ately after the filament supply was reconnected 
following the interruption. The plate supply line 
and the audio line between amplifiers remained 
connected; hence there was no “‘click,” either 
when the tone was interrupted or when it was 
started again. The operator first presented the 
tone to the observer at a level definitely above 
threshold. The observer signaled as long as he 
heard it. Then the operator interrupted the tone, 
shifted the attenuator setting to that for a less 


* The general appearance of the vibrator and its mode of 
application to the head are shown by Fig. 7 of reference 7; 
a schematic cross sectional diagram is given in Fig. 8 of 
the same paper. 

*The test room and control room are shown as room 2 
and room 1, respectively, in Fig. 4 of reference 7. 
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intense tone, and presented the fainter tone to 
the observer. (Attenuator changes were made 
during the time the tone was interrupted so no 
contact noises would reach the observer.) This 
process was repeated and the observer’s responses 
noted each time until he failed to signal at two 
intensities decreasing successively. Then the in- 
tensity was increased until he again signaled 
that he heard at two intensities increasing 
successively. The procedure was repeated until 
the operator found the maximum attenuator 
setting (faintest tone) for which the observer 
heard the tone a majority of times it was pre- 
sented. The effective amplitude of vibration of 
the vibrator button corresponding to this attenu- 
ator setting was called the threshold for that 
observer for the frequency tested. No regular 
sequence was followed in the time-length of 
tones or interruptions, as the responses might 
easily be affected by rhythmic times of tone 
duration and interruption. With all ‘‘clicks’’ and 
other extraneous noises eliminated and irregular 
tone intervals used, the operator could be sure 
that the tone corresponding to the threshold 
setting was the real threshold tone for the ob- 
server at the time of the test. 


EFFECTS OF VARIOUS FACTORS ON THE BONE 
CONDUCTION THRESHOLD 


Before measurements were made for the deter- 
mination of the average threshold of hearing 
curve, the effects on the threshold of various 
controllable factors such as the condition of the 
observer’s head, the position of application of 
the vibrator, the force of application of the 
vibrator and the contact area of the vibrator 
button were studied. It was found that acuity 
was greatest when the mouth of the observer was 
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Fic. 1. Block diagram of}the driving system for the electro- 
dynamic bone conduction vibrator. 
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closed and his teech touching or clenched, thus 
making the mandible practically an integral part 
of the vibrating skull instead of an accessory 
part, vibrating out of phase with it. Measure- 
ments taken with the external auditory canals 
occluded indicated a very marked lowering of 
the threshold below that for open canals. Thus 
it was decided to perform most of the measure- 
ments with open canals but later to make similar 
measurements for occluded canals and to com- 
pare the results. 

The threshold varied markedly with change of 
position of application of the vibrator button. 
For one observer, for example, the 1000 cycle 
threshold varied 18 db with change of position. 
The position for greatest acuity was the upper 
front teeth; that for poorest acuity was along 
the top midline of the head approximately six 
inches back from the forehead. With the vibrator 
button placed at the center of the forehead, 
just above the position of the frontal sinuses, 
the acuity was 10 db less than with the vibrator 
on the upper front teeth, and 3 db less than 
with the vibrator on the mastoid portion of the 
temporal bone.!° In spite of this, the center of 
the forehead was chosen as the position of appli- 
cation of the vibrator button because (1) it is 
the most convenient portion of the head for that 
purpose, (2) the frontal bone at that point is 
very uniform in density and elasticity over 
quite a large area, allowing comparatively wide 
latitude in the location of the button without 
introducing serious errors in the measurements, 
and (3) since hearing by bone conduction is 
essentially a binaural phenomenon, it is ad- 
vantageous to place the vibrator cenfrally to 
provide equal distribution of the vibrations to 
the two ears when testing persons of normal 
hearing. 

The effect of the variation of the force of 
application of the vibrator button was tested by 
taking thresholds with forces of application 
varying from 85-650 grams-weight for fre- 
quencies as low as 125 cycles, and as high as 
2500 cycles, which represented the highest fre- 
quency that could be used without interference 
from stray air radiation. The variation in 
threshold with force of application was found to 


10 Compare with the differential given by N. H. Kelley, 
J. Exp. Psych. 21, 211 (1937), Table III. 
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be small for all frequencies for forces greate 
than 370 grams-weight—that is, great enough to 
compress the skin sufficiently to allow it to 
vibrate in unison with the bones. When the force 
was decreased the variation in threshold was 
found to be small for high frequencies (within 
+3db). For low frequencies, however, the 
variation was greater; the threshold was raised 
10 db for 125 and 250 cycles by the reduction of 
the force of application to 85 grams-weight, 
On the basis of these tests the force of application 
of 370 grams-weight was chosen for further tests, 

The effect of variation of the area of vibrator 
button when applied with a force of 370 grams- 
weight was next studied. For both open and 
occluded canals, the variation of threshold with 
button area at low and medium frequencies was 
found to be practically negligible (a maximum 
variation of +3 db). For closed canals, however, 
the larger button areas gave much greater acuity 
(that is, lower thresholds) at high frequencies, 
For 3000 cycles, the threshold intensity in db 
decreased 10 db almost linearly with increase in 
button area up to 4.75 sq. cm. For 7000 cycles it 
decreased 24 db up to the same button area. 
Accurate measurements could not be made at 
the higher frequencies for open canals because 
of interference from stray air radiation. Buttons 
of area less than 0.5 square centimeter were too 
uncomfortable for the observer when a force of 
application of 370 grams-weight was used. A 
compromise among all these factors led to a 
choice of button area of 0.665 square centimeter. 

The unit chosen for measuring the threshold 
was the effective (that is, root-mean-square) 
amplitude of vibration of the vibrator button 
(0.665 sq. cm area) when pressed against the 
center of the observer’s forehead with a force of 
370 grams-weight. 


THE THRESHOLD OF HEARING CURVE FOR 
BONE CONDUCTION 


To determine the threshold of hearing curve 
for bone conduction, thresholds were taken for 
eighteen observers within the age range 20-28 
and of average age 23. Their air conduction 
hearing lay within +5 db of normal, as measured 
with an accurately calibrated audiometer. They 
were instructed to place their chins on the sponge 
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20 500 
FREQUENCY~CYCLES PER SECOND 


Fic. 2. Threshold curves for bone conduction (open ear 
canals). Curve 1, average threshold of hearing curve 
(present investigation). Curve 2, threshold of hearing 
curve (Knudsen and Jones). Curve 3, threshold of hearing 
curve (von Békésy). Curve 4, tentative ‘‘safe-upper-limit” 
curve (present investigation). (See area 1, Fig. 5.) 


rubber chin rest and to keep their mouths closed 
and their teeth touching, but not clenched. For 
each observer, thresholds were taken at fre- 
quencies of 80, 125, 250, 500, 1000, 1500 and 
2000 cycles per second. At frequencies greater 
than 2500 cycles the stray air radiation interfered 
to an extent which prevented the determination 
of accurate thresholds. For the frequencies listed 
above, the stray air radiation thresholds meas- 
ured with the vibrator just off the head, were at 
intensities at least 20 db above the bone con- 
duction thresholds, and thus did not interfere 
with them. The thresholds of the eighteen ob- 
servers at each frequency were averaged and 
each arithmetic mean corrected for frequency 
responses of the vibrator and its driving system. 
The corrected mean values are shown in Fig. 2 
by the circles of curve 1. The variations in 
threshold among the eighteen observers at the 
test frequencies are indicated by the vertical 
bars drawn above and below the circles. The 
variation is rather large, amounting to more 
than 20 decibels at some frequencies and more 
than 10 at all of them. Similar variations have 
been reported by other workers in bone con- 
duction research. Part of the large total variation 
may be attributed to inaccuracies of determining 
individual thresholds. How large a part may be 
estimated from the following : fourteen thresholds 
at 1000 cycles were taken for a trained observer 
over a period of a month and a half with all 
conditions as nearly constant as possible. Seven 
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EFFECTIVE AMPLITUDE - CENTIMETERS 


Fic. 3. Comparison of threshold of hearing curves for 
air and bone conduction (open canals). Curve 1, average 
bone conduction curve (present investigation). Curve 2, 
bone conduction curve for a single observer (present 
investigation). Curve 3, monaural, minimum audible 
pressure, air conduction curve for a single observer 
(Waetzmann and Kiebs). Curve 4, average binaural, 
minimum audible field, air conduction curve (Fletcher and 
Munson). Curve 5, monaural, minimum audible pressure, 
air conduction curve (von Békésy). 


of the thresholds were the same, and the extreme 
range was +3 db. The average error for a single 
determination was 1.1 db. Thus it is clear that 
most of the total variations arose from the indi- 
vidual skull differences. 

The average curve for this investigation (curve 
1) lies slightly above that of Knudsen and 
Jones® (curve 2) over most of the frequency 
range covered by the latter. This separation of a 
few decibels would be expected, for their meas- 
urements were taken on the mastoid portion of 
the temporal bone while those for curve 1 were 
taken on the forehead. Curve 1 lies some 30- 
35 db above von Békésy’s® curve (curve 3). 
He measured the effective amplitude of vibration 
of the forehead skin, compressed over the frontal 
bone and vibrating with it, rather than the 
effective amplitude of the vibrator. In another 
experiment" he found that only one to four 
percent of the vibrations imparted to the skin 
were transmitted to the skull; one percent trans- 
mission would correspond to a 40 db transmission 
loss and four percent transmission to a 28 db 
transmission loss. Thus the 30-35 db difference 
is what one would expect. 

In Fig. 3, the average bone conduction curve 
(No. 1) is compared with the bone conduction 
curve (No. 2) for a single observer (whose hearing 


1! Reference 6, Figs. 7 and 8. 
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is somewhat more acute than average) and with 
various air conduction curves. Curve 3 is a 
monaural air conduction curve for a single 
observer, determined by Waetzmann and Keibs,” 
using a method correcting for the absorption 
coefficient and the transmission factor of the ear 
drum. Curve 4 is the binaural minimum audible 
sound field curve of Fletcher and Munson.! 
Curve 5 is the monaural minimum audible 
pressure curve given by von Békésy® in terms 
of the effective amplitude of vibration of a small 
capsule placed in the ear canal as a source of 
sound. The separation of the air conduction and 
bone conduction curves is determined by the 
separation given by von Békésy® between his air 
and bone conduction curves at 800 cycles and 
by the distance of the bone conduction curve of 
this investigation from von Békésy’s. Fletcher 
and Munson’s! and Waetzmann and Keib’s® 
curves were arbitrarily made to coincide with von 
Békésy’s® air conduction curve at 800 cycles. 
Curve 1 has approximately the same slope as 
the binaural minimum audible sound field curve 
over the frequency range 80-500 cycles. Above 
that frequency the curves converge. The same 
is true for the monaural curves, although the 
convergence is not as great. Thus it appears 
that the minimum audible intensity for bone 
conduction is not exactly the same function of 
frequency as that for air conduction. 

The frequency range of threshold measure- 
ments was restricted by the interference of stray 
air radiation, in spite of the precautions taken in 
the design of the bone conduction vibrator to 
avoid it. To extend the measurements to higher 
frequencies it was necessary to locate the sources 
of the stray radiation and to quiet them if 
possible; if not, to eliminate the effects they 
produced. 

The sound reaching the ears by air conduction 
was found to have two sources: (1) the sound 
radiated to the air by parts of the bone conductor 
other than the button and either conducted as 
progressive waves to the ear by diffraction 
around the head or reflected so as to give rise to 
standing wave patterns in the room with maxima 
in positions to actuate the air of the external ear 


122 E, Waetzmann and L. Keibs, Ann. d. Physik (5), 26, 
141 (1936), Fig. 1. 
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canals ; (2) the sound radiated to the air from the 
back side of the vibrator button (or from the 
observer’s forehead) and reflected from the front 
surface of the vibrator housing so as eventually 
to reach the ears, either as progressive waves or 
the maxima of a standing wave pattern. 

It was possible to eliminate the effect of these 
unwanted sounds for any one frequency by the 
expedients of adjusting the positions of the 
vibrator and the head of the observer in the room 
so that his ears were in positions of pressure 
minima of the standing wave pattern, and of 
wrapping the vibrator (particularly the front 
face of the vibrator housing) with loose, highly 
absorptive material to reduce the sound radiated 
directly and to absorb the sound reflected from 
the front face of the housing after it was radiated 
from the back side of the vibrator button and the 
observer’s forehead. The position of the head and 
vibrator were adjusted for one frequency (10,000 
cycles) by putting an extension on the vibrator 
rod, raising the vibrator by leveling screws, and 
moving the observer’s head slightly to one side, 
yet keeping the vibrator button at the same 
position of the forehead and the observer's chin 
on the chin rest. Without the special precautions, 
the open canal bone conduction threshold at 
10,000 cycles coincided with the stray air radi- 
ation threshold; with the precautions, the bone 
conduction threshold was 15 db below the stray 
air radiation threshold. To obtain thresholds at 
other high frequencies it would be necessary to 
provide vibrator rod extensions of different 
lengths, and to reorient the vibrator and the 
observer’s head in the room for each frequency. 

(An alternative expedient would be to enclose 
the vibrator in an absolutely soundproof housing 
and surround the vibrator rod and button witha 
soft rubber cylindrical shell mounted on the 
housing so that it would not touch the rod or 
button but so that its free flaring end would 
impinge on the forehead in a circle surrounding 
the point of application of the vibrator button 
and thus prevent radiation from the back side of 
the button. However, such an arrangement 
would be an extra body pressing against the head 
and might interfere with the accurate determi- 
nation of thresholds. For clinica! testing, it 
would be desirable, and has, as a matter of fact, 
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FREQUENCY - CYCLES PER SECOND 


Fic. 4. Threshold of hearing curves for bone conduction 
(open and closed ear canals, single observer). Curve 1, 
open canals. Curve 2, closed canals. 


already been employed on the bone conduction 
vibrator of a recent audiometer.") 

Numerous experimenters have found that the 
acuity of hearing by bone conduction is increased 
by occlusion of the external auditory canals. 
They have found that at low frequencies the 
separation of the open canal and closed canal 
thresholds is of the order of 10-20 db, that the 
curves converge at medium frequencies, and that 
the separation disappears at frequencies above 
2500 cycles.!* Essentially the same results were 
obtained in the present investigation when no 
special precautions were taken to obviate the 
effects of stray air radiation, as is seen in Fig. 4. 
Curve 1 is the open canal threshold for a single 
observer, curve 2 the occluded canal curve for 
the same observer. There is a separation of over 
20 db from 50 to 200 cycles, a decreasing but 
definite separation to 2000 cycles. Over this 
frequency range there was no interference from 
air radiation. Above 2500 cycles the open canal 
threshold was apparently equal to or even lower 
than the closed canal threshold ; this was caused 
by the interference of stray air radiation with the 
open canal threshold. At 10,000 cycles, when the 
special precautions described above were em- 
ployed, the open canal bone conduction threshold 
was 15 db below the stray radiation air con- 
duction threshold and 6 db above the closed canal 
threshold. Thus the phenomenon of increased 
bone conduction acuity on occlusion of the 





81. H. Jones and V. O. Knudsen, Laryngoscope 46, 523 
(1936), Fig. 6. 

“See reference 5, p. 490; A. G. Pohlman, Ann. Otol. 
Rhinol. and Laryngol. 39, 927 (1930); and reference 10. 
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external auditory canals exists at the high fre- 
quency of 10,000 cycles. 


_ THE USABLE INTENSITY RANGE AND TOTAL 


FREQUENCY RANGE FOR BONE CONDUCTION 


For the determination of the usable intensity 
range for bone conduction for the audible fre- 
quency range, a large power amplifier was used 
to drive the bone conduction vibrator, and the 
intensity of vibrations increased to the point 
corresponding to the maximum undistorted out- 
put of the amplifier. The amplitudes produced 
under these conditions are shown in Fig. 2 by 
curve 4. There was no feeling of pain (similar te 
that elicited by intense air borne sounds) pro- 
duced for any of the observers at any frequency 
from 50-12,000 cycles when the vibrator was 
applied to the center of the forehead for short 
intervals. However, an application of the vibrator 
at the high intensity of curve 4 at 1000 cycles for 
several minutes produced for one observer an 
unpitched tinnitus (sounding like a rushing air 
stream) which lasted two days. Hence the 
intensities indicated seem to be near the limit of 
the usable intensity range, at least for the 
middle and higher frequencies. Experience with 
the testing of speech sounds in which most of the 
energy is concentrated between 200 and 300 
cycles, and with various massage vibrators (low 
frequency) indicates that the usable intensity is 
greater at the low frequencies. 

From the author’s experiments and those 
conducted by others, the audible frequency range 
for bone conduction may be set at 20-20,000 
cycles. 


THE AUDITORY SENSATION AREA FOR 
BONE CONDUCTION 


On the basis of the data just given, a tentative 
auditory sensation area was constructed for bone 
conduction for open canals. This is shown in 
Fig. 5 as area 1. It is compared with the monaural 
auditory sensation area (No. 2) based on the 
measurements of Sivian and White? and of 
Wegel ;3 also with the binaural free field auditory 
sensation area (No. 3) given by Fletcher and 
Munson.! The placement of the pressure level 
scale (in terms of which area 2 is plotted) with 
respect to the effective amplitude scale (in terms 


— 
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@ INTENSITY LEVEL IN 08 FROM 10° WATT PER SQ-CM. 


FREQUENCY ~ CYCLES PER SECOND 


Fic. 5. Auditory sensation areas. Area 1, tentative bone 
conduction (present investigation). Area 2, monaural, 
minimum audible pressure, air conduction (Sivian and 
White, and Wegel, extrapolated by Knudsen). Area 3, 
binaural, minimum audible field, air conduction (Fletcher 
and Munson). 


of which area 1 is plotted) was made on the 
relative placement of the Békésy air conduction 
curve and the Sivian and White minimum 
audible pressure curve in Sivian and White’s 
original article.!* The placement of the intensity 
level scale (in terms of which area 3 is plotted) 
with respect to the pressure level scale was 
determined by the relative placement of the 


6 Reference 2, Fig. 9. 


minimum audible pressure data of Sivian and 
White and the binaural free field data by Fletche 
and Munson." The bone conduction area is seen 
to be smaller than either of the air conduction 
areas, having approximately the same frequency 
range as the latter but a smaller usable intensity 
range. However, it is large enough to include the 
frequencies and intensities of speech sounds g 
that speech may be heard as well by bone 
conduction as by air conduction by persons with 
normal hearing!’ and so that persons with hearing 
impairments of the conductive type may utilize 
bone conduction hearing aids to good advantage, 

The author wishes to acknowledge his jp. 
debtedness to Professor V. O. Knudsen fo, 
suggesting this research and for invaluable help 
and guidance in carrying it out; to Mrs. Ruth 
B. Watson and Mr. R. B. Watson, who assisted 
in many of the tests ; to Mr. L. W. Sepmeyer and 
Mr. Henry Jung for assistance in the design and 
construction of some of the apparatus used, 
Acknowledgment is also made to the department 
of physics of the University of California at Los 
Angeles for use of its facilities and apparatus and 
for its financial support of the research. 


16 Reference 1, Fig. 2. 
17 Reference 7, Fig. 11. 
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Periods of Longitudinal Vibration of Steel Cones and Truncated Cones 


G. W. PIERCE AND ATHERTON Noyes, Jr. 


Cruft Laboratory, Harvard University, Cambridge, Massachusetts 


HIS paper presents the results of some 
measurements of the periods of vibration of 
solid steel cones and truncated cones, excited 
into longitudinal vibration by magnetostriction. 
A theoretical investigation is given together with 
a comparison of observed and theoretical re- 
sults. Harmonic as well as fundamental modes of 
vibration were studied. 


Modes of longitudinal vibration of the cone 


The complete solid cone, illustrated at the 
top of Fig. 1, was made to vibrate continuously 
in its natural fundamental longitudinal mode, 
with a loop of motion at both ends, and its 
frequency f; was determined. The same vibrator 
was then vibrated in its first dynamic overtone, 
with a loop of motion at both ends and a third 
loop of motion intermediate between these two 
end loops. The frequency of this mode of vibra- 
tion is designated fe. Likewise the same vibrator 
was made to execute its higher dynamic over- 
tones with frequencies f3, fs, etc. This cone is 
designated by specifying its length /, its diameter 
D,, with D,=0. 

The experiments extend also to the funda- 
mental frequencies and overtone frequencies of 
a series of truncated solid cones having various 
lengths / and various diameters D, and D, as 
shown in the lower illustration of Fig. 1. 


Materials, dimensions of the vibrators, and 
measured frequencies 


Two samples of steel were employed in the 
measurements: first a nickel-steel with a sound 
velocity of 461,200 cm/sec.; and second, a so- 
called stainless steel with a sound velocity of 
538,400 cm/sec. The first happened to be the 
better magnetostrictive vibrator and gave higher 
harmonics than the second, but was used up by 
injudicious cutting so that we did not obtain 
from it as many separate truncated specimens 
as we might have done with a better planned use 
of the material. With the second specimen we 
proceeded more carefully in successive reductions 
of length, but were limited somewhat in the 
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harmonic range by the quality of the material. 


The experimental results of our measurements 


are given in Table I. 


Method of measuring these frequencies 


The values of the frequencies of vibration of 
these vibrators were determined with the aid of 
the twin magnetostriction oscillator, described 
in a recent paper! by the present authors. The 
frequencies were obtained by measuring with an 
audiofrequency meter the difference between the 
frequencies of the conical vibrators and the 
known frequencies of a set of standard magneto- 
strictive rods. These measurements of fre- 
quency were not made with extreme accuracy but 
are generally better than 0.1 percent. 


A comparison of observed frequencies and 
theoretically computed values 


It is shown below by a theoretical investiga- 
tion that the frequencies are related to the 


diameters of the vibrators by the equation 


r(2fl/c) r2fl\?  Da/Ds 
) 


tan (2xfl/c) (D./Ds—1)?’ 


in which /=length of the vibrator, whether a 
cone or a truncated cone, 
D,=diameter of the large end, 
D,=diameter of the small end, 
f=frequency, whether fundamental or 
overtone, 
c=velocity of longitudinal sound waves 
in the cone. 


c 





Fic. 1. 
1 J. Acous. Soc. Am. 9, 185 (1938). 
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This quantity c is a function of the ratio of condition is not fulfilled is,” 
diameter to length of the specimen. The function 
is not known exactly for the cones, but for a 
uniform cylinder, if v=the velocity obtained 8? 

with a rod of diameter small alia miata with  * Ber euperimentat verification see David S. Muzzey Jr 
its length, the velocity c in cases where this Phys. Rev. 36, 935 (1930), which gives references to 


n2xo?2D? |} 
=o[1- approximately; (2) 


TABLE I. Materials, dimensions and observed frequencies. 








OBSERVED FRE- 
LENGTH | LARGE DIAMETER | SMALL DIAMETER | QUENCIES CYCLES 
MATERIAL CM Da CM Db CM PER SECOND 





Nickel-Steel, v=461,200 cm/sec. 14.08 1.261 0.139 21,030 
36,520 
51,800 


67,360 
83,200 


98,870 








Nickel-Steel, v=461,200 cm/sec. 22,810 
41,810 
61,110 
80,700 
99,960 

118,960 

137,990 








Nickel-Steel, v = 461,200 cm/sec. 27,150 
51,800 
76,690 


101,380 








Stainless Steel, v= 538,400 cm/sec. ‘ 27 24,160 
41.470 
58,450 





24,310 
42,100 








24,480 
42,680 
24,640 
43,170 














24,880 
43,970 





25,290 
45,170 





26,110 
47,500 











28,120 
52,400 








30,250 
57,380 





35,100 
67,940 





























52,910 

















ey, Jr., 
Ices to 


Da/Dp 
9.075 


3.68 


2.33 
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where 2=an integer defining the mode of vibra- 
tion, 

D=diameter of the cylindrical rod, 

«= Poisson’s ratio. 


The second term in the brace of (2) is small in 
comparison with unity. 

In the case of the cones and truncated cones 
we shall assume that we may replace D? of the 
cylindrical case by (D2+D,?)/2, giving 


nno?(Da2+D,2) )3 
16/? 


c=vi1— , approximately. (2a) 





This assumption is probably not quite correct 
over the whole range of our observations as 
given in Table I. 

Equation (1) is a theoretical equation giving 
2fl/c (theoretical) as a function of the diameter- 
ratio D,/D,; while from our observed values of 
Table I we have calculated 2f//v, and by the'use 
of (3) obtained 2fl/c (observed) for various 
given values of D./D». 

It is somewhat simpler to compare 2f1/nc, ob- 
served and calculated, as is done in Fig. 2. 
In the curves of Fig. 2 the diameter-ratio 
D,/D, is plotted as abscissae and 2f1/nc calcu- 
lated by Eq. (1) as ordinates, giving the lines of 
the curves, for »=1 (fundamental), »=2 (first 
dynamic overtone), etc. The corresponding ob- 
served values of 2f1/nc are plotted as crosses or 
circles. The crosses are for the nickel-steel 
specimens and the circles for the stainless steel 
specimens. 


Agreement between observed and calculated 
frequencies 


It is seen in Fig. 2 that for the fundamental 
mode of vibration (n=1) and for the first dyn- 
amic overtone (n=2) the observed frequencies 
fall on the calculated curves with the exception 
of the points at D,/D,=2.33 where the observed 
frequencies are lower than the calculated by 
about 0.2 percent for n=1 and n=3. This 
departure, while a little greater than the errors 
of frequency measurements, is probably not 
greater than the inhomogeneity of the material 





theoretical treatment by 


4 Pochhammer, Chree, Love, 
Rayleigh and others. 
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or the errors of measurement of mechanical 
dimensions of this specimen. 

On the other hand, for the higher dynamic 
overtones, »=3 to n=6, and especially for large 
values of D,/D, the observed values are lower 
in frequency than the calculated values by as 
much as 0.5 percent, whence it appears that 
the equations (1) and (2) though approximately 
correct are not quite exact for high overtone 
vibrations of truncated cones which have a large 
ratio of diameter of large end to diameter of 
small end. 


Frequencies of complete cone 


We attempted to prepare a complete cone; 
i.e., to make D,=0. In this case Eq. (1) becomes 


tan r2fl/c=n2fl/c. (3) 
The roots of this equation’ are 


2fl/c= 1.4303, 2.4590, 3.4709, 
4.4774, 5.4818, 6.4844 (4) 


corresponding to n=1, 2, 3, 4, 5, 6, respectively. 

With our attempted complete cone we ob-_ 
tained experimentally for n=1, 2, and 3, re- 
spectively, 


2f1/c(observed) = 1.384, 2.375, 3.355. (5) 


The specimen was made of the stainless steel 
and would not oscillate at the higher overtone 
modes. It is seen that the observed values (5) are 


3Love, The Mathematical Theory of Elasticity, third 
edition, (Cambridge University Press, 1920) p. 289. 
Rayleigh, Theory of Sound, second edition, Vol. 1, (Mac- 
millan, 1894), p. 334 has also these values (4) but with 
4.4774 misprinted 4.4747. 
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less than the computed values (4) by about 3.5 
percent. However, mechanical measurements 
showed that the small end of the specimen was in 
fact not a true point with a uniform taper to the 
end. It is found that if we compute from the 
observed frequency relations of (5) by Eq. (1), 
the values for the diameter-ratio D,/D,, we 
obtain D,/D,= 30.4, 28.6, 28.7, respectively, for 
n=1, 2, and 3 of Eq. (5). This means that the 
observed and calculated values of 2fl/c agree 
if the cone actually had an effective small-end 
diameter D, equal to about 1/29 of D,=0.044 
cm. Actually the tip diameter was certainly less 
than this, but calculations from other geo- 
metrical measurements showed that a complete 
cone having the general taper of this one should 
have had a length of 15.90 cm, whereas the 
actual length was only 15.42 cm. That is, the 
cone tapered somewhat more rapidly near the 
tip than elsewhere. Now, had the taper been 
perfect with a resulting length of 15.90 cm, the 
frequency would have been very nearly un- 
changed. See Table I for frequency changes due 
to cutting off the cone tip. Hence 2fl/c would 
have been about 3 percent greater than the 
values we obtained by using a value of 15.42 cm, 
and this would have brought the values of (5) 
into close agreement with (4) 


Theoretical derivation of Eq. (1) 


In a brief search of the literature we have not 
been able to find the general Eq. (1). Accordingly, 
we here give a theoretical treatment of the 
problem, based on parts of Chapter XIII of 
Love’s Mathematical Theory of Elasticity In 
that chapter on page 297, for propagation of 
waves in elastic solid media, Love gives the 
equation of motion in the form 


aA A dA 
+i) (=.=) +unu,2, 2) 
Ox dy dz 


( Cu dv ew 
=F , ’ 
0? of ar 


In this equation u, v, w are components of dis- 
placement (not velocity) and 


Ou ov dw 
A=—+—+— (7) 
Ox Oy a2 
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as pointed out on p. 281 and p. 299 of Love's 
book. By differentiating the three equations of 
(6) respectively by x, y, and 2 and adding the 
results Love obtains 


(A+2y)V2A = p(02A/d?*). (8) 


Love next points out that, if certain curl com. 
ponents vanish, this equation may be written 


CV" = 0o/d0?, (9) 
where 


c=((A+2y)/p)! (10) 


is the velocity of propagation of the disturbance 
and ¢=a potential function, not the usual veloc. 
ity potential of Rayleigh, but a function whose 
gradient is the particle displacement, with com- 
ponents designated in Love as u, v, w. 

Keeping Love’s definition of ¢ we shall change 
his notation for displacement and let =the 
particle displacement in the direction of r drawn 
from the vertex of the cone. 

Then for motion only in the direction r, 


£=09/0dr. (11) 


Love points out on p. 298 that, if the motion is 
irrotational and if ¢ is a function of 7 and ¢ alone, 
it satisfies the equation 


f(r—ct) F(r+ct) 
o= . eee 


r r 


, (12) 


with f and F denoting arbitrary functions. c is 
the velocity of the wave. 

If now as usual we restrict these arbitrary func- 
tions to sinusoidal functions of the time we have 


Asinw(t—r/c) Bsinw(t+r/c) 
¢=<—— + rg 


r r 


Love next discusses the components of traction 
per unit area designated by him —X,, — Y,, —Z, 
in any normal direction drawn toward the sur- 
face, and obtains at the bottom of his p. 300 


Ou Ov Ow 
( ’ ; =—-—(X,, Y,, 2»). 
ot ot at 


In our notation, since these normal traction com- 
ponents become pressure, we may write 


p=pc(d€/dt). (14) 
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Now employing Eqs. (11), (14), and (13) we 
obtain 


p=pc(d/dt)(d/dr) 
—A cosw(t—r/c) Aw sin w(t—r/c) 
= ca] +- 





9 


ia cr 





Becoswlitr/c) Bw sin w(t+r "| 
r? cr 


Expanding the sums and differences in this 
equation we obtain 





pcw A+B wr w or) 
| ——— cos ——(A+B)- sin a cos wt 
r r c c c 
A-B wr w wr 
-| sin ——(A—B)- cos sin wi (15) 
r c C c 


If now r=a at the large end of the cone, p must 
be zero at r=a; whence 
either 
I. A+B=0 and _ tan wa/c=wa/c 
or (16) 
I. A—B=0 and _ tan wa/c=—c/wa. 
If also r=b at the small end of the cone, » must 
be zero at r=b; whence 
either 
I. A+B=0 and _ tan wb/c=wb/c 
or (17) 
Il. A—B=0 and _ tan wb/c=—c/wb. 
We shall now combine (16) and (17), using the 


LONGITUDINAL VIBRATION 








OF CONES 





combination formula 


tan x—tan y 





tan (x—y) = ey 
1+tan x tan y 


obtaining from I of (16) and I of (17), 


wa wh wa/c—wh/c 
tan (—-~)-~——. (18) 
Cc 1+ (wad) /c? 


or, likewise, if we use II of (16) and II of (17), we 


obtain 
wt wh —c/wa+c/wb 
tas (— -—) NO cet 
= <a 1+c?/ (wad) 
wa/c—wb/c 


1+ (wab/c?)’ 


which is the same as (18). 
It only remains to change the form of Eq. (18) 
by use of the notation indicated in Fig. 1; 
b/a=D,/Di, a-—l=b, w=2rf. 
which substituted into (18) gives, on inversion, 


D./D> 


a(2fl/ w2fl\?2 
fil) 1+(—) ——_————. (19) 
{D./De—1}? 


tan (r(2fl/c)) 7 c 


This is the Eq. (1), which we have employed 
above and found (cf. Fig. 2) to be approxi- 
mately checked by the experimental values of 
frequencies f for different lengths / and diameter 
ratios D,/D, of the conical vibrators. 


Method of calculation of curves of Fig. 2 


The curves drawn in Fig. 2 were calculated 
from Eq. (1) using the value of ¢ given in (2a). 
It is seen that if D,/D», 1 and v are given the 
solution of (1) for frequency f is not simple. On 
the other hand if / and v are known we may 
introduce various values of f and compute D,/D». 
This latter method was adopted and 2f,//nc 
was then plotted vs. the computed values of 
D,/D, to give the curves of Fig. 2. 


For a given ratio of the end diameters the 
dynamic overtone frequencies are nearly 
linearly related to the overtone number n 


It is an interesting and curious result that 
values of f,/f1, plotted in Fig. 3 against the over- 
tone number ” give nearly a straight line for 
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each fixed value of D,/D,. Except in the case of 
the top curve for a cylindrical rod (D,/D,=1) 
these nearly straight lines do not pass through 
the origin. This means that with the cones and 
truncated cones as vibrators and for any given 
value of D,/D, the overtone frequencies are not 
integral multiples of the fundamental frequency, 
but that the difference of frequencies of suc- 
cessive overtones (fn4i1—fn) of a given cone or 
trunctted cone is nearly a constant quantity for 
emi, 2, 3,4.... 

The top and bottom curves of Fig. 3 were 
calculated from Eq. (1), with omission of the 
slight correction factor for the change of c with n 
as given in Eq. (2a). The top curve for D,/D,=1 
is accurately straight with slope unity; and the 
bottom curve for D,/D,=infinity is given by 
the values of Eq. (4), from which, on division by 
the first number 


fo/fr=1, 1.719, 2.427, 3.130, 3.833, 4.534 
with differences 
(frii—fn)/fr=0.719, 0.708, 0.703, 0.703, 0.701. 


The three intermediate lines of Fig. 3 are 
taken from the curves of Fig. 2 computed from 
Eq. (1), and the circles on the lines are the ob- 
served values with the nickel-steel set of trun- 
cated conical vibrators. Similar results were also 
obtained for the set of stainless steel vibrators. 
The worst departure from rectilinearity is 
in the set of computed values of the preceding 
paragraph. 

This approximately linear relation is inherent 
in Eq. (1) as well as in the observed vibration 
frequencies, as may be shown independently 
but with some difficulty by a mathematical 
examination of Eq. (1). In this examination it 
is best to employ the series expansion of an 
antitangent. The treatment of this problem is 
not here given. 


Comparisons of theory and results with theory 
of conical air columns 


In one of the sections above we have pointed 
out the absence of literature on solid conical 
vibrators. On the other hand there is consider- 
able theoretical material on the vibration of 
conical air columns. Lord Rayleigh,‘ assuming 


4 Rayleigh, Theory of Sound, Vol. II, (1894), p. 281. 
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certain limitations as to smallness of diameter 
of a truncated conical pipe filled with air, qwiy, 
both ends of the conical pipe closed, finds the. 
oretically for the air column an equation that 
leads to our Eq. (1). The condition imposed by 
Lord Rayleigh is that the velocity of the vibrating 
column is zero at the ends of the column, whereas 
with our solid steel vibrators it is the pressup 
at the two ends that is zero. In this case of 
pressure zero at the two ends Lord Rayleigh 
obtains for conical air columns the same resylt 
for frequency as for cylindrical air columns; 
namely (in our notation), sin (wl/c)=0. This 
is far from true of our vibrators, except in the 
special case of D,/D,=1, when our equations 
give the same results. Barton’ and Stewart and 
Lindsay® have arrived at the same results as 
Rayleigh. Barton also describes briefly some 
elementary experimental work in confirmation 
of the theory of Rayleigh. 

Of course there is a difference between our 
vibrators and the air-column vibrators of Ray- 
leigh, in that his air-columns were constrained 
against sideways velocities by the conical tubes 
containing the air, while our rods were unre- 
strained against sideways motions. It seems, 
however, that neither Rayleigh and others nor 
we have introduced the condition of restraint 
or unrestraint into the theoretical treatment. 


Summary 


We have checked by experiment the relation 
given by Eq. (1) for the vibration frequencies of 
cones and truncated cones, as a function of their 
end diameters. For the fundamental and first 
overtone the agreement is very good; for the 
higher overtones, for which the diameter of the 
vibrator is more important in comparison with 
the wave-length of the vibration, the theory 
is less closely checked by experiment, though the 
departure is not great. 

We have found that for cones, truncated 
cones and uniform cylinders the frequency 
differences between successive overtones 1s 
nearly constant, but that only in the case of 
narrow uniform cylinders are the overtone fre- 
quencies exact integral multiples of the funda- 
mental frequency. 


§ Barton, Textbook on Sound, (1919) pp. 251-262. 
6 Stewart and Lindsay, Acoustics, (1930) pp. 79-83. 
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The fundamental frequency of a truncated 
cone of fixed length is always higher than the 
frequency of a uniform cylinder of the same 
length and material, and is a function of the 
ratio of the end diameters ; the ratio of the funda- 
mental of a truncated cone to the fundamental 
frequency of a cylindrical rod of the same length 
and material changes from 1.000 to 1.4303 as the 
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ratio of the end diameters changes from 1 to 
infinity. The curves of Fig. 2 relate these fre- 
quency ratios, for fundamental and for overtones, 


_ to the diameter ratios. 


A theoretical derivation of the fundamental 
equations and a discussion of the departure of 
the results from theoretical results by Rayleigh 
and others for conical air columns are given. 
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The Role of the Speaker Impedance in Resonance in a Closed Pipe 


WituiaM D. PHELPS 
Research Division, RCA Manufacturing Company, Camden, New Jersey 


(Received August 31, 1937) 


An analysis is made of a standing wave resonance system consisting of an eight-foot iron 
pipe with an eight-inch dynamic speaker at one end and a movable iron piston at the other. 
Considering the speaker both as a source and as an impedance, expressions are derived for the 
phase shift, ¢:. of the pressure wave at the speaker, and for the pressure amplitude in the tube. 
The relation connecting / and \ at resonance is found to be /=(n+¢/27)d/2, placing the 
speaker at a distance (¢/47)\ from the first particle velocity node. (Contrast this equation 
with the relation /=mn)/2 for a tube rigidly closed at both ends.) Initial resonance frequencies 
and a pressure-frequency curve are computed for comparison with experimental data. The 
possibility of using the resonance conditions for measuring the impedance of diaphragms is 


pointed out. 


INTRODUCTION 


N order to obtain high sound pressures for 


investigating amplitude distortion in micro- - 


phones, a standing wave resonance system was 
built. It consists of an eight-foot iron pipe ten 
inches in diameter, and a five-sixteenths inch 
wall with an eight-inch dynamic speaker at one 
end and a movable iron piston at the other (see 
photograph, Fig. 1). The microphone is sup- 
ported from the piston. Use of the apparatus 
revealed that the relation between the length / 
of the tube and the wavelength \ at resonance 
is certainly not given by the formula /=n)/2, 
where 7 is a whole number 1, 2, 3, etc. Also, a 
rapid variation of pressure at resonance with 
resonance frequency in the region of the speaker 
resonance was noted. Present treatments of the 
problem! of forced vibrations in a tube have 
failed to properly take into account the effect 
of the source on the resonance phenomena. 


THEORY 


Before outlining a mathematical analysis of 
the problem, an intuitive picture of resonance 
in the tube will be considered. The speaker at 
any time /) sends down the tube a pressure pulse 
of phase 6. The pulse reaches the end of the tube 
closed with a rigid piston, suffers no phase change 
on reflection, and arrives back at the speaker at 
a time ¢ in the same phase @. If, on reflection at 


1 Stewart and Lindsay, Acoustics, p. 57. 
2Lamb, Dynamical Theory of Sound, p. 174. 
3 Crandall, Theory of Vibrating Systems and Sound, p. 103. 
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the speaker, it suffers a phase change ¢, it now 
is of the phase 6+ ¢. But in order for re-enforce- 
ment with the pulse now going out from the 
speaker and resonance to occur, the speaker must 
go through the phase n27+ ¢, n integral, while 
the original pulse is traveling the length of 
the tube and back. This reasoning leads to 
the conclusion that w(t—¢to) =2wl/c=n2r+¢ or 
l=(n+¢/2r)\/2 as the resonance condition, 
where c is the velocity of sound and the other 
terms have been previously defined. Further 
analysis is necessary to determine the function 
y, the pressure amplitude in the tube, and to 
verify the relation connecting length and wave- 
length at resonance. 

Assuming undamped plane wave propagation 
in the tube (see Fig. 2) let the particle velocity 
and pressure in the wave traveling to the right 
be represented by 


£,= Aye” ei", (1) 
Pi=pocAie**e'*'" (2) 


and the corresponding quantities in the wave 
to the left be represented by 

£o = Aget ikteivt | (3) 

a= pocAgetitzeiot, (4) 

One wishes to know the phase shift in the 

incident pressure ~2 when reflected at the speaker. 


Since the amplitude of p2 is reduced on reflection, 
and since* | p2| =| p1|, it follows that p; is made 


4 See Eq. (13). 
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up of two components. If p;’ represents the con- 
tribution to p; of the speaker as a source and 
pi’ the pressure reflected by the speaker as an 
impedance, then for a steady state 


pi'+pi" =p (5) 


and the phase shift is given by the expression 
(pi’’/p2)2-0. For convenience in evaluating p1''/p2 
one may, without loss of generality, assume for 
the moment that p;'=0, i.e., that the voice coil 
current is zero, then p;’’=(1),~0 and hence from 
Eqs. (2) and (4) 


(pi'’ / p2)2m0 = Ar” / Ad. (6) 


Also from the boundary condition on volume 
currents at x=0 


S2(A1”’ — Ao)e*#*=Siés, (7) 


where és is the cone velocity. The equation of 
motion for the speaker is 


bi’ +he= —ZaSiés, (8) 


where Za, is the speaker impedance in acoustic 
ohms including damping due to the electrical 
circuit. 
Again from Eqs. (2), (4), and (7) Eq. (8) 
becomes 
poc(Ay"’+Ae)e* = —ZaySo(Ar’ —Aa)e', (9) 
whence 


(pr / pe) 20 = (Za1— Za) /(Za+Zai) (10) 
or 
(pi / p2) 2-0 = Be’, (11) 


where 
g=tan-! (2Z4X 4;)/(Rav+Xar—Za’*) 


IMPEDANCE 


Fic. 1. The standing wave resonance system. 
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is the phase shift in pressure amplitude on reflec- 
tion at the speaker, Z4 = poc/ Se, Z4,;=Ra,+iXay, 
and 


B= [(Rai—Za)?+Xar }'/[(Rai+Za)l?+X a1}! 


is the ratio of the reflected to the incident pres- 
sure amplitudes. 

As the next step, the pressure amplitude in the 
tube as a function of frequency, speaker imped- 
ance, and pipe constants will be determined. 
From the conditions imposed on the volume 
currents at the two ends of the tube, one finds at 


x=0 

S2(A1—Ae)e@* = Sis. (12) 
At x=! 

(A ye~**! — A ge**!) iat = 6 
or Ag= Aye" *!, (13) 


One should note that while Eqs. (7) and (12) 
are of the same form, they are for different con- 
ditons; namely, for F=0, and for F= Foe’, 
respectively. Also, the equation of motion for 
the speaker becomes 


F= S?2(Z41—-1Z4 cot kl)és, (14) 


where —iZ, cot k/ is the input impedance of the 
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Fic. 2. Diagram of the standing wave resonance system. 
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Fic. 3. Determination of the initial resonance frequen- 
cies of pipe and speaker, Ra: <ZA. ¢ is the phase shift by 
the speaker. 


closed tube.’ From Eqs. (12), (13), and (14) one 
finds 


Ay =e! Fy/2S,So[4 sin kl Zait+Za cos ki]. (15) 


Hence the pressure amplitude at any distance x is 
| (Pitp2) |e 
poc cos k(l—x) Fo 
~ |$,SefiRay sin kl-+Za cos kl— Xa; sin kl] 


. (16) 


From Eq. (16) one may easily derive the 
resonance condition on k, /, and ¢. For a given 
value of k the pressure amplitude is a maximum 
at a pressure loop when / is such that the absolute 
value of the complex denominator is a minimum. 
Considering the denominator as a function of / 
and minimizing, the resonance condition is found 
to be 


tan 2k) =2Z4X 4,/(RaPv+Xa"?—Za?). (17) 
Comparing Eqs. (17) and (11) it is seen that 
(18) 


(19) 


tan 2kl=tan ¢ 
or 2kl=2nr+ ¢. 


Eq. (19) derived on the basis of maximum ampli- 
tude is identical with the relation arrived at 
intuitively on the basis of phase alone. It empha- 
sizes the importance of phase reenforcement of 
the reflected and generated pressure at the 
speaker for high pressures in such a system. 
Having derived equations showing the effect 
of the speaker on the resonance frequencies of 
the system, and on other aspects of the problem, 


5 Olson and Massa, Applied Acoustics, p. 35. 
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Fic. 4. Calculated pressure-frequency characteristic for 
the first four peaks. The resonance frequencies are accu- 
rately given by the phase change on reflection at the 
speaker. Note the peak value trend in the region of the 
speaker resonance frequency. 


it is instructive to graphically illustrate some of 
these by curves. If the pipe length and cross- 
sectional area and speaker impedance are known, 
then Eq. (19) and the expression for ¢ in Eq. (11) 
enable one to compute the resonance frequencies. 
In Fig. 3 these are given by the intersection of 
the curves ¢ and 2kl—2nz, both considered func- 
tions of kl, for the first four resonance frequencies 
of the system. Table I shows a comparison of the 
calculated and measured values of k/ at reso- 
nance. With the exception of the first one, the 
agreement is not bad.°® 

Using Eq. (16) the frequency characteristic 
(see Fig. 4) of the pressure amplitude at a loop, 
say at x=/, has been plotted for ki in the region 
of the speaker resonance frequency assuming 
constant voice coil current. The speaker reso- 
nates at kl=0.787r. It is to be noted that the 
magnitude of the peaks increase with k while 


TABLE I, Comparison of calculated and measured kl at 
resonance, kl=1.42rv107. 








MEASURED 


0.87 

1.387 
2.207 
3.242 


CALCULATED 


0.627 
1.307 
2.207 
3.167 





* In this and subsequent computations the contribution 
to Ra, due to damping in the cone suspension has been 
neglected. 
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SPEAKER 


the minima decrease. While this peak value de- 
pends on a number of factors as given by Eq. 
(16), the rise is explicable in the fact that the 
phase shift by the speaker decreases with in- 
creasing k, causing the cone to approach a par- 
ticle velocity node. The decrease in the minima 
is due to the increasing speaker impedance. In 
Fig. 5 is plotted, for comparison with Fig. 4, a 
pressure-frequency characteristic obtained with 
a pressure microphone placed at x =/—2.54 cm. 
No corrections have been applied. With the 
exception of the first peak the curves in the two 
figures are similar. 


DISCUSSION 


One is now in a position to answer the ques- 
tion: Is the speaker at a particle velocity node? 
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Fic. 5. Experimental pressure-frequency characteristic of 
the system described taken with a pressure microphone. 
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SPACE AMPLITUDE OF $ 


Fic. 6. Representation of the space amplitude of the 
particle velocity along the tube at resonance, »=2. The 
speaker is located at x=0. 


The answer is that it is not in general. The 
distance from the speaker to the first node is 
equal to (¢/41)\. This condition is illustrated 
in Fig. 6, for example, where n= 2. 

It should also be noted that the speaker ab- 
sorbs some of the incident wave as well as changes 
its phase on reflection. 

Since the speaker impedance enters into the 
resonance phenomena, it is to be expected that 
resonance might be used to measure the imped- 
ance of diaphragms. An examination of Eq. (17) 
substantiates this statement. By measuring R, 1, 
and Za, Ra; and Xa; may be computed. 


CONCLUSIONS 


1. The speaker shifts the phase of the reflected 
wave by an angle ¢ thereby making the reso- 
nance condition 1=(n+¢/2z)d/2 in contrast 
with the condition /=m\/2 for a tube rigidly 
closed at both ends. 

2. The speaker is at a distance (¢/47)A from 
the first particle velocity node. 

3. The pressure amplitude is dependent on 
speaker impedance and phase shift. 

4. The resonance Eqs. (17) and (18) offer a 
method for measuring speaker and microphone 
impedance. 
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Absolute Sound Intensity in Liquids by Spherical Torsion Pendula 


EviAs KLEIN 
Naval Research Laboratory, Washington, D. C. 


(Received August 18, 1937) 


These experiments demonstrate the use of spherical radiometers to determine the energy 
density in an acoustic field. It is shown that the theoretical work of King is applicable to the 
measurement of radiation pressure in liquids. Considering the torsion pendulum as an absolute 
instrument, pressure microphones are calibrated as secondary standards. 





INTRODUCTION 


HE object of this paper is to describe a 

method for absolute sound intensity meas- 
urements in liquids; and to indicate means of 
calibrating pressure microphones in terms of 
fundamental units. This investigation is based 
entirely upon the well-known phenomenon of 
acoustic radiation pressure which postulates a 
nonlinear character to the wave progressing in 
a fluid medium. 

In the realm of high frequency acoustics, 
particularly in liquids, the radiometer is almost 
indispensable as a measuring instrument. This 
device, in its various forms, responds to the 
steady continuous pressure existing in a sound 
field. The alternating pressures commonly asso- 
ciated with sound propagation phenomena are 
only indirectly involved in these measurements. 
The relation between the two types of pressures is 


(y+1) p 
R,= oe 
2 pe 





(1) 


where R, is the radiation pressure which a 
radiometer evaluates; y is the ratio of the 
specific heats; p is the r.m.s. oscillatory pressure ; 
and pc is the acoustic resistance of the medium in 
which the pressures are measured. In air R,+0.8 
X10-*p?, while in water R,+0.44X10-"p*. In 
other words, one dyne per square cm of oscillatory 
pressure in air yields less than a microdyne of 
radiation pressure; whereas, in water the same 
unit of r.m.s. pressure (1 dyne/cm?) produces less 
than half a ten-thousandth of a microdyne in 
radiation pressure. This indicates the magnitudes 
which are involved in radiation pressure ex- 
periments. 

Radiation pressure is associated with the 
momentum of wave motion and is measured by 


the impact of the wave upon an obstacle. One 
may regard radiation pressure as a rectified 
component (second order in magnitude) of the 
oscillatory motion which a vibrating source 
imposes upon the medium. This increment in 
static pressure in a sound field is due to the 
nonlinear relation existing between pressure and 
density. Indeed, if it were possible to generate an 
acoustic wave disturbance and so influence the 
medium that it would respond only to a linear 
law of force, then there should be no momentum 
propagated and no radiation pressure detectable. 
However, the results of many investigators! have 
demonstrated not only the presence of radiation 
pressure, but have verified its quantitive relation 
to the energy density per unit volume, E, of the 
medium, namely 


Ryp=(y+1)E/2. (2) 


THE PRESENT PROBLEM 


The general procedure of all radiometric de- 
terminations is to place an obstacle in the path of 
the wave train, which obstacle causes a change in 
the propagated momentum and itself experiences 
a reaction. Usually, this disturbs the sound field 
and makes absolute measurements unreliable. A 
primary requisite of any measuring device is that 
its influence upon the domain under observation 
shall be either entirely negligible or exactly 
calculable. Therefore, the obstacle must be of 
such geometric form as to introduce (a) no 
noticeable distortion into the sound field; or 
(b) its reflective and refractive effects should be 
definitely known. Both phases of this problem 


1 Dvorak, Pogg. Ann. 57, 42 (1876); Meyer, Phil. Mag. 6, 
225 (1878); Zernow, Ann. d. Physik 21, 131 (1906); Ray- 
leigh, Scientif. Papers, V, 41 and 262; Weaver, Phys. Rev. 
15, 400 (1920); Altberg and Holtzman, Physik. Zeits. 26, 
49 (1925). 
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PISTON SOURCE 


STANOARD 
SIGNAL 
GENERATOR 


Pn 


CRYSTAL UNIT 


Fic. 1. Radiometer assembly. 


are being investigated ; however, for the present,” 
the discussion is confined to spherical torsion 
pendula as related to (b). 

Accordingly, sound energy is allowed to fall 
upon solid spheres of various sizes comparable to 
the wave-length in the medium. The acoustic 
field is symmetrically distorted by the obstacle. 
The latter reflects some of the incident energy, 
while a portion of it may unite behind the sphere 
so as to form a shadow. An exact knowledge of 
the scattering and diffraction effects of spheres 
when subjected to acoustic radiation pressure has 
been made available in a theoretical development 
by L. V. King.’ Following this mathematical 
treatment precise calculations have been made 
herein for the total force exerted upon spherical 
obstacles in plane stationary wave systems. The 
results of these computations are compared with 
a series of experiments which show reasonable 
quantitative agreement with King’s analysis. 


EXPERIMENTAL METHOD AND APPARATUS 


The methods of producing plane stationary 
waves in liquids at high frequencies are quite well 


2 Correlation of the two methods is already in progress 
with encouraging results, as demonstrated at the Washing- 
ton meeting, May 3, 1937. 

3L. V. King, Proc. Roy. Soc. A147, 212 (1934). 


known. Nevertheless, certain of the precautions 
taken are noteworthy. Fig. 1 is a composite 
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Fic. 3. 
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Fic. 4. 


sketch of the principal parts of the apparatus. A 
nonresonant crystal quartz sound generator 
(piston type) fills one end of a rigid cylindrical 
tube, six inches’ in diameter. The other end, 
about 2 feet away is closed by a machined 
reflector. The latter is so arranged as to function 
either as a complete reflector or as a reflector 
minus a central plug. The microphone to be 
calibrated is placed in the hole left by removing 
the plug, about 1}’’ diameter (Fig. 2 and Fig. 3). 
The tube and associated parts. are completely 
immersed in a tank of distilled water. Acous- 
tically the length of the tube is adjustable by 
providing telescopic’movement between the pis- 
ton source and the rigid tube. Also the micro- 
phone position can be regulated in the reflector 
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Fic. 5. 


opening so that the pressure upon the microphone 
is a maximum. By this means, sharply defined 
nodes and antinodes are established within the 
cylinder. These are essential in absolute determi- 
nations of radiation pressures in standing wave 
systems. Moreover, the nodes and antinodes 
should be fixed in space relative to the closures of 
the tube. To maintain these conditions, the 
temperature of the water tank, which houses the 
apparatus is kept fairly constant; also, the 
frequency and voltage of the electrical oscillator 
which actuates the sound source are carefully 
controlled. 

A spherical torsion balance (Fig. 4 and Fig. 5) 
is employed to explore and evaluate the radiation 
pressure generated within the tube. The pressure 
is exerted only upon the measuring sphere which 
is kept at all times exactly on the axis of the 
cylinder. The sphere visible above the cylinder 
acts merely as a counterpoise. The center of each 
sphere is at a fixed distance from the axis of 
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suspension. A force causing the measuring sphere 
to turn is balanced by the torsion of the wire. A 
scale and mirror are employed to indicate the 
angular deflections. The photograph shows the 
1’ spherical suspended system ready to be placed 
into the sliding portion of the tube. Alongside are 
the 3’’ spheres to indicate the plug fitting and to 
show the only opening of consequence through 
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which energy from the inside of the tube may 
escape to the surrounding water. 

Before extending the description of these 
experimental details, it is reasonable to inquire as 
to the need for such elaborate precautions to 
minimize energy flow from the tube. The neces- 
sity for eliminating leakage and diffusion of 
sound through openings was proved by tests with 
a series of slots parallel to the axis of the tube. 
The slots‘ (one such is shown in Fig. 6) varied in 
width from 4/10 to /4. The results showed 
considerable distortion and asymmetry in the 
standing wave system inside the cylinder due to 
energy transmission through the apertures. This 
was evident from the radiation pressure curves 
when plotted in a manner similar to Fig. 7, 
Fig. 8 and Fig. 9. The ideal condition should 
present a sinusoidal distribution of R, along the 
axis of the cylinder. Unfortunately, the errors 








‘In this connection it is interesting to recall that the 
amount of sound energy which passes through an opening 
is considerably out of proportion to the area of the opening. 
rit see Wintergest and Knecht, V.D.I. 76, 777, 
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introduced into radiation pressure measurements 
by the presence of pseudo-stationary waves due 
’ to an apparent attenuation in the medium are 
not readily calculable. Therefore, apertures of 
any sort were carefully avoided in the design of 
the present cylindrical enclosure. Also a serious 
effort was made to have as nearly a perfect 
standing wave system as is practically realizable 
in distilled water contained in a tube.® 

The latter is over 1} inches in thickness and is 
composed of three concentric layers of different 
metals. The inner shell is of ;;’’ brass while the 
outer shell is 3’’ bronze. Hot molten lead was 
poured between these two cylinders so as to form 
a tubular sandwich of brass-lead-bronze which 
machined as a single cylinder about 3 ft. long. 
The sector shown in Fig. 6 as the sliding portion 

5 Radial oscillations may become very troublesome at 
certain frequencies despite the tube design and construc- 
tion. (Field, Can. J. Research 5, 131 (1931).) In these tests 
it was necessary to provide an inner lining for the main 


cylinder in order to break up at some frequencies supposed 
radial oscillations. 
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of the main cylinder was cut from the original 
length of tubular-sandwich, and is about 6” 
longer than the lower portion of the tube over 
which it travels on guides. The rectangular plug 
through which the spheres are joined (Fig. 5) is 
also part of the original cylinder. Bronze spheres 
were used in these experiments. A nickel tube 
about 0.05 inches in diameter joins the two 
spheres. The suspensions are 0.001” tungsten 
wire supporting the }’’ sphere and 0.002” 
carrying the 3’ spheres. 


MEASUREMENTS AND DATA 


In these experiments acoustic radiation pres- 
sure data serve a twofold purpose. The energy 
density of the incident wave may be evaluated. 
At the same time a pressure microphone is 
calibrated as a secondary standard. The two sets 
of observations are made simultaneously by inde- 
pendent observers. Deflections of the spherical 
torsion pendulum are noted at different points 


TABLE I. 
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along the axis of the cylinder. These readings, 
obtained either by using the torsion head or by 
the direct deflection method, are taken one to 
three mm. apart depending upon the frequency 
of the emitted sound. Examples of such data are 
shown in Fig. 7, Fig. 8 and Fig. 9. 

The voltages impressed upon the microphone 
located in. the reflector plate are measured while 
the deflections of the sphere are recorded in its 
various positions. Fig. 10 illustrates the general 
relationship existing between the received volt- 
ages and the current amplitude of the source at 
constant frequency. The three curves indicate 
three different locations of the radiometer, and 
suggest the likely effect of the suspended system 
position upon the r.m.s. pressure developed at the 
reflector. 

From Fig. 10 it will be noted that the piston- 
source amplitude varies directly with the alter- 
nating pressures as registered by the microphone. 
On the other hand, the same current-amplitude 
of the piston source is proportional to the square 
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P is the total force acting on the sphere of radius 
@; po and p; are the densities of the liquid and 
sphere respectively. |A| is the complex ampli- 
tude of the wave; h, the distance from the center 
of the sphere to the nearest velocity node; 
k=2n/d wave-length; and a=ka. F,, G, and H, 
are related 

H,2= F,2+G,?, (4) 


where /7,” is expressible in a Bessel series. (3) may 
be written as 


P=rpo|A|2A sin 2kh, (5) 


where A represents the numerical equivalents of 
the series in the square bracket [ _] for different 
values of (a=ka). 

Now the mean total energy density in the 
incident wave is 


E=}pok?| A’|. 
Therefore, P=(2rE/k*)A sin 2kh. (6) 











root of the radiometer deflection, shown by the 
experimental curve in Fig. 11, and implied in 
Eq. (1). 

Table I represents averages of considerable 
data gathered primarily for the purpose of 
testing King’s theory, and to ascertain the 
suitability of the spherical torsion balance for 
absolute measurements. Although a sufficient 
number of check determinations by other 
methods has now been made to establish these 
results quantitatively, yet they lack in ultimate 
precision and refinement when regarded as 
standard values. Therefore, the different columns 
merely indicate certain observed and calculated 
quantities. They also show estimates of other 
magnitudes dealt with experimentally. 


THEORY AND CALCULATIONS 


The general expression for the mean radiation 
pressure on a sphere in a plane stationary field, as 
developed by King is: 


Pi 
(n+1) (PasiPatGniiGn) 


ents 13 es 





jat—n(n+2)} (3) 





Expression (6) is used to evaluate experimental 
results. A was calculated in terms of Bateman’s or 
Bessel functions, and is represented graphically® 
in Fig. 12. Another form of (3) for the total force 
on a sphere is 


P=F=rxa’Ee sin 2kh, (7) 


where o is a function of (a2=ka) and the density 
ratio (po/p1). 

Calculations of o in Eq. (7) were carried out by 
Sivian and O'Neil’ for the measurement of 
aerial sound intensities. Comparing (6) and (7), a 
relation between o and A is obtained, viz. 


o=(2/a2)A. (8) 
Fig. 13 shows a comparison of the Bell Labora- 


6 The calculations were made and checked with the 
helpful cooperation of Dr. M. Sasuly and Mr. H. Hedinger. 

7 The author is indebted to the Bell Laboratories for 
permission to publish expression (7) and the solid curve in 
Fig, 13 for comparison . 
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tories calculations of o with its values derived 
from (8). It should be observed that the solid ¢ 
curve was computed with the assumption that 
po/p1=0, quite legitimate for aerial waves. In 
the A calculations, the experimental value of 
(po/pi=1/8.8) was used. Different series and 
distinct methods were employed in evaluating o 
and A, respectively. These curves bring to light 
the experimental fact that in a stationary wave 
much larger radiation pressure is exerted on a 
sphere whose value of a is relatively small. 
Theory also predicts that for the same maximum 
amplitude of sound, the measuring sphere of 
small a experiences a considerably greater magni- 
tude of radiation pressure in a standing wave 
system than it does in a progressive wave. For 
larger values of a, however, the conditions are 
reversed. That is, radiation pressure is greater in 
a progressive wave. 

To reduce experimental data, either equation 


(6) or (7) is suitable provided, of course, the 
values of A and o are available for different 
ranges of a. Although (6) is consistently used in 
the present calculations, the physical significance® 
of the factors involved is more apparent in (7). 
Rewriting the latter in terms of radiation pres- 
sure, R,, we have 

R,=F/xa?= Eo sin 2kh. (7’) 
This expression states first of all that in a plane 
stationary wave in which the mean energy 
density is E, the radiation pressure acting upon a 
spherical obstacle is generally less (depending 
upon ¢ which is always less than unity) than the 
energy density in the neighborhood of the sphere 
because of scattering and diffraction. 

Secondly, the radiation pressure is a periodic 
function of the position of the sphere relative to 
the particle velocity nodes and loops. Moreover, 
the periodicity is double that of the standing 
wave pattern. This dependence of radiation 
pressure upon the location of the sphere was 
theoretically predicted by King and amply veri- 
fied in these experiments. 

In reducing the data to numerical values, the 
average of the radiation pressure amplitude was 
used for a given series of deflections. Thus (6) 
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Fic. 13. Solid curve, Bell Laboratories calculation; broken 
curve. N.R.L. reduction from ¢ =2A/a?. 


8 For further explanation see Appendix. 
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P=(22rE/k?)A. (9) 
Now the torque acting upon the sphere due to 
radiation pressure is balanced by the twist in the 
suspension after the system turns through a 
resultant angle 6 or 










Pl=sé, (10) 










where / is the distance from the center of the 
sphere to the axis of rotation; and s is the torsion 
constant of the suspension. Combining (9) and 
(10) ; and solving for the energy density, we have 


E=s6k?/2xrAl=27256/WAl. (11) 


It is reasonable to expect that E would remain 
fixed while operating at one frequency with the 
same power into the piston source; and, main- 
taining all other experimental conditions con- 
stant (except the torsion balance size). That 
being true, the different diameter spheres should 
experience a torque according to 













$161/Ail1 = S262/Adle, (12) 


where the subscripts 1 and 2 correspond to the 
factors associated with the }’’ and }” spheres, 
respectively. Rewriting, there results 


(Sol1/S ile) 62/6; = Ae/ Ai. (13) 






The coefficient of the ratio of deflections is made 
up of predetermined experimental constants, the 
value of which is 8.50. Therefore, the validity of 
(3) may be tested for various frequencies by the 
aid of (13), provided that E does not change 
during the investigation of any one frequency. 
This procedure was followed for the frequencies 
recorded in the table. Columns 8 and 9 show the 
approximate agreement between theory and ex- 
periment. By a series of motional impedance 
measurements associated with this project, as 
well as by certain piezoelectric studies, the as- 
sumption that E remains constant was further 
justified. Column 10 suggests the limits of alter- 
nating pressures exerted upon the crystal micro- 
phone located in the reflector plate while column 
11 gives the approximate voltages developed by 
several microphones. 
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DISCUSSION AND CONCLUSION 








It is desirable to point out some likely sources 
of error which were not considered in this series 
of experiments. Their effects upon the absolute 
values of radiation pressures are still to be ascer- 
tained. Consequently, estimates of accuracy are 
omitted until a number of influencing factors are 
made more determinate. 

Although noticeable in most of the curves 
similar to Fig. 7, Fig. 8 and Fig. 9, the attenua- 
tion of the signal in the medium between the 
source and radiometer was neglected in the com- 
putations. Incidentally, this type of radiometer 
might be readily adapted to absorption measure- 
ments in liquids. 

Nor was any attempt made to correlate the 
size of the tube with the high frequencies used for 
optimum radiation pressures. The data available 
need more study. Furthermore, absolute results 
obtained from acoustic measurements in liquids 
contained in tubes with beamed sources require 
careful interpretations. 

The spheres used were assumed to be “‘rigid.”’ 
Preliminary calculations based upon the physical 
properties of the spheres showed that neither 
their spheroidal nor their radial modes of natural 
mechanical vibration frequencies were close to 
the excited acoustic frequencies. Nevertheless, 
the discrepancies observed between theory and 
experiment at 75 and 100 kc lead one to suspect 
that a nearby mechanical resonance in the }” 
sphere may have caused some trouble. 

These and other attending difficulties notwith- 
standing, the method described for absolute 
acoustic determinations in liquids gives a plaus- 
ible answer to the question ‘‘How Much Sound?”’. 


APPENDIX 


The physical significance of an equation, however per- 
tinent, is not quite acceptable as an explanation of a 
phenomenon which is vaguely understood. Thus, while dis- 
cussing the paper at the Washington meeting, there arose 
two questions concerning radiation pressure in stationary 
fields. 

(1) If a stationary wave is equivalent to two progressive 
waves of the same frequency and amplitude, traveling in 
opposite directions, why is the sphere (or any other sym- 
metrical obstacle), when placed in this sound field, not 
subject to equal and opposite forces? 

(2) How does a standing wave system of one frequency 
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Fic. 14. Relation magnitudes of pressure components on 
sphere in successive positions with reference to the station- 
ary wave system. 


and presumably free from harmonics cause the sphere to 
react with twice the frequency of the sound generated? 
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The following is an attempt at a graphical description 
of the phenomena fully portrayed in mathematical 
language by King. 

Consider two plane progressive (identical) waves moving 
in opposite directions. Let these constitute the resolved 
components of the actual stationary wave system. Each of 
these progressive longitudinal waves carried with it a 
certain amount of momentum which is always directed 
away from the source. This is sketched in Fig. 14 as a 
rectified component of the wave train. Each group of seg- 
ments may be regarded as an increment in static pressure 
of the medium. In magnitude, it fluctuates from point to 
point with the wave. 

Now suppose that an infinitely thin metallic disk (say 
a section of a larger sphere) of sufficient wave front area to 
be considered an obstacle were used as a radiometer. Equal 
streams of momenta would fall on each surface of the disk, 
due to the oppositely traveling waves, and the obstacle 
should remain at rest. In other words, the rate of change of 
momentum on either side of the disk is equal and opposite 
and no motion will result. However, if the same disk be 
extended to finite thickness, comparable to the wave- 
length, then there will be positions along the wave train 
where the two faces of the disk, normal to the direction of 
propagation, are subjected to unequal changes in momenta 
and motion will ensue. The occurrence of such points is 
depicted by the aid of a half-inch sphere in successive 
locations on the axis of a plane stationary field. The arrows 
represent the component magnitudes of radiation pressures 
exerted. Those on the left with two barbs, for the incident 
wave; and, on the right the single barbs show the reflected 
wave. The resultant effect of the two is plotted in the 
lower portion of the drawing. 

Following this graphical analysis it becomes apparent 
how the periodicity of the pendulum deflections is twice 
that of the standing wave system. The same conclusion 
might be drawn from considerations of the energy density 
in a stationary wave. For, let an incident plane wave, of 
sine-law distribution and energy density E, be reflected at 
the end of a tube to form a standing wave system. Upon 
plotting the energy density of the stationary wave against 
distance along the axis of the tube, it will be found that: 

(1) The energy density in the stationary wave fluctuates 
between zero and 4E. 

(2) The fluctuations have double the periodicity of the 
stationary wave. 





- 


62 h&a ith @Ooenmea 6 © aw ee ff. 





otion 
tical 


ving 
ved 
ch of 
it a 
=cted 
asa 
" seg- 
ssure 
nt to 


(say 
ea to 
“qual 
disk, 
itacle 
ige of 
osite 
sk be 
vave- 
train 
on of 
ienta 
its is 
ssive 
‘TOWS 
sures 
ident 
ected 
| the 


arent 
twice 
usion 
nsity 
re, of 
ed at 
Upon 
ainst 
that: 
ates 


f the 





APRIL, 





1938 





VOLUME 9 





Tones Produced by a Wire Placed in an Ignited Gas Jet 


J. J. Coop 
Indiana University, Bloomington, Indiana 
(Received August 16, 1937) 


It was found that a wire placed in an ignited gas jet 
would produce a tone when the velocity of the gas reached 
a critical value. The tone was amplified by use of a second 
wire and by use of a photoelectric cell. The relation between 
the diameter D of the wire, the diameter O of the orifice, 
the distance d of the wire from the orifice and the critical 
velocity U; for the initiation of the tone was found to be of 
the form OU;=kd/D+C, where & and C are constants. 
This equation is not linear for values of d less than one 
centimeter. It appears that a similar relation exists for the 
critical velocity Uy at which ‘“‘flaring”’ starts. The products 
OU;/v, where v is the kinematic viscosity, were found to be 
nearly the same for the two gases used. At any distance of 
a wire from the orifice the frequency of the tone was found 
to be a linear function of the efflux velocity of the gas. The 
expression for the frequency was found to be of the form 

=k/D(U— Ub), where k is a constant and Us is the 
velocity intercept. The quantity D(dN/dU) has an average 
value of about 0.047 for wire diameters between 0.04 and 


HILE’ studying the action of a sensitive 

flame the writer found that a wire placed 
in the gas flame would, under certain conditions, 
produce a musical note. A study of this phe- 
nomenon was made and the purpose of this 
paper is to describe the results of this investi- 
gation. 


THEORY 


The tones produced by an obstacle placed in a 
fluid stream are commonly known as Aeolian 
tones. In the ordinary Aeolian tone phenomena 
the wire is moving with respect to a fluid of 
infinite extent. In the case here considered the 
wire is in a stream flowing into a similar fluid 
at rest. 

In either Case vortices are shed alternately 
from the sides of the obstacle and are carried 
down the stream in two rows. This procession 
of alternate vortices was first noted by Mallock! 
and studied independently by Bernard? in 1908. 
The stability of this vortex system has been 
investigated mathematically by Karman.’ 


1 Mallock, Proc. Roy. Soc. 9, 262 (1907). 
2 Bernard, Comptes rendus 147, 839, 970 (1908). 
’ Karman, Ann. d. Physik 44, 801 (1914). 


0.1 cm. For a constant frequency the relation between the 
orifice velocity and the distance was found to be of the 
form U= Kd", where n is approximately 3; nO} is a constant 
and KO! nearly constant. An approximate relation between 
frequency, velocity, and distance is given by N=kU/d!+C, 
where & and C are constants. A thin metallic sheet placed 
against the wire on the downstream side prevents the 
production of the tone. As the distance between this sheet 
and the wire is varied the tone ceases at a critical distance, 
which is a function of the velocity and the diameter of the 
wire. For a fixed velocity this critical distance is propor- 
tional to the square root of the diameter of the wire. It was 
concluded that the tone is an Aeolian tone modified by the 
flow of one stream into a similar fluid at rest. An equation 
of the form N=kU/D was derived on the basis of the 
Bernard-Karman vortex theory. By comparing experi- 
mental with theoretical results it was concluded that for 
an ignited jet the density varies inversely as the distance 
from the orifice. 


In Fig. 1 let 1 be the distance between two 
consecutive vortices in the same row, and the 
separation of the two rows by h, and let a be 
the ‘‘stagger’”” of one row upon the other. 
Karman found that there were two conditions of 
equilibrium for the system, namely: a=0, and 
a=1/2. The system is unstable for the case a=0, 
while for the case a=//2, maximum stability 
exists for h/1=0.28. 

The application of the Karman theory to the 
case here considered is very similar to that for jet 
tones given by Kruger and Schmidtke.* The fre- 
quency of the tone is the number of vortices 
leaving one side of the wire in one second. This 
is given by N=u/l where N is the frequency, 
u is the velocity of the vortices relative to the 
wire, and 7 is the distance between consecutive 
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Fic. 1. Parallel rows of alternate vortices. 


4 Kruger and Schmidtke, Ann. d. Physik 60, 701 (1919). 












vortices in the same row. Now u must be pro- 
portional to the stream velocity U so that 
u=aU where a is a constant. Also the distance 
between consecutive vortices must be propor- 
tional to the diameter of the wire so that /=bD 
where D is the diameter of the wire and 3d is a 
constant. Then the frequency N is given by 


N=u/l=aU/bD, 
or ND/U=a/b=a constant. (1) 


This relation was found by Strouhal,> who 
rotated a stretched wire about an axis parallel to 
its length. The tones observed were harmonics of 
the wire. In the case to be considered here the 
wire does not vibrate and the frequency of the 
tone is continuously variable. 


DESCRIPTION OF APPARATUS AND METHOD 


The platinum wires were supported at one end 
in a micrometer, which permitted distances of 
0.0005 cm to be read. The efflux velocity of the 
gas was calculated from the relation U=Q/At 
where Q is the volume of gas passing through 
the orifice in a time ¢, and A is the area of the 
cross section of the orifice. The volume of the gas 
was measured by a standard 0.1 cu. ft. meter 
made by the American Meter Company. In 
general the time was taken for a complete 
revolution of the index, although the dial was 
divided into 100 divisions and thus could be 
read to 0.001 cu. ft. The meter was calibrated by 
the aspirator bottle method, after which the 
water level was adjusted to compensate for the 
water removed by the gas. A diagram of 
the apparatus is shown in Fig. 2. 

Between the gas meter and the gas main a 35 
liter reservoir was used to reduce the pressure 
fluctuations. The pressure was varied by means 
of a regulating valve placed between the reservoir 
and meter. As these experiments were conducted 


5 Strouhal, Ann. d. Physik 5, 216 (1878). 
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in a constant temperature room the error in 
volume measurements, due to temperature varia- 
tions, was very small, since the usual variation 
was less than 2°C. 

While velocities were not calculated from 
pressure data it was found convenient at times 
to obtain the velocities from a pressure velocity 
curve previously determined by direct measure- 
ment. The jet consisted of a brass pipe 12.5 cm 
long, with an inside diameter of 0.66 cm, into 
one end of which could be placed removable 
nozzles. All nozzles except the one marked (g) 
(glass) were made of brass. The design of the 
nozzles is shown in Fig. 3. The length of the taper 
L’ was 1.6 cm. The Jength of the channel, L—L’, 
was about 2 mm. The nozzle with diameter 
O=0.907 mm was made with a taper length of 
0.3 cm and channel length of 1.1 cm. 

The frequency of the note emitted by the wire 
was determined by obtaining zero beat with an 
audio oscillator. A General Radio oscillator 
type 377B was used. In practice it was found 
convenient to set the oscillator at a desired 
frequency and to vary the gas pressure on the 
jet to secure resonance. 


METHOD OF AMPLIFYING TONE 


It was found that the intensity of the tone 
emitted by the wire was greatly increased when 
a second wire was placed in the flame above the 
first, or sound initiating wire. Best results were 
obtained when the second wire was placed 
parallel to the first wire at a distance of about 
2 cm downstream and about halfway between 
the center of the burning jet and the visible 
boundary. By projecting the light from this 
wire, which was a bright red, onto a photo- 
electric cell further amplification was possible. 
An audio amplifier with a voltage amplification 
of 50,000 was used in connection with the photo- 
electric cell. 
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TONES FROM A WIRE 


A cathode-ray oscilloscope connected to the 
photoelectric cell through the amplifier showed 
that the oscillation was a sine wave and that 
the position of the second wire did not affect the 
frequency except to introduce, or to amplify, 
the first overtone when placed very close to the 
first, or sound initiating wire. 

When the light from the flame was projected 
onto the photoelectric cell little effect was pro- 
duced if the wire was less than about 1 cm from 
the orifice, while for greater distances a large 
output was obtained from the cell. When pro- 
jecting the light of the flame the second wire was 
removed from the flame. 


ACTION OF WIRE IN FLAME 


If a wire whose diameter is nearly that of the 
orifice is moved slowly across an ignited jet, 
having the proper efflux velocity, a critical 
point will be reached near the center of the jet 
at which a musical tone will be heard. The 
space in which the tone continues is of the order 
of one-tenth of the visible diameter of the flame 
and varies with the distance from the orifice and 
the orifice velocity of the gas. The frequency of 
the tone varies as the wire is moved across the jet, 
an effect which produces the largest source of 
error in this study. 

Since the width of this critical region decreases 
with the velocity the method used was to move 
the wire backwards and forwards across the 
center of the jet while reducing the velocity. 
The final setting was made for a clear tone near 
the center of the jet. When all measurements are 
carefully made the entire apparatus may be re- 
assembled and a reading repeated with an 
average variation of about two percent. Occa- 
sional variations of four or five percent were 


observed. 
CRITICAL VELOCITIES 


If after placing the wire in position in the flame 
the velocity is slowly decreased, by decreasing 
the pressure, a point will be reached where the 
tone ceases. By increasing the velocity again 
the tone starts almost exactly at the same 
velocity at which it stopped. The tone con- 
tinues at a velocity slightly less than that 
required to start it. In practice the average 
velocity at which slight fluctuations in gas 
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pressure would cause the tone to start and stop 
was taken to be the velocity of initiation U;. 

As the velocity is increased beyond the value 
for initiation the frequency increases to a critical 
point at which one or more subharmonics start. 
These subharmonics are usually accompanied by 
a characteristic change in the shape of the flame. 
If the velocity is increased beyond this point the 
changes in shape and the number of sub- 
harmonics become increasingly rapid until finally 
the separate tones merge into the characteristic 
“flare.” This critical velocity is called the flare 
velocity with wire and is designated by Uy. 
The results for the initiation and flare velocities 
are given in Table I and are shown in Fig. 4, 
where critical velocities are plotted against 
orifice diameter. In these tables D is the diameter 
of the wire, d is the distance of the wire from 
the orifice, and O is the diameter of the orifice. 
The data in Table I were obtained while using 
artificial gas. Table II shows the results of a 
similar experiment using natural gas from the. 
supply at Lexington, Kentucky. 

After the wire is removed from the flame a 
still further increase in velocity brings the jet 
into the condition in which it is sensitive to 
sound. A further increase in velocity causes the 
flame to break down into the familiar flaring or 
turbulent condition. The results for these critical 
velocities for sensitivity and flare are shown in 
Fig. 4. All data plotted in Fig. 4 were obtained 
using artificial gas. 

While the curves shown in Fig. 4 are not 
exactly hyperbolae this is very nearly the true 
relation for the initiation of the tone, using 
artificial gas. If it is assumed that the kinematic 
viscosity v is constant then OU/vr is constant, 
that is, this transition corresponds to a definite 


TABLE I. Initiation and flare velocity. 


D=0.0484 cm d=0.05 cm 


U; m/sec. 


7.15 
7.53 
8.80 
10.47 
13.90 


O mm 


1.779 
1.689 
1.420 
1.169 


Uy m/sec. 
18.14 
19.42 
23.50 
27.60 
34.45 


OU; cm?/sec. 


323 
330 
332 
325 
317 
307 


322 


OU; cm?/sec. 


127.0 
126.8 
125.0 
122.2 
126.2 


Average 
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Reynold’s number. Since the kinematic viscosity 
is a function of temperature it may be that each 
of these transitions correspond to a definite 
Reynold’s number. If the kinematic coefficient 
’ of viscosity were known for each of the gases 
used, this relation could be tested. Since the 
composition of the artificial gas used does not 
remain constant the value of the kinematic 
coefficient as calculated from the average com- 
position must be considered only an approxima- 
tion. The calculated value is 0.23 at 24°C. The 
corresponding value of v,, for methane, the main 
constituent of the natural gas used, was 0.167. 
For the same temperature the following relation 
is to be expected if the Reynold’s number for a 
given transition is constant: 


(OU)a/0.233 = (OU)m/0.167. 


From Tables I and II the products OU; for 
artificial and natural gas are, respectively, 126.3 
and 95.4. These values substituted into the above 
equation give 542 for artificial gas compared with 
565 for natural gas. 


INITIATION VELOCITY AS A FUNCTION OF ORIFICE 
DIAMETER AND DISTANCE OF WIRE 
FROM ORIFICE 


An examination of Fig. 5 shows that there is a 
linear relation between the distance of the wire 
from the orifice and the initiation velocity. 
When the wire is at distances greater than about 
one centimeter the equation which gives the 
relation between d and U; is 


d=mU;+K (2) 
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where m=dd/dU. Table III gives the values of 
m and the ratio O/m; for different orifice di- 
ameters. The fact that the quantity O/m; 
appears to be a constant is a very striking 
result. 

Using the average value 0.550 for O/m; the 
following equation is obtained, 


O/m;=O(dU/dd) =0.55. (3) 
Eq. (3) becomes, on integrating, 
OU;=0.55d+OU5. (4) 


While further data are needed on this point it 
appears that the d intercepts of the lines in 
Fig. 5 are constant. Using the average value of 
— 2.49 found in Table III Eq. (4) becomes, 


OU;=0.55d+1.37. (5) 


By use of Eq. (5) the velocity of initiation may 
be calculated for any orifice diameter and for any 
distance of the wire from the orifice, when using 
a wire whose diameter is 1.01 mm. The equation, 
of course, is limited to the linear portion of the 
curves shown in Fig. 5, and to the assumption 
that the d intercepts are constant. A similar 
equation appears to hold for the flare velocity, 
in fact, the initiating and flare lines for a given 
orifice are almost parallel, as is indicated in 
Table III. A flare line is shown in Fig. 10. The 
determination of the flare velocity becomes 
increasingly difficult as the distance of the wire 
from the orifice is increased. 






INITIATION VELOCITY AND DIAMETER OF WIRE 
AS A FUNCTION OF THE DISTANCE OF 
THE WIRE FROM THE ORIFICE 


In Fig. 6 the initiation velocities for wires of 
different diameters were plotted as a function of 
the distance from the orifice. The results of a 


TABLE II. Initiation and flare velocity—natural gas. 





D=0.084 cm d=0.05 cm 





O mm U; m/sec. | Us m/sec. | OU; cm?/sec. OW; OUfs 
1.779 6.20 9.44 111 82 168 
1.420 7.30 12.88 103 89.5 | 183 
1.169 7.60 14.70 89 82.2 | 171.5 
0.797 9.80 20.90 78.6 88 166 
Average 95.4 85.4 | 172.1 
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study of these curves are given in Table IV. 
In this table m,, is dd/dU, or the slope of the 
line. The intercepts on the velocity axis appear 
to be a constant for a given orifice, although the 
linear relation is not valid near the orifice. The 
shape of the curves near the orifice is shown 
more clearly in Fig. 5. The convergence of these 
initiating lines on the velocity axis indicates 
that near the orifice the velocity for the initiation 
of the tone is independent of the diameter of 
the wire. A further study of these critical 
velocities was made by placing the wires at a 
distance of 0.05 cm from the orifice. For this 
distance both the velocity of initiation and the 
flare velocity are nearly independent of the wire 
diameter. The initiating velocity tends to in- 
crease for wires smaller than about 0.4 mm. 

In addition to the common intercept indicated 
in Table IV another striking result is that the 
ratio D/m, is very nearly constant for a given 
orifice. These relations may be expressed by the 
equation, 


D/m,=D(dU/dd) =a constant, 
or DU;=kd+e. (6) 


For the particular orifice, O=1.169 mm. Eq. 
(6) becomes 
U;=0.415d/D+12.6. (7) 


The essential difference between Eqs. (4) and (6) 
is, that for different orifices the common inter- 
cept is on the d axis while for different wire 
diameters the common intercept is on the U or 
velocity axis. 


COMBINED EFFECT OF DIAMETER OF ORIFICE 
AND DIAMETER OF WIRE ON THE 
INITIATION VELOCITY 


The two equations, (5) and (7), may be com- 
bined into a single equation. Since the diameter 


TABLE III. Initiation and flare velocities vs. orifice diameters; 
D=1.01 mm. 




















Oc mi do (INTERCEPT) O/mi my (FLARE LINE) 
0.1779 | 0.313 — 2.23 0.568 0.29 
.1420 .275 — 2.68 .516 .26 
.1169 .210 — 2.57 557 18? 
.0907 .163 .558 
Average —2.49 .550 




















TONES FROM A WIRE IN A GAS JET 325 


cM 
w 








| ; 
INITIATION OF 5 


> 





ORIFICE IN 


=) 


» 





_ 





DISTANCE FROM 

















c :« ts «+ 2 2 oe Se 
Orrricr VeLociry M/see 


Fic. 5. 


of the wire used in obtaining Eq. (5) had the 
value D=0.101 cm it follows that for any orifice 
diameter O, the constant k in Eq. (6) is given 
by the relation 


(D/mw)o,=0.101/m4. (8) 


For example if the ratio D/m.,, is desired for an 
orifice whose diameter is 0.097 cm the value of 
m; for this orifice taken from Table III gives the 
ratio 0.101/0.163 or D/m,=0.620, which is in 
good agreement with the average value of 0.602 
given in Table IV. 

By substituting the corresponding value of m; 
from Eq. (3) into Eq. (8) the following equation 
results 


OD/m.=OD(dU/dd) =0.055. (9) 
Eq. (9) becomes, after integrating, 
ODU;=0.055d+C. (10) 


In evaluating C in the above equation it is to be 
observed that when d=0 then OU;=OUp=a 
constant, independent of D the diameter of the 
wire, with the exceptions mentioned. The average 
value of OU») was found to be 1.41, where U is 
measured in m/sec. and O in cm. (The value of 
the intercepts for the orifice whose diameter is 
0.0907 cm was not included in this average due 
to its construction.) With the substitution of the 
value of C, which is C=1.41D, Eq. (10) becomes 


OU;=0.055d/D+1.41. (11) 


This equation is surprisingly simple in view of 
the fact that the actual velocity and viscosity 
are both unknown functions of the distance d. 
There is some reason to expect that the con- 
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stants in the above equation would depend on 
the construction of the jet as well as they would 
depend on the nature of the gas. 


FREQUENCY OF THE TONE 


The frequency of the tone depends on the 
velocity of the gas, the diameter of the orifice, 
the diameter of the wire, and the distance of 
the wire from the orifice. It has been found that 
whatever the distance of the wire from the 
orifice, the frequency is a linear function of 
the orifice velocity of the gas. This feature is 
clearly shown in Fig. 7, where the frequency is 
plotted as ordinates against the orifice velocity, 
in meters per second, as abscissae. It is observed 
that the slope decreases with the orifice distance 
of the wire so that the equation of these lines is 
of the form 

N=f(d)U-—g(d), (12) 


where f(d) represents the variation of the slope 
with distance and g(d) stands for the variation 
of the velocity intercept with distance. The 
direct determinations of these functions is made 
difficult by the variation of the slope depending 
on the position of the wire with respect to the 
center of the jet. 

An examination of these curves shows that 
lower frequencies are observed at greater dis- 
tances from the orifice. While the entire range of 
frequencies between the initiation of the tone and 
the beginning of flaring is not shown in Fig. 7 
it is a general rule that the frequency range for 
any distance of the wire from the orifice is about 
two octaves. The most noticeable feature of 
Fig. 8 is the abrupt change in slope, which 
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appears at about the same frequency at the 
various distances. A study of these slopes is 
given in Table V. 

For O=1.169 mm and D=1.01 mm an average 
of four values of m’/m gave 1.70. For O=1.420 
and D=1.01 mm the average of five such ratios 
of m’/m was 1.45. While no definite conclusions 
can be drawn, it appears that when a change in 
the slope occurs the ratio m’/m is approximately 
constant. The cause of these breaks is unknown 
and their appearance is quite erratic. It was 
noted that when breaks occurred in a family of 
curves like those shown in Fig. 7 they appeared 
in consecutive curves. 


FREQUENCY OF TONE AS A FUNCTION OF 
VELOCITY AND DIAMETER OF WIRE 


When the wire is placed very near the orifice 
it may be assumed that the velocity of the gas 
passing the wire is equal to the efflux velocity. 
The frequency velocity curves for wires of 
different diameters placed at a distance of 0.05 
cm from the orifice are shown in Fig. 9. These 
curves resemble those obtained by Strouhal® for 
Aeolian tones from wire of different diameter, 
and they are almost identical with the curves 
obtained by Kohlrausch’ for jet tones from 
circular orifices. The results of a study of these 
and similar curves are given in Table VI. In this 
table the slope (dN/dU) is measured in vibra- 
tions per cm/sec. 

A comparison of the values of dN/dU for 
D=0.101 cm suggests that the value 0.736 in 
column II of Table VI is in erorr and the product 
D(dN/dUO) =0.074 was not included in the 
average. Data from wires smaller than 0.04 cm 
are not included in the averages given in Table 


TABLE IV. Initiation velocity and diameter of wire. 








O=1.169 mM O=0.097 mM 





U (intercept) 


22.5 
23.3 
22.9 


D/myw | U (intercept)| mw | D/my 


0.240 |0.420 0.159 | 0.636 
182] . .143 | .568 
.158 .104| .602 
.128 


Dcm Mw 


0.101 
.081 
.063 
.050 








Average 22.9 


6 Strouhal, Ann. d. Physik 5, 216 (1878). 
7 Kohlrausch, Ann. d. Physik 13, 545 (1881). 
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VI. Kohlrausch’? measured the frequency of the 
tone produced by a jet of air streaming into air 
and found the products O(dN/dU) to lie between 
0.04 and 0.06. Also, he found that these products 
decreased with the diameter of the orifice, the 
values being almost identical with those given 
for wires of the same diameter, as given in 
Table VI. 

These results tend to show that for the same 
efflux velocity the frequency of a tone from a 
wire of diameter D is almost the same as that 
for a jet of the same diameter. It must be noted, 
however, that no tones were observed without 
placing some obstacle in the stream. Even when 
the light from a burning jet was focused onto a 
photoelectric cell using an amplification of 50,000 
no tone was heard. 

The equation for these lines is the same as 
that given by Kohlrausch for jet tones, 


N=k/D(U—U)), (13) 


where O is replaced by D. In this case Uo, the 
velocity intercept, is nearly a constant for the 
lines in Fig. 9. However, other data tends to 
show that Uy») becomes larger for very small 
wires. The approximate equation for the lines 
shown in Fig. 9 is 


N=(0.048/D) U—660, (14) 


where U is measured in cm/sec. and D is meas- 
ured in cm. 

While not enough data were taken to show the 
general relation of the quantity D(dN/dU) to the 
diameter of the orifice, it was found that for 
D=0.0484 cm at a distance d=0.05 cm the same 
velocities produced corresponding frequencies 
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for the four jets whose diameters were 0.797, 
0.907, 1.009, and 1.169 mm. The common 
frequency range for these jets was 1000-2000. 
For larger wires other data tended to show that 
the quantity D(dN/dU) is slightly greater for 
smaller orifice diameters. A difficulty encountered 
in making this comparison is that a given wire 
diameter gives consistent results for only a very 
small range of orifice diameters. 


APPROXIMATE RELATION BETWEEN FREQUENCY, 
VELOCITY, AND DISTANCE 


A simple method for studying the relation of 
frequency to the distance of the wire from the 
orifice is to move the wire up or down in the jet 
and to determine the frequency at various 
distances while keeping the orifice velocity con- 
stant. The difficulty here is that the new setting 
required for each reading introduces new errors 
of position with respect to the center of the 
stream. While no simple relation was found, a 
rough approximation is Nd=a constant. 

In an attempt to overcome the difficulties 
mentioned above, a frequency velocity curve was 


TABLE V. Slopes of frequency velocity curves from Fig. 8. 





O =1.779 MM D =2.50 MM 


m’ (upper slope) 


35.7 
37.0 
31.2 
29.9 
24.3 
19.0 


m (lower slope) 


24.2 
22.2 
19.2 
15.6 
15.4 
10.3 


Average 
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determined at various distances from the orifice, 
as shown in Fig. 7. From these curves distance 
velocity curves at constant frequencies were 
constructed for four orifice diameters. For each 
_ orifice it is obvious that all possible constant 

frequency curves must be bounded by the 
initiation and flare lines for that orifice, as 
shown in Fig. 10. 

In studying these curves log U was plotted 
against log d, as shown in Fig. 11. For values of 
d greater than about 1 cm the resulting curves 
are linear. In this region the relation between 
velocity and distance for a constant frequency 
may be expressed by 


Un = Kd", (15) 


where Uy is the orifice velocity required to give a 
frequency JN at a distance d. In this same manner 
four equations for four different orifice diameters 
were determined. Values of ” varied from 0.395 
to 0.629. In most cases the value was not far 
from 0.5. It was found that when different 
orifice diameters produced the same frequency 
the constants in Eq. (15) were related to the 
orifice diameter O by nO!=a constant and by 
KO'=a constant. By substituting for K its 
value in terms of N an approximate relation was 
found to be of the form 
N=U/(Bd')+(100—A/B), (16) 
where A and B are constants for a given orifice 
diameter. The slopes of the lines similar to those 
shown in Fig. 7, should then be given by dN/dU 
=1/Bd'. While this equation shows the general 
nature of the change of slope with distance, 
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the actual decrease of the slope with distance is 
less than that predicted by this formula. 


NATURE OF TONE 


The variation of the frequency of the tone as 
a function of the distance, at constant orifice 
velocity, is somewhat similar to the relation 
which exists in the case of edge tones. In the 
edge tone phenomena there is no interaction 
between the vortices after they strike the edge. 
If the wire acted as an edge in the jet, then a 
thin sheet or wedge placed behind it so as to 
form a continuous wedge with the wire as the 
apex, should not prevent the production of a 
tone. 

To test this condition a metal sheet 0.025 cm 
thick and 1.9 cm wide was placed behind the 
wire. It was found that when the sheet was 
slowly moved toward the wire or away from it 
in the direction of the stream the tone ceased 
and started again at a definite distance from the 
wire. This critical distance varies with the 
velocity and the diameter of the wire and perhaps 
with the diameter of the orifice, though this 
relation was not examined. A graph of the 
relation between the distance for the cessation 
of the tone and the efflux velocity for a given 
wire and orifice diameter is shown in Fig. 12. 
The equation of this particular curve is given by 


2=19/U"1, (17) 


where z, the edge distance, is measured in 
millimeters and U is measured in meters per 
second. 

At a constant velocity this critical distance 
varies with the diameter of the wire. The values 
of z for various wire diameters at a constant 
velocity of 22 meters per second are given in 


TABLE VI. Diameter of wires and slopes of N-U curves. 
































O=1.009 mm (g) d=0.05 cm O=1.169 mu O=1.169 mm d=0.247 cm 
Dem |dN/dU| D(dN/dU)|dN/dU| D(dN/dU)| D_ |\dN/dU| D(dN/dU) 
0.101 0.736 0.074 0.538 0.054 0.1010 | 0.542 0.055 
.082 734 .060 558 .046 .0814| .604 .049 
.072 841 .046 .0635 | .755 .048 
.050 897 .045 823 041 .0509| .838 .043 
048 1.010 .049 
040 1.025 041 973 .039 .0407 | 1.050 .043 
.032 1.22 .038 .0330 | 1.057 .033 
025 .0254 | 1.305 033 
Average .048 045 .047 
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TONES FROM A WIRE 


Table VII. If the sheet is placed in contact with 
the wire and then moved in a horizontal plane 
the tone starts again at a definite distance X 
from the center of the wire. These values are 
given in Table VII. 

An examination of Table VII shows that for 
wire diameters larger than 0.40 mm the relation 
between z and D may be expressed by 


D‘/z=a constant. (18) 


Under the corditions given in Table VII the 
value of this constant is 1.73, hence Eq. (18) 
may be written as z=0.576D}. The ratio X/D 
increases as D decreases. The abrupt increase in 
the value of z as the diameter of the wire ap- 
proaches the thickness of the sheet is very 
apparent. 

These results tend to show that the tone is 
produced by an oscillation in the stream immedi- 
ately behind the wire, as is the case in the 
Aeolian tone phenomenon. In general, the in- 
vestigators studying Aeolian tones in gases have 
had to depend on the tone produced by the wire, 
which was in resonance with the shedding of the 
vortices. In these experiments the wire does not 
vibrate and the tones studied were observed 
only when the gas was burning. 

When the wire was heated to a bright red 
glow by an electric current a very faint tone 
could be heard. It is believed that the oscillations 
are present in the unignited gas and that the 
tone in an ignited gas is due to the miniature 
explosions produced by the formation of vortices 
behind the wire. The formation of an eddy 
behind an obstacle in the flame would certainly 
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bring about a mixture more suitable for com- 
bustion. It appears that the action of the 
amplifying wire is due partly to increased 
mixing of the gas as an eddy strikes it and 
partly due to the change in temperature of the 
wire produced by this agitation. 

The action of the amplifying wire suggests 
that the tone produced by the first, or sound 
initiating, wire may be due to its amplifying 
action on the vortices emitted by the orifice. 
However, in moving the amplifying wire in the 
stream there is no change of frequency with the 
distance from the source, while there is a change 
of frequency if the distance of the first wire from 
the orifice is changed. 
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DISCUSSION OF RESULTS 


Several points could be cleared up if the actual 
velocity of the gas passing a wire at a distance d 
from the orifice could be determined accurately. 
If it be assumed that the frequency for a given 
wire and orifice depends only on the velocity of 
the gas as it passes the wire, then the factor d” in 
Eq. (15) should give the rate of decrease of 
actual velocity with distance. It may be of 
value to consider the nature of flow in a jet. 

For the small pressures used it may be 
assumed that the total pressure of the gas 
remains constant. It follows that the decrease in 
velocity is due to frictional stresses alone. While 
fresh masses of air will be dragged into the jet 
the momentum I of the stream must be constant, 
or 


I=f fdrU*dA =a constant, (19) 


where 7 is the density of the gas and A is the 
area of the cross section of the stream. From 
this relation 


r,;U7Ai=fe U2A 2. (20) 


Observation shows that a jet is conical in shape, 
so that 


A,/A2=d;2/d.?. (21) 
With this substitution Eq. (20) becomes 
U2/U1= (di/d2) (11/12). (22) 


TABLE VII. Variation of edge distance with wire diameter. 














O =1.169 MM U =22.0 M/sEc. d =0.05 cm 
| 
Dmm zmm p} D)/z X mm X/D 
1.01 0.576 1.004 1.74 1.03 1.020 
.826 559 .908 1.63 .872 1.055 
.76 .510 .872 1.71 
50 | .397 707 | 1.78 595 | 1.190 
40 .340 .632 1.80 495 1.235 
32 750 .566 0.75 530 1.655 
25 760 .500 0.66 580 2.320 
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For an unignited jet 7;=72 and the velocity 
decreases inversely as the first power of the 
distance, or Eq. (22) may be written, 


U' =k/d. (23) 


Taking the approximate value of m in Eq. (15) 
to be 3 then Eq. (16) indicates that for an 
ignited jet the rate of decrease of velocity with 
distance may be expressed by 


U'=k/d', (24) 


a result which is obtained from Eq. (22) if the 
density of an ignited jet varies inversely as the 
distance. It was noticed that the temperature of 
the wire, as judged from its color, appeared to 
increase with the distance from the orifice in the 
region studied. 

It is generally considered that viscosity has a 
very small effect on the frequency of an Aeolian 
tone, but that the initiation of the tone is very 
definitely a function of viscosity. An examination 
of Fig. 10 shows that the increase of the orifice 
velocity with distance is less for the initiation of 
the tone than for the maintenance of a given 
frequency. For the initiation of vortices behind a 
cylinder Richardson® gives an average value of 
the quantity UD/v the number 30. This is to be 
compared with the values 542 and 565 found for 
artificial gas and natural gas, respectively. Unless 
this can be accounted for on the basis of increased 
viscosity due to the increase in temperature then 
it must be concluded that higher velocities are 
required for the initiation of a tone with a jet 
than for the ordinary Aeolian tone. 

In conclusion, the author wishes to thank Dr. 
A. L. Foley, the director of this research, for his 
helpful comments and suggestions during the 
progress of this work. The author wishes also to 
thank Professor W. S. Webb, of the department 
of physics at the University of Kentucky, for the 
use of apparatus and facilities required for 
making measurements with natural gas. 


8 F, G. Richardson, Sound, p. 153. 
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Sound Absorption and Attenuation by the Flue Method 


HAWLEY O. TAYLOR AND CHALMERS W. SHERWIN 
Wheaton College, Wheaton, Illinois 


(Received January 21, 1938) 


HE method of measuring the absorption of 

sound as herein outlined is the same as that 
developed in a former paper! in which the sound 
flue was used. The term “‘flue”’ was originally 
proposed in order to suggest a pipe of rather 
larger cross-sectional area than is usual in the 
case of tubes for the study of sound interference. 
Several writers of previous papers had noted the 
effect of absorption on the interference pattern in 
tubes, but as absorption was not of primary 
interest the coefficient was not developed. In the 
work of Eckhardt and Chrisler,? the effect of 
attenuation in the flue is taken into account as a 
damping factor in the wave equations, but as 
“exact solutions of the equations do not seem 
practicable’”’ the attenuation was not evaluated. 
An approximate solution however “appears to 
be adequate,” from which the following expres- 
sion for the coefficient of absorption, a, is 
developed : 


m$—ny)+ (n2!—ny3)/27? 
on -[ | 


m+n yi— (n2i—n,')/2 


where m is the intensity of the first maximum of 
the interference pattern in the sound flue 
(measuring from the material closing the end), 
and m, and m, are the intensities of the first two 
minima of this pattern. When attenuation is 
practically absent, this expression reduces to the 
simpler form developed in the former paper! to 
which reference has been made: 


m$—n*7? 4 
a=1-| |-— ee (2) 
m\+n} (m/n)'+(n/m)'+2 





The purpose of the present paper is to show 
how the sound flue may be used to determine the 
fraction of sound attenuated in a unit distance 
along the axis of the flue by the material lining 


1H. O. Taylor, ‘A Direct Method of Finding the Value 
of Materials as Sound Absorbers,” Phys. Rev. 2, 11 (1913). 

2E, A. Eckhardt and V. L. Chrisler, ‘‘Transmission and 
Absorption of Sound by Some Building Materials,’’ Sci. 
Papers of National Bureau of Standards, No. 526, April, 
1926. 


its sides, and to translate the change produced in 
sound intensity given by this coefficient into 
decibels per square foot of the material when so 
placed in the flue. 


THE SOUND FLUE 


The research utilized several sound flues, three 
of which had square cross sections with inside 
measurements of 5, 8 and 11 inches on a side 
and about 6 feet long. Closing one end of each 
flue was a cover clamped to a flange at the end. 
Sound absorbing material was attached to the 
inside of this cover and also as lining to the side 
walls of the flue. Hair felt, 1 inch thick, and 
Acousti-Celotex C1, C3 and CB were used for 
these linings. The source of sound at the open end 
of the flue was usually a telephone receiver 
connected to an audio-oscillator. A Magnavox 
connected to the oscillator was sometimes used 
for the frequency of 256~ because of the pure 
tone produced. The tone used must be practically 
free from harmonics for dependable results, also 
all cracks in the flue and end cover must be 
stopped with some such material as wax. 

The explorer for measuring the intensity of 
sound along the axis of the flue was a brass tube 
} inch in diameter and about 8 feet long, open at 
one end and connected to the cap of a telephone 
receiver at the other end. The intensity of the 
sound at the open end of the tube excited the 
telephone receiver, and the voltage generated in 
the receiver was measured with an amplifier and 
a suitable rectifier and ammeter. The exploring 
equipment was calibrated by comparison with a 
Rayleigh disk over the range of intensities 
encountered in the flue measurements. 


COMPLETE ATTENUATION IN THE SOUND FLUE 


If the flue is small in cross section and lined 
with material of high absorbing power, the 
standing wave interference pattern may be 
inappreciable, leaving only diminishing intensi- 
ties as indicated by the explorer when it was 
advanced from the open end of the flue inward. 
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In this case there was practically no sound wave 
reflected from the material at the closed end of 
the flue—the sound energy was practically 
consumed before it reached this end. 

A case of this type is shown in Fig. 1 which was 
obtained at a frequency of 1024~ for a 5-inch 
flue lined with }3-inch C1 Celotex. This curve 
shows practically the same fractional reduction in 
intensity, or in other words, a uniform attenuation 
per unit length throughout the length of the flue, 
amounting to 0.22 or about 63 decibels per foot. 


EFFECT OF ATTENUATION ON STANDING WAVES 


If the sound energy is not too highly thinned 
out before the waves reach the closed end of the 
flue, then a standing wave pattern results, the 
regions of maximum or minimum intensity being 
one-half wave-length apart. Suppose that a 
sound wave p proceeds down the flue in the 
direction of the arrow (Fig. 2). Let its intensity at 
position 0 be 7. (By intensity may be understood 
the rate at which sound energy would flow away 
if opportunity were given, as through the ex- 
ploring tube.) Then at positions 1, 2 and 3 (the 
reflecting surface), its intensity will be by, 6,2 
and 6,*7, respectively, where 0}, is the attenuation 
constant for a distance along the flue equal to a 
quarter wave-length. At the reflecting material 
closing the end of the flue, the effective surface of 
which seems to be embedded beneath the actual 
surface, sound energy will be absorbed, the 
fraction left being given by p. The reflected wave 
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W. SHERWIN 
r will return through the flue, its intensities at 
positions 3, 2, 1 and 0 being b;*zp, b,*ip, b;5ip and 
b; ip, respectively. The progressive and reflected 
waves combining will produce maximum and 
minimum regions of sound intensities, as shown 
in Fig. 2, the first and second minima from the 
reflecting surface having intensities of (say) n, 
and me, respectively, and the intervening maxi- 
mum having an intensity of m. The distance 
along the flue from ,; to m (that is, from 2 to 1) 
is a quarter wave-length \/4, as is also the dis- 
tance between each of the neighboring positions 
marked. Since intensity measurements made at 
these positions are used to determine the attenu- 
ation of the material lining the flue, it will be con- 
venient to use a quarter-wave-length as the unit 
distance for computing the attenuation constant. 

From the theory of sound waves interfering in 
the flue,! the maximum intensity m is pro- 
portional to the square of the sum of the ampli- 
tudes of the progressive and reflected waves; and 
the minimum intensity 1 is proportional to the 
square of their difference. Then taking the 
intensity of sound as proportional to the square 
of the amplitude of the waves: 


my (bri)!-+(br8ip)! _ (b1)!+ ('p)! 
-| =——_—__---—______=C€ (9) 


Ne (2) te (b,°ip)} 1— (b,°p)? 


and 
my}! (byt) §+(b1*ip)! 1+ (di‘p)} 
(b:*)*— (bi %ip)} 


(b1)'— (dip)? 
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SOUND ABSORPTION 
By eliminating p' and the factor (6;,7+1) from 
Eqs. (3) and (4), the following quadratic equation 
in (b;)! results: 


b,CD+(b;)}(D—C)—CD=0, 


from which 
: pecs ee 
isi 2CD 


n,i—no! ny3—n2)\ 2 142 
-| +((" ) +1) |. (5) 
2m} 2m} 

The data of Fig. 2 were taken from the 8 inch 
flue lined with Acousti-Celotex C3, 7 inch thick, 
and the end of the flue was closed with the same 
material. Sound of frequency 256~ was used. 
The intensities were: ;=13.8 at 26 centimeters 
from the closed end; m=450 at 55 cm; and 
n2=45 at 78 cm. The half-wave-length is seen to 
be 78-26 or 52 cm, that is, 104 cm for the wave- 
length. In the open air, the wave-length for this 
frequency would be about 133 cm. This shorten- 
ing of wave-length for sound in flues lined with 
absorbing material has already been noted.' By 
substituting values of ,, m2. and m in (5) we 
find 0.869 for the value of 0,, the fractional 
reduction of sound energy per quarter-wave- 
length. 








COEFFICIENT OF SOUND ABSORPTION 


The coefficient of sound absorption, a, is equal 
to 1—p. Using the expressions for p obtained 
from Eqs. (3) and (4), we get from (3): 


C—(b:)} 


a eee 
. aia ou ' 
and from (7): 


: , D(b;)'-1 ] (7) 
oP Ly) (by! D) I 


If there is practically no attenuation of sound 
in the flue, then b;=1, and C=D=(m/n)}. These 
values placed in Eqs. (6) and (7) reduce them to 
formula (2). On the same assumption, m; and m2 
of Eq. (1) will become equal, and the expression 
for a will also reduce to formula (2). 

Substituting the values indicated on Fig. 2, we 
get a=0.356 from (6), 0.354 from (7) and 0.348 
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from (1). These values are all in good agreement. 
On the other hand, formula (2), in which the 
attenuation is neglected, yields a value equal to 


0.634. This is about 80 percent too large. 


SoUND ATTENUATION PER SQUARE Foot 


The attenuation constant 6 for a square foot of 
lining material in position in the flue would be 
determined for one foot along the flue if the 
perimeter of the flue measured just one foot on 
the inside. For perimeters less than one foot, } 
must be found for a greater distance along the 
flue, and for greater perimeters a correspondingly 
less distance. If f is the inside width of one side 
of the flue, then a distance along the flue of 
1/(4f) would be required to make one square 
foot of flue lining, where f is measured in feet. 
The attenuation coefficient B for this length 
along the flue would be 


B=buen, (8) 
The relation between } and ), is 
b=b,4!, 


since a quarter-wave-length is contained in a 
length of one foot 4/d times, where X is expressed 
in feet. Then B in terms of }, is: 


Bad, (9) 


Eq. (9) gives the attenuation coefficient per 
square foot of material in position in the flue in 
terms of the attenuation constant b, determined 
from measurements of maximum and minimum 
intensities (m, m; and m2) of standing waves in the 
flue. 


ATTENUATION IN DECIBELS 


The decibel is a unit of change in sound 
intensity. If the intensity changes, because of 
attenuation, from i; to i2, then, expressed in 
decibels, db, this would be 


db = — 10 log (i2/11). (10) 


If this change occurred when sound swept over 
one square foot of material, then Eqs. (9) and 
(10) may be combined to give it in decibels. 
That is, 


log B=(1/(fd)) log b1=log (42/71), 
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which, when multiplied by — 10, gives the change 
in sound intensity in decibels: 


(—10/(fdA)) log bs = —10 log (i2/i1:) =db per sq. ft. 


Using the data of Fig. 2, db for material C3 at 
256~ would be 


db =(—10/3.04) log 0.869 =0.203. 


If the attenuation 0 is taken for one foot as the 
unit, then the expression for decibels per square 
foot would, from (8), be 


log B=(1/(4f)) log 6; 
db = (—10/(4f)) log b. 


If decibels per square foot are plotted as 
ordinates against f? (the inside area of the square 
cross section of the flue) as abscissas, the resulting 
graph (Fig. 3) approximates a_ rectangular 
hyperbola. It is reasonable, as the graph indi- 
cates, that, as the area of the flue diminishes, the 
attenuation in decibels per square foot of flue 
lining would approach infinity; and, as the flue 
area approaches infinity, the attenuation would 
diminish to the constant, dbo. The equation of 
the curve is: (db—dbo)f?=q (a constant). 


from which 
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Plotting decibels as ordinates against the recipro- 
cal of f? would give approximately a straight line. 
The equation of this line is 


db=q(1/f*) +dbo. (11) 


Table I gives some average results from many 
observations made on Acousti-Celotex C1 and 
C3. \ is the wave-length of the sound in feet; f is 
the inside width of the square flue in feet ; }, is the 
attenuation of sound per quarter-wave-length 
along the flue; and db is the attenuation of 
sound per square foot of flue lining. 

Db and 1/f? in the table were used for the 
graphs of Fig. 4. Knowing the slope and y- 
intercept for a material, its attenuation per 
square foot when used as a lining in a given size 


of flue may be calculated by means of Eq. (11). 


The given flue area in square feet is substituted 
for f?, and the equation is solved for db, the 
required attenuation per square foot. The length 
of flue required to make a lining area of one 
square foot would be 1/(4f), where 4f is the inside 
perimeter of the flue. This is the length of flue 
that would produce the attenuation found. 

Let us suppose that a flue of square cross 
section is lined with Acousti-Celotex C3. What 
flue area would be required to suppress 10 
decibels per foot length of flue for a sound of 
frequency 512~? 

To apply formula (11), decibels per square foot 
must be found and substituted in the formula for 
db. If the required area of the flue is f? square 
feet, then 4f is the length of its perimeter, and 
1/(4f) is the length along the flue required to 
give the C3 lining an area of one square foot. 
Since the attenuation in decibels is proportional 
to the length of flue, then, if there are to be 10 
decibels drop in sound intensity per foot, for a 


TABLE I. 





11-INCH FLUE 




















~ r f 1/f? bh db 1/f? bi db f 1/f2 hi db 
Cl 256 | 3.94 | 0.825 | 1.47 | 0.86 | 0.201 | 0.575 | 3.03 | 0.77 | 0.501 | 0.325 | 9.47 | 0.81 | 0.715 
C3 760 | 1.73. | .62 .692 522 | 3.67 73 .665 272 | 13.5 58 | 2.208 
Cl 512 | 1.97 825 | 1.47 85 434 575 | 3.03 81 .808 325 | 9.47 67 | 2.713 
C3 760 | 1.73 75 835 522 | 3.67 59 | 2.223 272 | 13.5 46 | 6.30 
C1 1024 985 | .825 | 1.47 84 932 1.246 325 | 9.47 68 | 5.23 


1.73 765 








8-INCH FLUE 5-INCH FLUE 















length of 1/(4f) feet, the drop in decibels would 
be 10/(4f)—which is the db of formula (11). 
From the graph (Fig. 4) for C3 512~, dbo=0.6 
and g=16/39. Therefore formula (11) becomes: 
10/(4f) =16/(39f2)+0.6, from which f=0.171 
feet or about 2} inches. The quadratic equation 
gives also a second and larger value for f. This is 
true also when the formula is applied for the 
value of f using the data of Table I, a second 
value is obtained beside that given in the table, 
indicating two sizes of flue for the same attenu- 
ation per foot along the flue. Data on flues 
having these larger areas are not yet available so 
that the discussion of this paragraph is in the 
nature of a progress report. 

The authors acknowledge with gratitude the 
valuable assistance of students who have partici- 
pated in the experimental work and of others who 
have shown interest in the work, particularly Mr. 
Jacob Mazer who furnished the audiofrequency 
oscillator and Dr. B. E. Eisenhour who loaned 
the Magnavox and other equipment. 


SyMBOLs USED 


a, coefficient of absorption of sound reflected from a 
material. 

B, attenuation constant per square foot of sound flue 
lining. 

b, attenuation constant per foot along sound flue. 

bi, attenuation constant per quarter-wave-length along 


sound flue. 

C, (m/nz)'. 

D, (m/n,)4. 

db, drop of sound intensity in decibels per square foot of 
flue lining. 


dbo, drop of intensity in decibels per square foot of material 
lining a flue of very large cross-sectional area. 
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f, inside width of square section of sound flue. 

1, t1, 22, intensity of sound under various conditions. 

A, length of sound wave. 

m, maximum sound intensity of standing wave in the flue. 

n, 1, M2, Minimum sound intensity of standing wave in 
the flue. 

bp, amplitude of progressive sound wave in the flue. 

aw, number of radians in 180° (3.1416). 

q, product of db drop per sq. ft. of flue lining and area of 
flue section. 

r, amplitude of reflected sound wave in the flue. 

p, coefficient of reflection of sound. 

~, cycles or complete vibrations per second. 





th 

‘i 

i 

: 
4 








APRIL, 1938 | 


S. A. VOLUME 9 


The M.K.S. System of Units Applied to Electroacoustics 


A. E, KENNELLY AND JACKSON H. Cook 
Harvard University and The Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received August 20, 1937) 


PURPOSE AND SCOPE 


HE subject of electroacoustics deals with 
the transformation of electric energy to 
sound energy and vice versa by means of various 
electroacoustic transducers—microphones and 
loudspeakers. The quantitative investigation of 
such transfers requires a definite system of units 
for expressing the magnitudes involved. At pres- 
ent two sets of units are commonly employed; 
electric quantities are expressed in the inter- 
national series of practical electric units, the 
volt, ohm, ampere, joule, etc.; while mechanic 
quantities are stated in the centimeter, gram, 
second (c.g.s.), system. These two sets of units 
are not homogeneous, one with the other, and 
their coincident use in electroacoustic work 
makes necessary the employment of unwieldy 
transformation coefficients. The joule, unit of 
energy in the practical series, is equivalent to 
10’ ergs of the c.g.s. system; the practical unit 
of power, the watt, is equal to 10’ ergs/second, etc. 
A great saving in labor on the part of electro- 
acoustic workers would be effected if a system 
of mechanic units, homogeneous with the prac- 
tical electric unit series, were employed. Such a 
system, the Giorgi meter, kilogram, second 
(m.k.s.) system of units, is available for use and 
its recent adoption by the I.E.C. (International 
Electrotechnical Commission) endows it with 
sufficient importance to warrant its immediate 
application.! It is the purpose of this paper to 
express commonly used acoustic quantities in 
the m.k.s. system, to note their magnitudes 
relative to the corresponding c.g.s. units, and 
to point out the advantages thereby obtained. 


Acoustic UNITs 


The unit of length in the m.k.s. system is the 
meter (m) (=10® cm), the unit of mass is the 


1A, E. Kennelly, “I.E.C. Adopts M.K.S. System of 
Units,” Elec. Eng. Dec. (1935). I.E.C. Document R.M. 
118. Minutes of the I.E.C., E.M.M.U. Meeting, 
Scheveningen-Brussels, June (1935). Gino Sacerdote, 
“L’Applicazione Delle Unita M.K.S. Elettromagnetiche 
(Giorgi) Nel Campo Dell’Elettroacoustica,” Alta Fre- 
quenza 9, 570 (1936). 
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kilogram (kg) (=10* gm). Time is measured in 
terms of the mean solar second (s) in both 
systems. 

Sound wave-length is measured in meters (m) 
(=10? cm). The wave-length of a 1000 c.p.s. 
sound wave in air, at ordinary barometric pres- 
sure intensity and temperature, is 0.344 (m). 

The velocity of propagation of a sound wave 
is stated in meters/second (m/s) (=10? cm/s). 
A free sound wave in air, under normal condi- 
tions, (20°C, 1 standard atmosphere) is propa- 
gated at a velocity, 344 (m/s). (This unit is 
almost universally used at the present time). 
Particle velocity is likewise expressed in (m/s). 

Sound pressure (intensity) is measured in new- 
tons/square meter (newtons/m?) (=10 barye or 
dynes /cm?) (Newton is the tentative name for the 
m.k.s. unit of force, one joule/meter, equals 105 
dynes). The present c.g.s. reference pressure (0.2) 
millibaryed) is 2X10-> newton/m? or 20 micro- 
newtons /m?. Normal barometric pressure is 1.0132 
X 105 newtons/m?=0.10132 mega-newton/m?. 

The unit of sound energy flow is the watt 
(=107 ergs/s). For an effective sound pressure, 
p (newtons/m?) in a plane or spherical free wave 
propagated with a velocity, c (m/s), in a gas of 
density, p (kg/m*), the sound energy flow, P 
(watts), through an area, a (m?”), is: 


P=(p'a/pc) cos 0 watts, 


where @ is the angle of incidence of the sound 
wave on the area, a. 

Sound intensity, the sound energy passing 
through unit area in unit time, is expressed in 
watts /square meter (watts/m*) (=10* ergs/s-cm?*). 
For an effective pressure, p (newtons/m?), in a 
plane or spherical free wave propagated with a 
velocity, c (m/s), in a gas of density, p (kg/m’), 
the sound intensity, I (watts/m?*) in the direction 
of propagation may be written: 


I=(p2/pc) watts/m?. 


Sound energy density has the units of joules/ 
cubic meter (joule/m*) (=10*! ergs/cm*). 
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Acoustic IMPEDANCE 


The subject of acoustic impedance is confused 
by the existence in the literature of conflicting 
definitions of this important quantity. Crandall, 
Kennelly, Morse, and others define the acoustic 
impedance of a sound medium to a surface lying 
in a wave front as the product of the surface 
area by the ratio of sound pressure to particle 
velocity in the medium. 


Z=(p/q)a. 


From this definition it follows logically that the 
motion of a given piston in a sound medium is 
impeded more than is the motion of asimilar piston 
having a smaller surface area. Furthermore the 
dimensions of acoustic impedance according to 
this definition, are those of a force divided by a 
velocity, identical with the dimensions of me- 
chanic impedance. 


Z=(p/q)a, 
[Force/Area | 
[Velocity ] 
Z=([Force ]|/[ Velocity ]. 


[Area ], 


Therefore in determining the behavior of a me- 
chanic system in a sound medium (a vibrating 
piston in air) the acoustic impedance offered by 
the medium introduces a directly additive modi- 
fication of the mechanic impedance of the system. 
In short, by this definition, acoustic impedance 
is mechanic impedance.” 

The fundamental dimensions of this kind of 
acoustic impedance are a mass divided by a time. 


’ [Force] [MLT~] 


= $$ = ———-= [MT]. 
[Velocity] [LT] 
The unit in the m.k.s. system is the kilogram/ 
second (kg/s) (=10* gm/s). It would be advan- 
tageous if the name mechanic ohm should be 
used for the unit of mechanic and acoustic im- 
pedance having this magnitude. 
The tentative “American Standards for Acous- 
tical Terminology”’ define the acoustic impedance 
of a sound medium on a given surface lying in a 


_ 7A. E. Kennelly and J. H. Cook, “Mechanic Impedance 
in M.K.S. Units.” Elec. Eng. 56, 1062-63 (1937). 


wave front as ‘‘the complex quotient of the sound 
pressure (force per unit area) on that surface by 
the flux (volume velocity or linear velocity multi- 
plied by the area) through the surface.’” 


Z=(p/q)(1/a). 


According to this definition the motion of a 
given piston in a sound medium would be less 
impeded than would be the motion of a piston 
having a smaller surface area. However a defini- 
tion of an acoustic impedance of this type is of 
use in problems in which the actual flow of sound 
energy through orifices is dealt with and a con- 
cept of impedance somewhat analogous to the 
resistance of a conductor to electric current flow 
is needed.‘ 

The ambiguity introduced by the above dual- 
ity of definition does not affect the concept of 
“specific acoustic impedance,’’ defined for a plane 
sound wave as the ratio of the sound pressure to 
the particle velocity in the sound medium. 


Z.=p/¢. 


The unit in the m.k.s. system is the kilogram/ 
second square meter (kg/s-m?) (=10- gm/s-cm? 
(Suggested name: mechanic ohm per square meter). 

For the case of a plane sound wave propagated 
with a velocity, c (m/s), in a sound medium of 
density, p (kg/m), the specific acoustic imped- 
ance reduces very simply to: 


Z.=pc kg/s-m?. 
The density of air under normal pressure at 20°C 


’ American Tentative Standard—Acoustical Termin- 
ology 24.1 (1936). American Standards Association. 

‘It is evident that these two definitions of acoustic 
impedance differ because they deal with two entirely 
different concepts; one the impedance offered to a plate 
vibrating in a sound medium, the other the impedance to 
the flow of fluid offered by an orifice in a sound medium. 
These two concepts might best be discriminated between 
by the use of different unit names such as piston impedance 
and orfice impedance, Z)» and Za, respectively. If this 
latter definition (Zz) is used exclusively for ‘‘acoustic 
impedance” it is necessary, when studying the motion 
of a vibrating piston in a sound medium, to multiply the 
acoustic impedance of the system by the square of the 
piston area in order that the contribution of the sound 
medium to the mechanic impedance of the system may 
be determined. The fundamental dimensions of acoustic 
impedance as defined in the “Tentative Standards” are 
[ML~“T~*}. 

_ £1 [Mer yir) 2 os 
ee my) a) 
The unit in the m.k.s. system is the kilogram/second- 
meter-to-the-fourth-power (kg/s-m*) (=10-> gm/s-cm*). 
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is 1.205 kg/m’; the velocity of propagation of a 
sound wave in air under these conditions is 
344.0 m/s. Thus: 


Z.=415 mech. ohms/m? (kg/s-m?) for air at 
20°C, 1 atmosphere. 


ARCHITECTURAL ACOUSTICS 


Reverberation time is measured in terms of 
the second (s), a unit common to both systems. 
Sabine’s equation for the reverberation time (T), 
in which the energy density is reduced by a 
factor 10°, is 


0.16 V 
penne: i, 


aa 





where V is the volume (m*) of a room having a 
total surface area (walls, floor and ceiling) (m7), 
with an average absorption coefficient, @. 

The steady state energy density, E(joules/m*), 
in this room when a sound source, P (watts), is 
operating in it is: 


E=— joules/m’. 
caa 


The coefficients of viscosity and of kinematic 
viscosity of a sound medium are used in studies 
of absorbing materials and of sound transmission. 
The coefficient of viscosity is measured in units 
of kilograms/meter-second (kg/m-s) (=10! gm/ 
cm-s). The coefficient of kinematic viscosity has 
the unit, square meter/second (m?/s) (=10! 
cm/s). For air at ordinary temperature and 
pressure: 


Coefficient of viscosity =18.410-* kg/m-s 
Coefficient of kinematic viscosity 
= 1.3210 m?/s, 


AND J. H. COOK 


CONCLUSIONS 


Electroacoustics embodies two branches of 
applied science, namely, (1) mechanics, (2) elec- 
trics including electromagnetics, electrostatics, 
piezoelectrics and thermoelectrics. 

If mechanic quantities are expressed in c.g.s. 
units it is also necessary to measure electro- 
magnetic quantities in c.g.s. electromagnetic 
units if homogeneity is to be maintained. It is 
however, almost impossible to avoid using the 
practical series of electric units (volt, ohm, am- 
pere, etc.) because these are invariably used in 
electric measurements as well as in commercial 
specifications and descriptions. This fact leads to 
the use of a dual system of units, mechanic 
measures being made in the c.g.s. system, elec- 
tromagnetic magnitudes being expressed in the 
practical series. These two systems are not homo- 
geneous, a given quantity having units of dif- 
ferent sizes in the two systems, (e.g. the unit of 
energy in the c.g.s. system is the erg and is equal 
to 10-7 joule, the m.k.s. unit of energy). Thus 
the use of dual systems makes necessary un- 
wieldy transformation coefficients that obscure 
the physical facts being investigated and com- 
plicate the numerical calculations 

When, however, the m.k.s. unit system is used 
for both mechanic and electric quantities all 
transformation coefficients disappear. 

A further advantage of m.k.s. mechanic units 
is their size which is suitable for engineering 
work.5 

At least one textbook® has thus far been pub- 
lished using the m.k.s. system of units exclu- 
sively, and several others are in preparation. 
6A. E. Kennelly, “The M.K.S. System of Giorgi as 
Adopted by the International Electrotechnical Commission 
(I.E.C.) in June 1935,”’ J. Eng. Ed. Dec. (1936). 

6 P. Vigoureux and C. E. Webb, Principles of Electric and 


Magnetic Measurements (Blackie & Son, Ltd., Glasgow; 
Prentice-Hall, Inc., New York, 1936). 
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SYSTEM OF UNITS 


APPENDIX I 
TABLE I. A comparative table of mechanic, acoustic, and electric units of use in electroacoustics. 








c.g.s. UNITS m.k.s. UNITS 


QUANTITY SYMBOL m.k.s. UNIT c.g.s. UNIT IN IN 
1 m.k.s. UNIT | 1 c.g.s. UNIT 
Mass M, m | kilogram, kg gram, gm 103 10-3 
Length L,/2 | meter, m centimeter, cm 10? 10 
Time T,# | second second 1 1 
Density kilogram per cu. meter gram per cu. centimeter, 10-3 103 
. kg/m! gm/cm! 
Velocity V,v | meters per second, m/s centimeters per second, 10? 107? 
cm/s 
Force ‘ ; F,f | newton, kg-m/s? dyne, gm-cm/s? 108 10-5 
Pressure (intensity) p kg/s-m barye, gm/s-cm 10! 107 
Sound Energy Flow - watt erg/s 10? 10-7 
Sound Intensity I watt /m? erg/s-cm? 103 10-3 
Sound Energy Density E joule/m$ erg/cm? 10! 107 
— Acoustic Im- | Z, m.k.s. mechanic ohm, kg/s | c.g.s. mechanic ohm, gm/s 10 10-3 
ance 
“Orifice” Acoustic Imped- | Za kg/s-m‘4 acoustic ohm, gm/s-cm‘ 10-5 105 
ance 
“Specific” Acoustic Imped- | Z; kg/s-m? gm/s-cm? 10-4 10! 
ance 
Viscosity kg/s-m poise, gm/s-cm 10! 1071 
Kinematic Viscosity m?/s cm?/s 104 10-* 
Work W joule erg 10? 10-7 
Power P,P | watt erg/s 10? 10-7 
Elastic Coefficient S newton/m dyne/cm 108 10-3 
Mechanic Impedance Z mechanic ohm, kg/s mechanic abohm, gm/s 108 10-3 
Force Factor a newton/ampere dyne/abampere 10° 10-6 
Electric Current I,i | ampere abampere 10-1 10! 
Electric Resistance R ohm abohm 10° 10-9 
Electric Inductance L henry abhenry 10° 10-° 
Magnetomotive Force F ampere turn 4rX107 =x 10+! 
Magnetic Intensity dC ampere-turn per meter oersted 4xX10-? Exe 
Magnetic Flux Density B® weber/square meter gauss 104 10-4 
Magnetic Reluctance R ampere-turn/weber 4xX 107° x 10° 
Magnetic Permanence e weber /ampere-turn r x 10° 4xX 107° 
Tv 








* J is sometimes used for this quantity. 


APPENDIX II 


An electroacoustic problem illustrating the use of the 
m.k.s. system of units 


The diaphragm of an electrodynamic loudspeaker is 
assumed to act as a rigid piston. We shall consider the 
behavior of this device operating first in a vacuum and 
secondly in air. We shall further subdivide each of these 
cases into operation with an applied mechanic vibromotive 
force and operation under the influence of an alternating 
electromotive fgrce applied to the voice coil terminals 
inducing a mechanic vibromotive force. Mechanic friction 
in the diaphragm and its supports is neglected in the 
analysis. The following constants of the loudspeaker are 
assumed: 


Radius of diaphragm 0.1m 
Area of diaphragm a=0.012 m? 
Mass of diaphragm (including at- 

tached voice coil) M =5 X10-3 kg 
Elastance of diaphragm support S =5 X10? newton/m 
Force factor of electromagnetic 

drive ( =2 newton/amp. 
Electric resistance of voice coil Re =20 ohms 
Inductance of voice coil Le = 107 henrys 


Case Ig Loudspeaker in a vacuum . 
Mechanic vibromotive force applied 
F =2 newtons r.m.s. 


The mechanic impedance of the system is: 
Za=0+j(Mw—S/w) mech. ohms (kg/s). 

With an r.m.s. force, F newtons, applied the r.m.s. 

velocity of the diaphragm is: 
V=F/Z m/s. 

This system is purely reactive, no active power being 
absorbed from the driving agency. (See Table II.) 

We must first determine the impedance to electric 
current flow offered by the moving coil. This impedance 
includes the electric impedance of the coil, when stationary, 
and the motional electric impedance of the device. The 
motional electric impedance is given by: 

Zu=Q?/Z ohms. 


The electric impedance of the voice coil, when stationary, 
is: 
Z.=R.+jlw ohms. 


The total electric impedance of the moving voice coil is: 
Ze=Zu+Z. ohms. 
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TABLE V. 











Case Ia. Loudspeaker in a vacuum. 
Mechanic vibromotive force applied; F =2 new- 


Case IIb. Loudspeaker in air. 
Alternating e.m.f. applied to voice coil E=10 
volts r.m.s. 





tons r.m.s. 
w Z V 
rad/s mech. ohms =kg/s m/s 
100 —j4.5 jA45 
316 0 Fe) 
1000 j4.5 —j.445 
10000 j50 —7.040 


Case Ib. Speaker in a vacuum. 
Alternating electromotive force applied to voice 
coil E= 10 volts r.m.s. 








TABLE III. 








Case Ib. Loudspeaker in a vacuum. 
Alternating e.m.f. applied to voice coil E=10 
volts r.m.s. 








w Ze ZM ZE 
rad/s ohms ohms ohms 
100 20+ 1.0 j.89 20+j1.89 
316 20+ 73.16 0 rv) 
1000 20+710.0 —j.89 20+ 9.11 
10000 20+7100.0 —7.08 20+7100 
I F Vv 
- amp. newlons m/s 
0.496 —j.047 0.992 — 7.094 0.021 + .221 
0 0 - 
414—.089 .828—j.178 | —.0396—j.184 
.019—7.096 .038 —7.192 — .0038 — 7.00076 


Case IIa. Loudspeaker in air. 
Mechanic vibromotive force applied F=2 new- 
tons r.m.s. 





TABLE IV. 





Case IIa. Loudspeaker in air. 
Mechanic vibromotive force applied F =2 new- 
tons r.m.s. 


w Z V P 
rad/s mech. ohms =kg/s m/s watts 

100 13.0—74.5 0.137 + 7.048 0.272 
316 13.0 .154 .308 
1000 13.0+j4.5 137 —j.048 272 
10000 13.0+-j50 .0097 —j.0375 .0195 


Case IIb. Loudspeaker in air. 
Alternating electromotive force applied to voice 
coil E=10 volts r.m.s. 





Hence the electric current that flows in the voice coil 
may be calculated: 


I=E/Zz amp. 





ZM Ze I 
ohms ohms amp 


@ . 
rad/s mech a 




















100 | 13.0—j4.5 |0.275 + 7.095 | 20.275 +71.095 |0.491 —j.0267 
316 | 13.0 308 20.308 +73.16 480 —j.0748 
1000 | 13.0+ 74.5 | .275 —j.095| 20.275 + 79.905 398 —j.195 
10000 | 13.0+750 | .0195 —j.075 | 20.0195 + 799.925 | .0193 —7.096 
F V P n 
newtons m/s watts Q% 
0.982 —j.0534 0.0688 +-7.0197 0.0667 1.36 
960 —j.149 .0738 —7.0115 .0726 1.51 
796 —7.390 .0454 —j.0457 0539 1.36 
-0386 —j.192 — .00341 —j.00166 .000187 097 








This current flowing in the voice coil gives rise to a me- 
chanic force acting on the diaphragm: 


F=QI newtons. 


The diaphragm, vibrating in a vacuum, dissipates no 
energy by its motion. The power supplied to the circuit is 
absorbed in the electric resistance of the voice coil. (See 
Table III.) 

The computations for this case are carried through in 
the same manner as those delineated in case Ia except 
that the mechanic impedance of the vibrating diaphragm 
contains a dissipative term accounting for the tendency of 
the system to supply acoustic energy to the surrounding 
medium, air. The magnitude of this mechanic or acoustic 
resistance is calculated from: 


R=pca mech. ohms (kg/s), 
p=density of air=1.205 kg/m' at 1 atmos., 20°C, 
c=velocity of propagation of sound in air, 
= 344.0 m/s at 1 atmos., 20°C, 
a=area of diaphragm =0.0314 m?, 
whence R=13.0 mech. ohms (kg/s). 


where 


The power dissipated in the mechanic system in this 
case may be determined from: 


P=V?R_ watts. 


This power is supplied to the air as acoustic power. 
(Table IV.) 

The computations for this case parallel those of case Ib 
with the exception noted above that the mechanic im- 
pedance of the system is altered by the air load on the 
diaphragm. 

The efficiency of transfer of electric power to acoustic 
power is the ratio of acoustic power (=mechanic power 
dissipation for this case of frictionless supports) to electric 
power input to the voice coil and is given in the column 
headed 7 in Table V. 


220 OL LR A hl AE MS 





AP 


loc 


ze 





no 
ee 
‘pt 
ym 


ng 
tic 


his 


stic 
wer 
tric 
mn 





PN GO Ss es lt 





APRIL, 






1938 J. A. 


- A. VOLUME 9 


The Position of the Vibrator in the Experiments of Melde and Kundt 


B. J. MILLER Anp L. O. OLSEN 
State University of Iowa, Iowa City, Iowa 
(Received November 10, 1937) 


It is shown that the point of application of the obligatory motion in the experiments of Melde 
and Kundt is at a displacement ‘‘node” in the driven medium. The word ‘“‘node” mzans a point 
of minimum amplitude in the standing wave in the tube or string. This result agrees with that 
of Rayleigh but is obtained in a simpler manner. The discussion is presented briefly to make 
the results more available as current textbooks indicate a nearly universal misunderstanding 


of these two well-known experiments. 





N examination of a number of widely used 
college physics texts shows a quite general 
misunderstanding of the position of the vibrator 
in the Kundt’s tube experiment. Many writers 
locate the end of the vibrating rod at a displace- 
ment loop of the standing wave in the tube. 
Others locate it “‘near a node.”’ Still others avoid 
any statement but indicate graphically an agree- 
ment with one of these locations or a position 
midway between a node and a loop. 

The purpose of this discussion is to point out 
the correct view, namely, that at resonance the 
distance from the piston to the nearest displace- 
ment node is zero except for terms involving the 
second and higher powers of the damping con- 
stant. Here the word ‘‘node’’ implies a point of 
minimum amplitude in the standing wave in the 
tube. These minimum amplitudes differ from 
zero except at the fixed end of the air column. 

Discussions of Melde’s experiment in these 
texts do not usually include the point here con- 
sidered; the illustrations presented are, for the 
most part, correct in this respect. It is probably 
not generally realized however that the end of 
the fork is located so precisely at a displacement 
node when the string is tuned for resonance. 
Thus it is worth while to point out that the 
treatment here presented applies equally well 
to the transverse form of Melde’s experiment. 
(Fork vibrating at right angles to the length of 
the string.) 

We make the following assumptions and con- 
ventions: 

(1) The fixed end (of the air column or string) 
is located at x=0. 

(2) The vibrator is at x=/ and executes a 
vibrator described by 


341 


u=a sin wt, 


where u is the displacement. The reaction of the 
driven medium is thus neglected. 

(3) A damping force proportional to the par- 
ticle velocity is assumed. 

Under these assumptions the wave equation 
takes the form: 


Ou 2Kou 1 du 
—-—_=__, (1) 
ox? 0c «Ct CO «(OF 


where c is the (7'/m)! in the case of the string 
and (yp/p)' for the air column. In order to 
simplify the subsequent algebra, the constant 
of proportionality is written in the form 2K/c. 
This equation is to be solved subject to the 
boundary conditions. 


At x=0, u=0 for all #, (2) 
At x=1, u=a sin wt. (3) 


A solution which satisfies these conditions can 
be derived from the solution 


u=Ae—™ sin k(x— Bt) + Be-* cos k(x — Bb) 
+Ce* sin k(x+t)+De™ cos k(x+ Bt). (4) 
c 
Here a=——————_-_ and s=-——_——_- 
(1+K?/k?)! (1+K?/k?)! 


k=2n/d, where \=27(8/w); and the positive 
square root is used for @ and 8B. 
Application of the boundary conditions gives: 


A=C, B=-D, 
A tan kl=B tanh al, 


a/2=A cos kl sinh al+B sin kl cosh al. 
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This gives u in the following form: 
u=a/(cos* kl sinh? al+sin? kl cosh? al) 
cos wt(sin kx cosh ax cos kl sinh al 
—cos kx sinh ax sin kl cosh al) (5) 
+sin wi(cos kx sinh ax cos kl sinh al 


+sin kx cosh ax sin ki cosh al). 


This can be further reduced if we set the coef- 
ficient of cos wt equal to —Rsin y and the 
coefficient of sin wt to R cos y. This gives 


aR sin (wt—vy) 


05 IR eter emeerenmeeeenemenenmeemeoe, (6) 
cos? kl sinh? al+sin® kl cosh? al 
R reduces to 
(sinh? ax cos? kx-+cos? ax sin? kx)! (7) 


~ (sinh? al cos? kl--cosh? al sin? ki)! 
and y is given by 


tanh ax tan kl—tanh al tan kx 
tan y =—_________———_ (8) 
tanh ax tanh al+tan kx tan Rl 


Finally u becomes! 


a(sinh* axcos*kx+cosh?axsin*kx)* 
.=—$—$—$$—$— sin (wi—¥). 
(sinh? al cos? kl+cosh? al sin? kl)! (9) 
In this form the physical interpretation is 
evident. Each particle of the air column or 
string executes a simple harmonic vibration with 
the frequency of the vibrator and with an 
amplitude 


a(sinh? ax cos? kx-+cosh? ax sin? kx)! aF(x) 





(sinh? al cos? k/+-cosh? al sin* kl)! F(l) 


and with a phase given by y. 

The solution differs from that for an undamped 
medium in that the dependence of the displace- 
ment on x is not sinusoidal and the phase of the 
vibration is a function of x. 

For any frequency of the vibrator there is a 
vibration in the string or air column having 
minima and maxima. The points of minimum 


1The result agrees with that of Rayleigh, Theory of 
Sound, second edition, Vol. I, page 198. 
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amplitude are the minima of F(x) given by 
a sinh 2ax=—k sin 2kx. (11) 


If a is small 
—2a’x=k sin 2kx. 


If we set x=m7/k+e, where 2ke <1 
€:= (—n7a’)/k'. (12) 


Thus the nodes are at x=n7/k if a?/k? is very 
small compared to unity, and are uniformly 
spaced to an even closer approximation.’ 

The location of the vibrator (x=/) relative to 
these minima is as yet not fixed as we have 
not used the condition of resonance. We define 
the resonant frequency as that which leads to a 
maxima displacement amplitude. This can be 
determined by minimizing F(l) with respect to 
k. This gives 


sin 2k] = (2a*/k*) sinh 2al (13) 


or for small a 
sin 2kl=2a‘l/k® 


setting kl=nr+kee (14) 
2keo= (nx/k+ €2) (2a*/k*) 
€.=nra'/k', (15) 


Thus at resonance 
k=(nr/l)(1+a*/k'), (16) 


which expression is identical for all practical 
purposes, with that calculated by more ele- 
mentary means. 

The distance of the vibrator from the nearest 
node can now be exhibited. From (14) 





nx nat a 
l=— » or, if—<1, 
k k® k? 
nT 
ik 


From (12) the nearest node is at x=mz/k 
—nro?/k® 
nro na’r® 


(17) 





kB Ba 


if we neglect higher powers of a. 





l—-x=d= 


2 To a closer approximation 


2 2 24-2 
a=" [145 (= -1)]- 


Thus we have a linear dependence of « on if a*/k*C1. 








AO er SS A MME - 





vi 
nt 


oo 
pc 





5) 


16) 


cal 
le- 


est 


17) 


ee ee ee ms 


VIBRATOR 


Thus the node at resonance is in front of the 
vibrator a distance proportional to a? and to the 
number of nodes. 

An estimate of the value of a can be obtained 
from the following considerations. If we insert 
“k” from (16) in Eq. (9), and neglect higher 
powers of a/k, we may write 


aF(x) 
“= 





sin (wt—y). 18 
sinh al (18) 


Also, to an approximation close enough for this 
computation (first power of a) 


sin kx=1 


at a loop. At such a point the amplitude is given 
by a/al. The amplitude of the vibrator is a and 
the ratio of the amplitude of the driven medium 
to that of the driven is thus 1/al. 

Suppose we use as an example some actual 
data on Melde’s experiment. The string used 
was a high grade of fishing tackle, but inhomo- 
geneities were apparently not entirely absent. 
However it gave more consistent results than 
various other strings and wires tried. The am- 
plitude of the string at resonance was con- 
servatively estimated to be ten times that of the 
fork. For a length of 200 cm this gives a=5 
X10-*. In the most adverse case \=100 cm 
n=4. From Eq. (17) then ‘‘d’’=0.013 cm. 

Thus the distance from the fork to the nearest 
node is theoretically 0.013 cm. This is entirely 
negligible in view of the fact that the various 
internodal distances differed from each other by 
one or two cm. 

In the case of Kundt’s experiment the am- 
plitude ratio is much greater and the customary 
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wave-length much shorter. leading to even 
smaller values of d. 

In Melde’s experiment it is usually found that 
the resonance point is not sharply determined by 
the amplitude which is essentially constant over 
a fairly wide range of load. Experimentally it is 
found that the resonance load can be much more 
sharply determined by requiring that the final 
node coincide with the end of the fork. 

The phase factor y is of some interest. In the 
vicinity of a node, from Eqs. (8), (12) and (15), 
tan y= ©. That is, the vibration at the nodes is 
x/2 radians out of phase with the vibration of 
points remote from the node. This was first 
noticed and discussed by C. V. Raman.’ 


SUMMARY AND CONCLUSION 


(1) An expression has been developed for a 
standing wave in a dissipative medium with an 
obligatory simple harmonic motion imposed at a 
point. 

(2) It has been shown that the points of 
minimum amplitude in the driven medium are 
equally spaced if powers of the damping constant 
higher than the second are neglected. This spac- 
ing is equal to that for the undamped case if the 
square of the damping coefficient is negligible. 

(3) At resonance, the mean position of the 
vibrator coincides with that of a minimum 
if a?/k®<1. 

(4) A method of estimating damping con- 
stants and a method of tuning to resonance have 
been suggested 

The authors wish to acknowledge their in- 
debtedness to Professor G. W. Stewart who called 
their attention to the general misunderstanding 
and the correct view. 
~ #C. V. Raman, Phys. Rev. 32, 309 (1911). 
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The Sound of Bells 


Jottings from my Experiences in the Domain of the Sound of Bells 


Jan Arts 
“Huize Ruwenberg” St. Michielsgestel, the Netherlands 


(Received September 25, 1937) 


THE STRIKE NOTE OF BELLS 


HE strike note of bells has been discussed by 

A. T. Jones,! and it may be assumed there- 

fore that the reader is acquainted with the strike 
note problem. 

It is some ten years ago now that I wrote to 
the late P. Griesbacher as follows: “‘I wonder 
why you don’t mention the seventh partial in 
Glockenmusik. Yet experience of the last few 
years taught me that this apparently unsig- 
nificant matter is of great importance for the 
exact determination of the strike note.’” 

From the moment that I noticed this phe- 
nomenon, originally strange but explainable, I 
have always, whenever an opportunity offered, 
noted the seventh partial together with the fifth 
partial with a view to be able to make out 
whether the mysterious strike note might be 
explainable as difference tone. But soon I verified 
that this could not be true, the tone-relations of 
the 5th and 7th partials being usually of such a 
nature as to exclude the possibility of the strike 
note being a difference tone. 

In order not to take up too much space I 
submit three tables, which will sufficiently prove 
and explain my allegation. I confine myself to 
these three though I could add a dozen similar 
cases. 

Table I shows the octaves below the fifth 
partials and the frequencies of the fifth and 
seventh partial of two bells in the tower at 
Maestricht, Limburg, the Netherlands; when 
these bells sound together, the two strike notes 
are heard in the interval of a minor third. 
Careful experimenting showed that the 5th and 
7th of bell I are represented as d’’ 8/16-a” 10/16, 
respectively, the 5th and 7th of bell II as 
f” 7/16-c’" 5/16. . 

This table clearly shows that if the strike 
note should be originated by the fifth-seventh, 
the strike note as difference tone is excluded; 
should this be the case, not a minor third would 


1A. T. Jones, J. Acous. Soc. Am. 8, 199-202 (1937). 
2 P. Griesbacher, Glockenmusik (Nachtrag, 1929), p. 43. 


be heard when the two bells rang together, 
but the major second d’ 13/16-f'0/0. The 
octaves d’ 8/16-f’ 7/16 below the fifth partials 
on the other hand are quite concordant with the 
minor third interval heard. 

Table II shows the octaves below the 5th 
partials and the tone relations of the 5th and 7th 
of two bells in the tower of Vlokhoven, Noord- 
Brabant. When the two bells are rung together, 
a good major second interval is heard. Were the 
difference tones ep’ 7/16-e’ 13/16 to decide of 
the strike notes, not a major but a minor second 
would be heard. So here too the octaves below 
the fifth partials are in favor of the strike note 
hypothesis. 

The tone relations of the 5th and 7th in 
Table III are of two bells, of Breda, Noord- 
Brabant. When rung together, here too the 
major second is heard. Again, were the differ- 
ence tones bb’ 7/16—-db’’ 1/16 to decide of the 
strike notes, the minor third and not the major 
second would be heard. 

Conclusion: ‘“The results of this test appear to 
favor the octave hypothesis rather than the 
difference tone hypothesis.” 

I for one am of opinion that the fifth partial, 
both in the case of small bells and of big ones, 
originates the strike note. In defining the ob- 
jectively nonexisting strike note, my only guide 
is the fifth; a guide fully to be relied on even in 
the most critical cases. 


























TABLE I. 
BELL I BELL II 
CycLes/ CyYcLes/ 
SEC. Note* SEC. Norte* 
Seventh 901.6= | a” 10/16||1053.6= | c’” 5/16 
Fifth 597.6= |d” 8/16]| 708.2=|f" 7/16 
Difference tone 304.0= |d’ 13/16)| 345.4=|f' 0/0 
Octave below the 
fifth partial 298.8= |d’ 8/16}| 354.1=/f’ 7/16 





* The fraction beside the name of the note is part of a half-step of an 
equally tempered tuning fork. The fraction on the left of the note- 


name points to a minus, the fraction on the right to a plus. 
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TABLE II. TABLE III. 
BELL I | BELL II Bett I ‘ Bett II 
CYCLEs, CyYcLes/ CyYcLes/ CycLes/ 
SEC. NoTE SEC. NoTE SEC. NOTE SEC. NoTE 
Seventh 935.2= | bb’’4/16||1042.0= |c’” 2/16 Seventh 1421.6= | f’" &/16)|1618.8= |g’ 12/16 
Fifth 619.6= | eb” 2/16|| 700.6= |f” 4/16 Fifth 948.8 = | bb’ 8/16}/1068.8= |c’”’ 9/16 
Difference tone 315.6= |eb’ 7/16|| 341.4= |e’ 13/16 Difference tone 





Octave below the 
fifth partial 











309.8= |eb’ 2/16)| 350.3= jf’ 4/16 








472.8= | bb’ 7/16)| 550.0= |db” 1/16 





Octave below the 
fifth partial 





474.4= | bb’ 8/16)| 534.4= |c” 9/16 














The late P. Griesbacher once asked me? 
“Just tell me: How do you test bells higher 
than a’’? As our tuning forks reach only to a’”, 
the method with the fifth leaves you in the 
lurch already at b’’. And is it your opinion too, 
that at that height the strike note somewhat loses its 
influence on the sound and that often the second 
and oftener still the first partial is predominant?” 

My answer: Last year I was sent two small 
bells for examination, bells of exactly the same 
diameter (pl.m. 40 cm) but formed after different 
profiles. The two bells were before me with 
mouths turned upwards. When struck respec- 
tively with a mallet there was no doubt about 
the strike note of the A bell being decidedly a 
quarter of a step below that of the B bell. 
After careful experimenting I obtained the tone 
relations as shown in Table IV. As the 5ths of 
the two bells were both of the same frequency 
the strike note concordance must be all right, 
I thought. At my request I had the two bells 
mounted for ringing and lo: when the two bells 
rang together it was as if only one bell was 
ringing, so perfect was the concordance between 
the two strike notes. So originally the distinct 
coming forth of the two first partials, differing 
about a quarter of a step, was misleading.‘ 


Is A BELL TONE, OBTAINED WITH A FORK, 
TUNED AN OCTAVE LOWER, AN ABSOLUTELY 
RELIABLE PICTURE OF THE FRE- 
QUENCE RELATION 1 : 2? 


For many years I have been using a set of 
tuning forks (Edelmann) with a total range of 


3 Letter: 20, X (1930). 

4My Dutch readers are informed that the firm of 
Petit & Fritsen, bell-founders at Aarle-Rixtel, N.B. found 
the case so interesting and convincing that they keep the 











TABLE IV. 
A BELL B BELL 
Fifth ec’ 12/16 c’’”’ 11/16 
Fourth g’”’ 12/16 g’”” 12/16 
Third eb’”’ 10/16 eb’” 12/16 
Second ec sre 2/16 c’”’ 
First 2/16 c”’ ee \2sne 








F# to a’” included, to investigate bell tones, of 
both church bells and carillon bells. As I could 
not afford to acquire the whole set at once on 
account of the great expense, I bought only one 
fork, the fork with a range of d#’ to a’. As the 
range of the fork was rather small, the fork was 
of little use to me; the range of the tones of 
church bells and carillon bells being very ex- 
tensive. Therefore I asked myself the question: 
Is it possible to bring forth the bell tones lying 
outside my tuning fork, as a harmonic overtone. 
And behold, by my first experiment I found that 
with a fork tuned an octave or duodecimo lower, 
the bell tone, which was an octave or duodecimo 
higher, reacted clearly and distinctly upon it. 
But a second question arose: Is the result arrived 
at in this way to be trusted; is the aliquot tone 
produced in this way, an absolutely reliable 
picture of the frequence relation: 1 : 2, 1 : 3? 
Gradually I came into the possession of the 
whole set of nine forks. I thought now to be 
able, to prove the reliability of my method. 
I experimented on the bell in the belfry of our 
college, with a diameter of about 38 cm. This bell 
has the 4 following partial tones: first: c’’ 8/16; 
second: b’”’ 15/16; third: eb’’’ 7/16; fourth: 
gb’”” 15/16. 
two bells mounted in their show room in order to give 


those interested at any time an opportunity to verify 
the fact. 


a —— 











TABLE V. The tones in column 1, were acquired with the 
normally pitched tuning fork. The fork tones in 
columns 2, 3, 4 and 5, bring forth as first aliquot 
tones, the tones of column 1. 








Fifth ec’ 6/16| c’” 6/16 | f” 6/16] c” 5/16 | ab’ 8/16 
Fourth || gb’” 15/16] ¢” 0 b’ 15/16 | gb’ 15/16 | eb’ 0/0 

Third ||eb’” 7/16| e ab’ 8/16 | eb’ 7/16|b 9/16 
Second ||}b’” 15/16} c” 0/0 e’ 15/16 | ce’ 0/0 
First ec” = 8/16| & + 8/16 | f 10/16} c 10/16 | Ab 10/16 


When setting the fork an octave lower to the 
first, so to c’ 8/16, the first partial very dis- 
tinctly was brought out. The second b” 15/16 
reacted very weakly on b’ 15/16 but distinctly 
came out when the fork was tuned on c”’ 0/0. 
The third eb’” 7/16 clearly answered the setting 
eb’ 7/16. The fourth gb’ 15/16 on the other 
hand, only clearly came forth on a g’’ 0/0 
setting. Now I tried whether I could elicit the 
‘same partials (plus the fifth) also as 3rd, 4th and 
5th aliquot tones. In this I succeeded perfectly, 
but, as Table V shows the results are not all in 
the exact frequence relations: 1:2, 1:3, 
£24, 1:5. 

In all my subsequent experiments I observed 
the peculiar phenomenon, that sometimes a fork 
tuned a pure octave lower, produced the pure 
over-octave and sometimes not. How to account 
for it, that the results are so different? Is the 
cause to be found in the physical nature of the 
bell or of the tuning fork? (Table VI.) 


AN UNKNOWN (?) BELL TONE 


I am in the possession of a small bell (weight 
about 12 kg). One day as I was experimenting on 
this bell I suddenly noticed a tone I had not 
heard so far. Since that time I regularly meet 
with this unknown (?) tone in small bells. This 
tone has the peculiarity (a) it is very modest in 
its appearance, (b) it does not react upon the 
tuning fork. It is possible, however, to make it 
free of beats to respond to the tuning fork and 
consequently its pitch exactly definable. Its 
singing-out time is of the same length or nearly 
so as that of the 3rd partial and therefore I 
presume that it is closely related with the third 
partial. 

5 In order to arrive at a wholly reliable result especially 
in the case of the lower bell tones, such as the first and the 


second, it is desirable to put a piece of thin cloth on the 
soundbow, to prevent the fork from jumping up. 
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TABLE VI. The partial tones in column I were acquired with 
the normally tuned fork; the fork tones in column 2 
bring forth (with a fork tuned a pure (?) 
octave lower) the tones of column 1. 























a - 2 
Third fH" 865/16 | f# 5/16 
Bell I Second d#” 5/16 | d# 3/16 
First d#’ § =67/16 | d# 8/16 
Fourth gH” 11/16 | g#” 11/16 
Third e” 8/16 | € 5/16 
Bell II Second c#’""" 7/16 | c#” 4/16 
First c#” 7/16 | c# 7/16 
Fourth g” 8/16 | ¢” 8/16 
Third eb’” 10/16 | eb” 8/16 
Bell III Second 4/16 c/” 2/16 c”’ 
First Se” 6/16 | c’ 8/16 
Fourth g’” 11/16 | g” 10/16 
Third eb’” 9/16 | eb” 7/16 
Bell IV Second ce” #8735 | <” 4/16 
First 4/16 c” 1/16 c’ 
Fourth r" i796 1 1/86f" 
Third db’” 7/16 | db” 4/16 
Bell V Second bb” 1/16 | bb’ 2/16 
First bb’ 1/16 | bb 2/16 
Fifth +" 0/0 gz” 0/0 
Fourth d#’” 0/0 6/16 d#” 
Bell VI Third b” 0/0 b’ 0/0 
Second 1/16 g#” 3/16 gz’ 
First 1/16 g#’ 2/16 gz 


I should be inclined to call it the ‘‘acoustic 
tone’”’ on account of its objective nonexistence. 
The first and second partial of my study- 
object are affected with beats. The third too has 
this defect. Now it is a well-known fact, that, if a 
partial sings out with beats, there are two sector 
groups, which each of them send out separate 
tones, only slightly differing in pitch. 

As it happens, with my bell the sectors of the 
ist and 2nd are quite congruent.® If I make the 
bell vibrate by means of the clapper in one of 











TABLE VII. 
| 

SINGING- SINGING- 

OUT PARTIALS OUT 
SECONDS A SECTORS ACcousTIC B SECTORS SECONDS 

0.5 f’’""" 15/16 Fourth 18/16 0.5 

3. d’”’ 4/16 Third er" 3/% a 

Si b’’” 1/16 Second ye” S/% 5. 

2 g’”” 2/16 Acoustic g#’” 0/0 3. 

15. 5/16 b’”’ First 8/16 b’”’ 15. 


6 Incidentally I wish to draw the reader’s attention to 
the possibly unknown fact that the position of the nodal 
meridians and specially that of the nodal circles is equally 
well definable by means of the adjustable tuning forks as 
with the resonators or any other instruments adapted to 
the purpose. 





act, rim FO lat 





SOUND OF BELLS 





; TABLE VIII. 
Bell 1 Bell 2 Bell 3 Bell 4 Bell 5 Bell 6 Bell 7 
! 
P.A. N. SS. N. S.S. N. S.S. N. S.S. N. S.S. N. SS. N. SS. 
Third ce” 0/0 9 ||eb’”” 10/16] 7 || eb’’” 9/16] 7 |leb’” 1/16| 7 ||f"" 0/0 | 6 |] db’ 7/16| 11 || 0/0 b” 11 
Second 1/16 a” 12 || 4/16 c’”" 11 |] ce’ 4/16] 9 ||5/16¢'" 10 ||d’’” 7/16] 10 || bb” 1/16 | 13 |] 1/16g#” | 14 
- Acoust. ” 0/0 8 jig” 1/16 | 7 || ab” 4/16] 6 ile” 10/16] 7 ||ab” 0/0 | 6 || f#” 2/16] 11 || e” 1/16 | 10 
First a’ 1/16 | 27 ||c” 6/16 | 25 || 4/16c” 26 | 3/16’ 30 || 12/164’ | 22 || bb’ 1/16] 47 || 1/16g#’ | 36 





The initials P.A., N. and S.S. respectively indicate: partial/acoustic, note and singing-out seconds. 


the A sectors, the elusive tone comes forth as 

g’’’ 2/16; if the bell is struck in the middle of 
° one of the B sectors another tone presents itself, 
namely g#’’"’ 0/0. In the beginning you have to 
strain your ears to discern it from the other bell 
tones; once its acquaintance thoroughly made, 
it is rather easily discernable. 

A synoptic table of the tone relations and 
singing-out-time is given in Table VII. For 
clearness sake of the above-said, I shall call the 
two sector groups A and B. Now I began experi- 
menting on a 50 kilogram bell and discovered 
between the ist partial c’ and the second 
partial c’” a tone of analogous character ¢”’, 
which took up 4 seconds in singing-out. Since 
that time I have often been observing this 
“super” but have not succeeded so far in de- 
fining its origin. I met it not only in bells of 
; Netherlands make, but also in many bells, 

found abroad. 
I give by way of elucidation (Table VIII) the 
| ist, 2nd and 3rd partials, besides the acoustic 


we 


tones with the singing-out times proper to them, 
j of seven bells. May I have stimulated hereby to 
; the study of this problem, of which, up till now, 


I have not been able to give a satisfactory ex- 
planation. 
. TEMPERATURE INFLUENCE ON BELLS 
The late Griesbacher came by experiments’ to 


7 If one wants to make this experiment by means of the 
tuning forks, it is necessary to handle it carefully as it is 
rather sensitive to temperature changes. 








the conclusion that the bell is not subject to 
temperature changes.® 
Some years ago—it was an exceptionally hot 

August—I took in the evening a bell into a 
room and tuned the forks carefully i.e., free of 
beats to the ist, 2nd and 3rd partial. The next 
day in the forenoon, the thermometer in the 
room where the tuning forks and the bell had 
spent the night, pointed to 18° Celsius. I con- 
trolled the fork tones together with the bell 
tones and the agreement was complete as the 
evening before. The forks were left in the room 
(18°C) and bell and thermometer were exposed 
to the sun for an hour. After an hour the ther- 
mometer pointed to 31°C. Now I brought the 
forks into contact with the bell one after another 
and counted 

for the ist partial 14 beats per 10 seconds 

Oo eee eee ee be 
Ra | | eee ” 
Some time afterwards I came to the same result 
on another bell: 

for the 1st partial 10 beats per 10 seconds 

a Teo” —_- " 
a | hme. Se " 
These results clearly prove the influence of the 
temperature on the bell body and consequently 
on the tones produced. Incidentally I wish to 
draw my readers attention to the ratio of beats: 
14, 29, 33, and 10, 20, 24 which is in agreement 
with frequence proportion: 1:2 (first-second) 
1 : 2, 4 (first—-minor third). 


8 Reference 2, pages 46-47. 
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The Velocity of Sound in Single Crystals of Bismuth 
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Brown University, Providence, Rhode Island 


(Received December 10, 1937) 


The velocity of compressional waves in single crystals of bismuth of two orientations (P; and 
P;) of the principal axis has been measured. The specimen to be tested (around 10 cm in 
length) is inserted between two longer steel rods, one of which is driven by a loudspeaker coil 
while the other has a free end. The velocity is computed from the observed distance between 
two nodes when all or part of the crystal is included between them. For the P; and P; orienta- 
tions the velocities are found to be 2028 meters/sec. and 1541 meters/sec., respectively, and 
the corresponding Young’s moduli come out to be 4.04310" dynes/cm? and 2.24510" 


dynes/cm?. 





INTRODUCTION 


HE present paper contains the results of the 
measurement of the velocity of compres- 
sional waves in single crystals of bismuth of 
’ various orientations of the principal axis. The 
measurements not only possess acoustical interest 
but yield a valuable method for the calculation 
of Young’s modulus and its dependence on 
crystal orientation. The method is believed to be 
particularly adapted to acoustical wave velocity 
determinations in small specimens of soft and 
brittle substances which are able to stand but 
little strain. 


THEORY OF THE METHOD 


Consider a solid rod of length J, and cross- 
sectional area S» inserted between two long rods 
of cross-sectional areas S; and 53, respectively 
(see Fig. 1). The middle rod is the crystal to be 
tested, while for convenience the other two 
rods are of some standard easily prepared metal, 
e.g., steel or brass. If the rods are the seat 
of compressional waves traveling lengthwise 
through them, the particle displacement £ and 
excess tension T at the boundary between 5S; 
and S2 will be given in terms of these quantities 
at a point 1, distant /, from the boundary 4A, 
by the relations! 





Fac. 1. 


1Cf., e.g., Lindsay and White, J. Acous. Soc. Am. 4, 
155 (1932), particularly Eqs. (8) with necessary change in 
notation. 
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£:o= & cos kili+(T1/piciw) sin kyl, 
Tw2= T; COs Rili— pice sin Rily, 


(1) 


where the quantities £; and 7, are taken at point 
1 and £12. and Tj. are taken at the boundary A 
but zm the first rod. The density of the first rod 
and the velocity of sound in it are given by p; 
and ¢1, respectively. As usual ky=w/cy. 
Denoting quantities at the boundaries A and 
B of the specimen rod by the subscripts 21 and 
23, respectively, we have the similar equations 


£o3= £1 Cos Rolo +(T21 ‘pol aw) sin Role, 
T:3= T2 cos Role — poCowéas sin Rolo, 


(2) 


where ko=w/cz and pz and ce refer to the speci- 
men being tested. Denoting by the subscript 3 
the displacement and excess tension at a point 
distant /; in the rod of cross section S; from the 
boundary B, we finally have the third pair of 
relations 


= E30 Cos Rils+(T32/piciw) sin Rils, 
T3= T 32 cos Ril3— pycywese sin Rls. 


(3) 


The boundary conditions at A and B, re- 
spectively, are 


f1o= £21 ’ 


fos=Es2 5 
4 
SiTi2= S2T 21; ( ) 


SoT 23 _ S37 32. 


With the use of (4) we can express £3 in terms of 
£, and 7,. The result is 


f= cos bss( cos Role cos Ryls 
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9 
~ 


Suppose now that the points 1 and 3 are 
nodes. Then £;=£=0 and hence the coefficient 
of T, in (5) must vanish. The condition that 1 
and 3 are nodes is therefore (taking the special 
case where S,;=5S3) 


tan k,l,+tan k,l; 


tan kel. = —————_—_—__—_ 
1/A—A tan kj]; tan kyl; 


with A = Sopece/S1p1c;. Let us now introduce the 
parameters 7; and 3 where 


tan mkiWl,=A tan kil, (7) 
tan n3kils=A tan Rls. 

This enables us to write (6) finally in the form 
Rolo+ mikilitnskils=nr, (8) 


where is an integer and in particular is equal 
to unity if 1 and 3 are successive nodes. Eq. (8) 
can be used to calculate ke and hence ce from a 
knowledge of the physical quantities already 
specified and measurement of /;, /2, 3. 


EXPERIMENTAL PROCEDURE 


The bismuth single crystals used were grown 
by a method described by one of the authors in a 
previous paper? and hence no discussion of this 
need be given here. Two types of crystal orienta- 
tion were employed, viz., (1) the P; orientation 
in which the principal axis of the crystal is 


2Cf. A. B. Focke and J. R. Hill, Phys. Rev. 50, 179 
(1936). 
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FIELD SUPPLY 
6 VOLTS 


CLOCK 


perpendicular to the length of the rod, and (2) 
the P; orientation with the principal axis parallel 
to the length of the rod. The crystal rods used 
were all of the order of 10 cm in length and 0.04 
cm? in cross section. 

The actual experimental arrangement is indi- 
cated in Fig. 2 which shows the crystal rod AB 
soldered between two sections of steel, S and S’. 
The former section S is fixed by a collar C toa 
specially constructed voice coil attachment, V, 
of a dynamic speaker unit L. S. This was driven 
by a 15-watt amplifier, A, with a practically 
constant output over the audible range. The 
General Radio oscillator F gave a continuous 
variation of frequencies in the desired range 
with fair frequency stability. 

As vibration detectors two piezoelectric (Ro- 
chelle salt) phonograph pick-ups (P in Fig. 2) 
were employed. These have a very nearly uni- 
form response up to 12,000 cycles. Since their 
output is small it must be amplified. This was 
done by feeding the current into the amplifiers 
T, the first stage of which was located in a metal 
shielded can on top of which the pick-ups were 
mounted. In this stage a type 6C6 tube was 
used. The output of the amplifiers was impressed 
on the vertical plates of a cathode-ray oscillo- 
graph, O, through resistance coupling. By means 
of the double pole, double throw switch Y, the 
pick-ups could be connected in turn to the 
oscillograph. To measure the frequency of the 
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vibrations of the crystal, the output of a beat 
frequency oscillator, B, was placed across the 
horizontal plates of the oscillograph. The fre- 
quency of this oscillator was varied between 700 
and 1500 cycles until a simple Lissajous figure 
was observed on the fluorescent screen. This 
made possible the determination of the ratio of 
the frequencies of the crystal and the oscillator. 
At the same time the output of the oscillator was 
fed into a two-stage amplifier, D, which drove a 
1000 cycle synchronous clock, N. If the time 
required by the second hand to make a complete 
revolution is denoted by ¢, the frequency of 
the oscillator is 


v’ =60,000/t. 


If uw is the Lissajous ratio the frequency of 
vibration of the crystal v is then given by 


v=pv'. 


In the experimental procedure the rod S was 
set in vibration by the oscillator F through the 
voice coil V, the pick-ups P were moved along 
the rod and the nodal points noted. With the 
indicated arrangment it was possible to locate 
the nodes to within less than one millimeter. 
It was found that the nodes could be most 
easily located with the rod vibrating in a reso- 
nance mode. 

The velocity of sound in steel enters into the 


TABLE I. 
Crys OrI!- 
TAL |LENGTH|ENTA- |FREQ. n(IN c2IN 
No. | IN cM | TION vy |licm|lecm|ls; cm/EgQ. (8))| M/sEc. 
1 12.3 Pi 6860 | 16.8 | 12.3 | 22.3 2 2140 
6900 | 16.7 | 8.1 0 1 1848 
6950|17.2 | 7.5] 0 1 1775 
2 11.3 Pi | 6860] 16.55} 7.8] 0 1 1753 
7420 | 20.7 | 11.3 | 14.5 2 1900 
3 10.8 Pi {10220} 10.2 | 10.8 | 14.5 2 2130 
11520 | 17.1 | 10.8 | 41.1 3 2070 
8330 9.9 9.4 1 1910 
6850 | 16.7 | 9.2} 0 1 2040 
4 11.1 Pi 7490 | 0 6.9 | 17.0 1 2025 
11320 | 15.3 | 11.1 | 18.6 3 1960 
11840 | 19.2 | 11.1 | 14.0 3 2035 
4860/ 0 9.6 | 27.5 1 2050 
5 11.6 Ps | 6850| 16.7 | 7.3] 0 1 1525 
8520 | 17.65) 11.6 | 27.0 3 1490 
8460 | 56.9 | 11.6 | 17.4 4 1460 
9270} 2.9 | 11.6 | 13.5 2 1560 
6 9.6 Ps | 9270] 44] 9.6 | 18.4 2 1583 
10120} 8.7 | 63] 0 1 1630 
7 11.0 Ps | 6300] 11.2 | 10.5] 0 1 1500 
6780 | 15.7 8.0 0 1 1517 
11700 | 0 6.85) 0 1 1610 
11280 | 0 68 0 1 1532 
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formula (8). This was measured by the simple 
method of determining the distance between two 
successive nodes in the vibrating steel rod and 
using the relation c=c,=),». This simple method 
cannot, of course, be used successfully in general 
in the bismuth crystal because of the short 
length and great fragility of the latter. 


RESULTS AND CONCLUSIONS 


The averaging of twelve measurements rang- 
ing in frequency from 3300 cycles to 8700 
cycles gave for the velocity of the compressional 
wave in the steel rod used the value c,=5141 
meters/sec. All the individual values were 
within 1 percent of the mean. 

The following data were used in the reduction 
of the experimental observations. For the steel 
pi=7.753 gram/cm’, S,=0.0445 cm?; for the 
bismuth crystal pe=9.83 g/cm’, S2=0.04 cm’. 

In the application of the experimental meas- 
urements to Eq. (8) it was found desirable to 
use linear interpolation. To render this prac- 
ticable without undue error it was necessary to 
take the readings in such a way that k,l, and 
k,l; are not close to 0 or 7/2 (or higher multiples 
of 7/2). In general it was found advantageous to 
have tan k,l; and tan k,/; differ in sign. Moreover 
the choice of A <1 (by proper choice of the cross- 
sectional areas) helps to avoid large variations 
in C2 due to small errors in the measurements. 

The probable error of the frequency measure- 
ments was 0.25 percent, and that in the length 
measurement approximately 1 percent. 

The bismuth crystals are extremely fragile and 
twin lamellae were apt to be produced by strain- 
ing the crystal. Any twinning of a crystal in the 
P, orientation changes the orientation in the 
direction of the P3; orientation and vice versa. 

Seven crystals were tested, four in the P; 
orientation and three in the P3. Table I contains 
the results of the measurements in the order in 
which they were actually carried out for each 
crystal. The occurrence of J=0 in the table 
means that in the particular measurement 
concerned no node was used in the rod in ques- 
tion. The fundamental equation (8) still remains 
valid in this special case. The last two measure- 
ments listed are of interest because in them both 
nodes were located in the crystal, and the formula 
(8) thus reduces to a particularly simple form, 
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for /, in these cases is simply a half-wave-length. 
It will be observed that /, is equal to the actual 
length of the crystal only when »>1. 
Examination of the table discloses that the 
velocity c2 for the P; orientation is definitely 
greater than that for the P; orientation. More- 
over, for all crystals save 1 and 2 the results 
for various frequencies and / values are reason- 
ably consistent. In the case of crystal 1, the 
great discrepancy between the second and third 
values and the first (which agrees fairly well with 
the average of 3 and 4) can be accounted for by 
the assumption that the strain on the crystal 
produced twinning which would result in a 
lowering of ce in the direction of the value for 
the P; orientation. This assumption was verified 
by actual etching and microscopic examination 


IN BISMUTH CRYSTALS 


TABLE II. 





ORIENTATION AVERAGE VELOCITY C2 


P, 2028 m/sec. 
P; 1541 m/sec. 


YOUNG'S MODULUS 


4.043 X 10" dynes/cm? 
2.245 X10" dynes/cm? 





of the crystal. A similar situation was found to 
exist in the case of crystal 2. With care, twinning 
was to a large extent avoided in the remaining 
cases and the velocities may be taken to repre- 
sent those for relatively pure orientations. 
Table II presents the final results with average 
velocities computed using crystals 3-7 inclusive. 
The corresponding values of Young’s modulus 
Y=pc? are also included. These results are con- 
sistent with the known lattice structure of 
bismuth in the two orientations in question. 
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An Improved Method for Obtaining Photographs of Ripple Wave Actions 


M. E. RAQUET AND F. R. WATSON 
University of Illinois, Urbana, Illinois 
(Received February 21, 1938) 





























HE use of ripple waves to show reflection, 
diffraction and other actions of waves is 
well known.! Light shining upward through the 
glass bottom of a tank containing a shallow layer 
of water will be modified by the waves on the 
water surface and will show a pattern of the 
waves when striking a frosted glass surface, of 
which photographs may be taken. 
Figure 1 pictures the arrangement of the 
apparatus. Puffs of air directed perpendicularly 





Fic. 1. Apparatus for photographing ripple waves. 
1 Vincent, Phil. Mag. 43, 411 (1897); 45, 191 (1898). 
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against the surface of the shallow water in the 
tank set up waves that proceed outward in 
concentric circles. These are made visible strobo- 
scopically by flashes of light that pass up through 
the glass bottom of the tank and cast a shadow 
of the waves on the frosted glass above the water. 
The puffs of air and the flashes of light are regu- 
lated by the number of holes in the rotating 
wheel. If both pass through the same ring of 
holes, the waves appear to stand still; that is, 
the waves are shown in one position by the flash 
of light and then all is dark until each wave 
moves forward to the position of its predecessor 
when another flash illuminates them. If the 
flashes of light come less frequently than the 
waves by passing the light through another ring 
with fewer holes, the waves appear to move 
forward slowly. 

Figures 2 and 3 show the distorting effect of 
surface tension on the pattern of waves. In Fig. 2, 
two obstacles placed on the bottom of the tank 
to reflect waves project slightly above the water 
surface so that the water meniscus formed by 
the surface tension at the sides of the obstacles 
makes a sort of concave lens that distorts the 
direct light, while in Fig. 3, the water level is 
slightly above the top of the obstacles with a 
rounded edge of water, since the oily top of the 
obstacle prevents the water from covering it. 
Figs. 4 and 5 show clear outlines of the obstacles 
and no distortion of the waves; these being 
obtained by the method described later. Fig. 6 
shows a photograph of waves in a miniature 





Fics. 2 (left) ANp 3 (right). Photographs showing distortion of water waves 
by surface tension. Fig. 2, Water too shallow. Fig. 3, Water too deep. 
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Fics. 6 (left) AND 7 (right). Photographs of auditorium models. 


Fig. 6, distortion due to surface tension. Fig. 7, surface tension 
distortion eliminated. 





Fics. 4 (top) AND 5 (bottom). Photographs of water 
waves without distortion. Fig. 4, megaphone with short 
waves. Fig. 5, diffraction with long waves. 


model of an auditorium. The irregular outline of 
the model is due to the surface tension action 
which draws the water into irregular drops. 
Water waves reflected from this irregular surface 
are distorted. Fig. 7 shows the same model with 
the surface tension effects overcome, with a 
clear outline of the model and a regular undis- 
torted reflection of the waves. 

The method followed in eliminating the surface 
tension distortion was as follows. The reflecting 
obstacles were made plane parallel. These were 


placed on the plate glass bottom of the tank 
which had been carefully leveled ; then the water 
which completely covered the obstacles was 
slowly siphoned off until the layer of water over 
the obstacle was so thin that the waves could 
not pass but were completely reflected. (This 
suggestion was made to F. R. Watson by F. V. 
Hunt.) A series of photographs was made while 
the water was being lowered to secure the best 
condition,—too much water or too little water 
gave distorted photographs. It was necessary to 
clean the tank of all traces of oil, to use dis- 
tilled water and to make the obstacles scrupu- 
lously clean. This latter condition was obtained 
by the application of a special cleaning solution 
of acid-alcohol (95 percent alcohol, 5 percent 
glacial acetic acid) using a clean cloth and not 
touching the objects with the fingers, which 
always seem to have some oil traces on them. 
This procedure was followed by washing the 
obstacles with plain alcohol and wiping them 
with a clean cloth, after which they were placed 
in the tank, covered with water, and photographs 
taken as described. 











APRIL, 1938 J. A. 


Letter to 


A “Foaming Tone” Obtained from Vibrating 
Strings and Rods 


In a contribution to the editor, the author states that a 
novel “foaming tone” may be obtained by two types of 
instrument, and suggests several different arrangements 
for generating the tones. 

In one arrangement, a sonometer box is used that is 80 
cm long, 30 cm wide, and is equipped with 154 strings of 
metal or silk with a space between adjacent strings of 
about 2 mm. When an operator scrapes a finger nail or 
small metal object across the strings, a foaming-like tone 
should be obtained, depending on the velocity of scraping. 
The strings can either be tuned in unison or adjusted to 
give 22 octaves, each octave being the same as every other 
octave and covering the same frequencies. As a substi- 
tute for the finger nail, a stream of air from a tube may be 
used. 

In another type of construction, a row of teeth in a 
copper bar is drawn across a single string to generate the 
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the Editor 


tones. The arrangement calls for 350 teeth in a copper bar, 
1 meter long, 3 cm X3 cm in cross section, the teeth being 
confined in a space of 70 cm of the bar; the string being 
80 cm long and is to be stroked at right angles at a point 
8 cm from one end. For ease of operation, the bar should 
be equipped with a handle and a supporting guide. The 
teeth can be tuned in unison or in octaves. Various tones 
will depend on the velocity of scraping and on the segments 
of the string which are adjusted by hand. 

The author states that he has obtained the tones by 
apparatus roughly similar to the suggested instruments; 
but adds that he is not acquainted with the musical 
principles, and that the details of arrangement await the 
advice of musicians. He ‘“‘hopes that the problem will be 
investigated and completely developed.” 

Tar CuHuNnG-LING 
(formerly named Tai Chen) 
Soochow University, 


Soochow, China, 
October 15, 1937, 


Acoustical Society News 


HE Spring meeting of the Acoustical Society is to be 

held in Washington, D. C. on May 2-3, 1938. Those 
who attended the last Washington meeting will remember 
what an excellent program we had and will be glad to 
know that Mr. V. L. Chrisler is again in charge and is 
already laying plans for the Spring meeting. Address all 
correspondence pertaining to the meeting to Mr. V. L. 
Chrisler, Department of Commerce, National Bureau of 
Standards, Washington, D. C. The present plans call for 
a symposium on acoustic standards. On Monday evening a 
trip is planned through the Department of Justice Building, 





where we will hear about the ‘‘“G-Men”’ and their work. 
The Wardman Park Hotel will be headquarters. 

On March 14-15, a symposium on “Physics in the 
Automotive Industry,” arranged by the American Insti- 
tute of Physics, was held in Ann Arbor, Michigan, the 
department of physics of the University of Michigan 
acting as host. 

The fall meeting of the Acoustical Society will be held 
in November at Harvard University, Cambridge, Massa- 
chusetts. Professor F. A. Saunders is Chairman of the 
Program and Arrangements Committee. 
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Current Publications on Acoustics 


F. A, FIRESTONE 


147 East Physics Building, University of Michigan, Ann Arbor, Michigan 
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ain Reviews of Contemporary Papers 

int EMBERS who have occasion to translate important papers from foreign journals will be 

ane doi 1 service to th bers who d d i ily, 4 

The : oing a real service to those members who do not rea the foreign languages easily, if they 

pind will prepare a review of such articles for publication in this section. 

ins These reviews should be of the nature of a lengthy abstract (500 to 1500 words) rather than a 
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i author's contribution as possible. One or two figures may be included if desired. 
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On the Vibration Insulating by Means of Springs and Putting 
Damping Materials. E. Meyer, L. Keme.. Zeits. f. 
en) techn. Physik 18, 299-304 (1937).—To insulate vibrating wo=1/(MF)}, a =wo/w, B=R/(M/F)} 
apparatus, springs, layers of rubber or cork, etc., are used. 


Fig. 1 shows the equivalent diagram and the electrical this can be written 


U=10 logio (a26?+ a") /[a*6?+(1—a?)?]. db. 


When a=wo/w=1, the transmission curve shows a reso- 
nance peak, and for w/wo=1/a=v2, k=K, everything is 
the same as if the system would not be insulated. For 
Foc. 1. Macheatest end lectitenl eusteay of ws greater values of w/wo, the transmission decreases at a rate 
insulated vibrating system. which is the greater the smaller R. For most of the ma- 
terials used in practice, R is not a fluid resistance, but is 
due to hysteresis in the material and consequently decreases 
as 1/w approx. Then a8=7 is a constant, and for w/wo> v2 | 
the transmission becomes practically independent of R, 
decreasing 12 db for the frequency twice as great. Fig. 2 





ha ala 


analogy of an insulated vibrating system. “k’” is the 
periodic force producing the vibration, M the mass of the 
system, F its compliance, R the dissipation resistance and 
K the force exerted on the ground. Using a logarithmic 











/ scale, the transmission U will be given by : : ; 
represents the experimental arrangement. A moving coil 
. K? R+1/°F in the air gap of a pot shaped magnet serves as generator. 
U=10 logio ik =10 logio R?+(wM—1 loF)2 The force exerted on the base plate is registered by means 
rk, of a quartz crystal and a logarithmic voltmeter. The 
mechanical parts have been cemented together to ensure 
the close contact, and great care had to be taken to avoid 
sti- resonances. Leaving out the test piece (short circuit test) 
the the force exerted on the quartz crystal should be inde- 
yan ! pendent of the frequency, which gives a means of con- 
| trolling the working of the apparatus. 
eld 
sa- | 
the db 


VW 


Fic. 2. Experimental arrangement. A, suspension; B, pot shaped 50 100 200 500 1000 2000 /sec. 
magnet; M, moving coil, vibrating mass (aluminum or iron); c. 
amplifier cog valve voltmeter; D, recording instrument; E, case plate; Se ; 
G, quartz crystal; H, heterodyne oscillator; J, test piece. Fic. 3. Transmission curve for a spring insulating. 
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Figure 3 gives a record obtained with a screw-spring 
insulating. The great number of resonances greatly de- 
creases the vibration insulation. They are still present, if 
the fundamental period of the system is lowered by means 
of a heavier mass. Only a specially constructed annular 
spring showed the exact insulation as demanded by the 
theoretical curve. Also layers of cork and rubber have 
natural periods of their own, which will fall into the 
working range if the layer is thick. The lowest natural 
frequency of such a system is proportional to the thickness 
of the layer, whilst the first resonance of the layer occurs, 
when its thickness is equal to \/2, \ being the wave-length 
in the particular material. The transmission curve for 
hair felt shows a very damped resonance peak, and the 
insulation is not good, even for high frequencies. Ap- 
parently the dissipation mechanism of felt is different from 
the usual hysteresis loss. Fiber board shows a nonlinear 
behavior but its insulation is good (nonlinear resonance 
curve). A mechanical low pass filter consisting of four 
sections gave a very great decrease of the transmission for 
the higher frequencies. Each section was made up of a 
mass (M=600 grams) and a compliance (rubber disk, 
4 ohms dia., thickness 10 mm). 

In conclusion, it can be said that above the fundamental 
resonance, springs or layers of rubber or cork insulate 
equally well (12 db increase in the insulation for a frequency 
twice as great), but if the material has any additional 
resonances within the frequency range under consideration 
the transmission is greatly increased. 

E. J. SKUDRZYK 


The Theoretical Derivations of the Laws of Reverbera- 
tion. H. CREMER, L. CREMER. Akustische Zeits. 2, 225-241, 
296-302 (1937).—This paper gives a critical survey of the 
various methods of deriving the reverberation formulae, 
pointing out the approximations and the limits of validity 
in every case. 

(1) Te obtain Sabine’s law on statistical considerations, 
it is necessary to assume that every ray of sound crosses 
all points of space, and that it crosses every one of them 
in all directions. This last assumption is not necessary if 
the so called cosine law holds for the absorber: 


a'(3) = —log [1—a(d) ]=constant/cos #, 


where a(#) is the absorption factor of the wall and # the 
angle of incidence. a! is generally called the absorption 
exponent. Let us prove this for a rectangular chamber 
the edges of which are J;, 12, 13. Let cos 31, cos 32, cos 33 
be the direction cosines of the sound ray with regard to 
the three faces. (The same ray always strikes the same 
face under the same angle of incidence.) Then the number 
of times the ray will strike the threé walls during one 
second will be given by 
c cos 3;/21;, 


C COS 32/22, c cos 33/2]s3, 


respectively. Herice we obtain for the energy at the time ¢ 
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T= To(1 —d0)° cos 31 t/2y (4 —4a1)° cos By t/2l; 
++ «(1 —asi)¢ ©08 9stls, (1) 


where do, @11-++@39, @3; are the absorption factors of the 
six walls. This can be written 


T:=To:-%* with 26=(c/4V)22a,(8,) cosd, Fn». (2) 


n 
Different rays will thus be damped differently, according 
to their angle of incidence. Only if (1) holds 6 becomes 
independent of # and Sabine’s formula holds exactly. 
Most of the absorbers show an increase of absorption with 
the angle of incidence so that the cosine law frequently 
will be a fair approximation. 

(2) If our calculations are based on the wave equation, 
reverberation is identified with the dying out of the 
damped modes of vibration. Let us first consider the 
natural functions of an undamped rectangular chamber. 
They are 


cos kx cos key cos k;32/'*,* (Ry=nyt/ly). 


Every one of these can be decomposed into eight waves 
of the form 
Jib hyrt hoy kg zt) 


The resonance frequencies are given by 
f= (c/2) [(m1/21;)?+ (2/21)?+ (n3/213)? ]}. 


They can be obtained graphically as the distance from the 
origin to a point 1, #2, 23 of a rectangular network, the 
units of which are ¢c/2l;, c/2lo, c/2l3 in the x, y, and z 
direction, respectively. The direction of the resonance in 
space is the same as that of the ray connecting the origin 
to the corresponding point of the network. If the frequency 
is large, the number of resonances between the frequencies 
fand f+Af asymptotically is equal to the volume contained 
in an octant of the spherical shell with the radii f and 
f+Af, divided by (c/2)3(1/1:)(1/l2)(1/l3), the volume of 
one cell of the network. Similarly the number of resonances 
between f and f+Af in the direction AQ is asymptotically 
given by 

= ——. = SV eafan 
(c/2)§(1/h)(1/l2)(1/ls) 8 


and this is independent of the angle 2, no particular 
direction being favored in the limit. Thus every mode of 
vibration leads to a uniform distribution of energy and 
the sound field is diffuse, if within a certain frequency 
range the modes are excited with equal amplitudes. These 
considerations are not limited to a rectangular chamber, 
but as can be shown apply to any shape of space. 

(3) M. J. Strutt introduces a resistive force r(xyz)v, “‘v” 
being the particle velocity. Using Rayleigh’s method for 
small perturbations, his results are obtained as follows: 
The energy dissipated each second in the whole space is 


An 


N= JSr(x, y, 2)v*(xyz)dV 
and the total energy will be given by 


EV= JS pr*(x, y, 2)dV, 
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where E is the average energy density. Putting 


vy? = 792/268 
N=-—(d/dt)(EV) =26EV. 
26= frvdV/p fvdV. 


we get 
Hence 


If the energy density is uniform and r(xyz) spread over a 
number of wave-lengths, v? can be canceled out in numer- 
ator and denominator, giving Strutt’s formula 


256= SfrdV/pV. : (3) 


For porous absorbers of thickness d, and area F,, we get, 
performing the integration, 





26= — a | = 


It follows from above that d, should be a few wave-lengths 
but in no case less than \/4. Also, the energy calculation 
was based on the undisturbed vibration which is only 
allowable if (r/pw)<K1. Under these conditions (3) may 
even be used to calculate the frequency variation of the 
absorption of porous absorbents. There is no difficulty in 
showing that the cosine law 


a'(3) =2rd/pce cos 3 


applies to such absorbers. 

(4) Schuster and Waetzmann introduce a complex wall 
resistance and Van den Dungen derived the reverberation 
formula for great values of wall impedance, with the aid 
of a perturbation method. The exponential law is found 
to hold, if the wall resistance is independent of the angle 
of incidence. This takes the place of the previously men- 
tioned cosine law, and always holds for porous absorbents 
of large streaming resistance at the low frequencies. For 
other values of streaming resistance and frequencies 
deviations from Sabine’s law are to be expected. 

(5) The rectangular chamber lends itself to calculation, 
even if any wall resistance is admitted. The general 
solution of the wave equation can again be split up into 
eight waves of the type 


J (iky—8y zt (ikg—sq)yt (ikg—sg)2+(io—6)E, 


By a few calculations we get (2) again. The only assump- 
tion we require is that 


(s2+s2+53)IKX(k2+k2 +3?) i, 


which is true in all practical cases. Previously we have 
derived Eq. (2) for very short wave-lengths; we now find 
that this equation is still valid even if the wave-length is 
of the same order of magnitude as the dimensions of the 
chamber. 

The authors also treat the case of a prismatic chamber 
the cross section of which is a right angled bilateral 
triangle. Such a chamber results if one of square section 
is divided into two halves along the diagonal. The partition 
wall has no effect on the sound field, as the particle velocity 
perpendicular to it is zero (the two resonances of the 
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square chamber which have the same frequency but 
different directions being excited equally). Hence such a 
triangular chamber has the same natural frequencies as 
the square chamber. The triangular chamber is of practical 
importance, as it is free from flutter echoes. 

(6) The rectangular chamber could be used to determine 
the absorption factor at low frequencies and various 
angles of incidence. This would be done by exciting one 
mode at the time, measuring the reverberation time and 
calculating the corresponding angle of incidence. 

E, J. SKUDRZYK 


Construction and Investigation of a Damped Room. 
W. JANovsky, F. SPpANDGcK. Akustische Zeits. 2, 322-331 
(1937).—The cross section of the damped room is irregular, 
approaching the shape of a triangle. Its sides are about 
8 m long and the distance floor to ceiling is 3 m. The walls 
are lined with Heraklith plates, then comes a 10 cm thick 
layer of slag wool, a 10 cm layer of cotton wool and 
finally a 30 cm layer of cotton wadding suspended edge- 
wise. The slag wool is held in position by means of a wire 
grating which at the same time acts as electrical screen. 
Perpendicular to the wall, strips of wire gauze are fixed at 
distances of one meter. Their object is to prevent the 
rapid spread of fire. They are supposed to act similar to 
Davy’s safety lamp. 

A measurement of the propagation along the axis of 
the room gave the same results as have already been 
described by E. H. Bedell in this journal.! The free air 
propagation could be obtained up to a distance of 2} m, 
the point source being placed in the middle of the room. 
For greater distances the pressure decreased at a greater 
rate than 1/r. E. H. Bedell suggested that along the 
damped surface the wave velocity is reduced, and the 





Fic. 1. Difference between the acoustical and optical angle of incidence 
plotted against the height of the loudspeaker from the floor. 


wave front bent or diffracted towards the absorbing 
surface, more energy being absorbed than if the wave 
remained spherical. To prove this, the height of a loud- 
speaker was varied from 15 to 150 cm whilst the angle of 
incidence at which the sound was striking the floor was 
determined, at a distance of 4 m from the speaker. This 
was done by means of two microphones one close to the 
floor, the other movable in a plane 30 cm from the ground. 
The output of the two separate microphone amplifiers 
was coupled together in antiphase and the movable 
microphone adjusted to the minimum of loudness. This 
gave the planes of equal phase and hence the direction of 
the incident sound. The result is represented in Fig. 1. 
A quantitative idea of the equivalent damping due to the 
bending of the waves could be obtained by measuring the 
intensity of sound whilst the height of sound source and 
receiver were varied. The result depends to a great extent 
on the directional characteristics of the sound source. On 
the other hand a loudspeaker of uniform direction char- 
acteristics showed marked directivity when placed close 
to the damped wall. 

Finally the absorption coefficient of the damped wall 
was measured. The loudspeaker was placed at the floor 
and the microphone suspended at a mean height. A short 
electrical impulse was sent through the loudspeaker coil, 
and the microphone output, corresponding to the direct 
impulse and its first reflection, was oscillographed. 14 
percent of the pressure amplitude was reflected, giving an 
absorption coefficient of 98 percent. 

E. J. SkuDRzYK 


1E. H. Bedell, J. Acous. Soc. Am. 8, 118-125 (1936). 


Magnetic Strip Recording with Tapes and Circular Re- 
cording and Reproducing Heads. H. Luseck, Akustische 
Zeits. 2, 273-295 (1937).—The object of this work was to 
investigate the circular recording and reproducing heads 
in conjunction with nonmagnetic tapes covered with 
ferromagnetic powder. 

The circular recording head (recording horseshoe mag- 
net) consists of a ferromagnetic ring with a very narrow 
air gap (0.03 mm), the ring carrying a coil. The ring is 
usually made from isoperm, but it is shown in this publi- 
cation that much better results can be obtained with 
substances that have a very large permeability, remanence 
being reduced to a negligible value by means of an addi- 
tional air gap. The flux leaking out from the air gap, and 
only that portion, magnetizes the tape which is drawn in 
front of it. Such an arrangement leads to a very sharp 
definition of the magnetic field in front of the air gap, and 
with quite simple means tape velocities as low as 60 
cm/sec. still give a linear reproduction of the frequency 
range 50-10,000 c/sec. Nonmagnetic tapes covered with 
ferromagnetic powders are of advantage as they are 
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cheaper and lighter than steel wire and at a junction 
they are simply pasted together. Even after playing a tape 
10,000 times wearing could hardly be noticed. 

The recorded sound is erased by magnetizing the tape 
to about its saturation value. By means of an additional 
d.c. current in the recording head, when the next recording 
takes place, the new working point is brought back to the 
straight portion of the magnetization characteristics. 

A consideration of the magnetic intensity following one 
particle of the tape leads to an expression for the highest 
frequency recorded without great distortion. This fre- 
quency is proportional to the tape velocity and inversely 
proportional to the a.c. amplitude of the recording current, 
and is the greater the sharper the magnetic definition of 
the air gap. It is therefore necessary to reduce the a.c. 
amplitudes when recording high frequencies. These con- 
clusions are confirmed by measurements. A consideration 
of the erasing and recording process leads to an analysis 
of the most favorable working point and shows that the 
minimum of distortion corresponds to a maximum of 
intensity on the record. It is found that the hysteresis 
loop of the tape is a direct characteristic of the recording 
process. 

The receiving head is similar to the recording head. 
For an understanding of the receiving process the tape is 
best considered as a series of elementary magnets. Those 
facing the air gap send a magnetic flux through the 
induction coil, those opposite the solid portion of the ring 
are short circuited by the latter. A short calculation gives 
an approximate expression for the e.m.f. induced in the 
receiving head and extensive measurements show how to 
correct it. It is found that the induced e.m.f. essentially 
is determined by the wave-length recorded on the tape 
and the tape velocity, also that there is a certain amount 
of demagnetization taking place, increasing exponentially 
with the frequency, at a rate which is the greater the higher 
the permeability of the tape. For small slits the frequency 
response becomes independent of the length of the air gap. 

The ferromagnetic powder on the tape is not quite 
uniform, therefore, there is a certain amount of background 
noise, but it is hoped, that improvement in this respect 
soon will be obtained. 


E, J. SkUDRZYK 


CORRECTION 


Disturbance Produced by Nonlinear Distortion. H. J. 
V. BRAUNMUHL AND W. WEBER. Akustische Zeits. 2, 
135-147 (1937). In the review of this paper published in 
this journal (9, 170 (1937)) it was stated that the sounds 
were picked up from high quality records, whereas the 
authors inform us that the sounds were picked up directly 
from the studios. 
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